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ABSTRACT 


In  this  work  an  investigation  is  made  of  uniform  approximations 
to  the  solutions  of  linear  differential  equations  with  variable  coefficients. 

The  ordinary  differential  equations  are  replaced  by  an  appropriate  set  of 
partial  differential  equations  that  determine  the  unknown  function  in  terms 
of  a  set  of  independent  "time  scales.  "  The  time  scales  are  determined  so 
as  to  obtain  uniformly  valid  approximations.  The  partial  differential 
equations,  in  conjunction  with  the  requirement  of  uniformity  of  the  approximation 
in  a  given  interval,  determine  the  time  scales  through  a  Bet  of  "clock  functions" 
k.  ,  which  may  depend  on  the  interval  of  interest.  It  is  essential  for  the 
success  of  the  approximation  that  the  clock  functions  be  nonlinear  functions  of 
time,  in  addition  to  being  complex  quantities.  The  constant  coefficient  case 
arises  as  a  natural  limit.  Thus  the  present  approach  generalizes  earlier  time 
scale  analyses.  With  this  generalization  we  recover  for  second  order  systems 
the  Liouville -Green  (or  WKBJ)  approximation.  The  difference  between  the 
present  approach  and  the  PLK  method  irj  emphasized  with  examples. 

Bounds  on  the  errors  committed  are  established  for  the  second  and 
third  order  equations.  The  use  of  two  time  scales  (with  nonlinear  clocks) 
enables  us  to  obtain  approximations  to  the  amplitude  and  phase  of  each 
of  the  modes  of  n**1  order  equations. 

The  prototypes  that  are  of  interest  are  the  linearized  equations 
governing  the  motion  of  VTOL.  aircraft.  These  equations  constitute  a 
system  of  rather  high  order  in  the  time  derivative  (third  or  fourth  order  for 
motion  in  the  plane  of  symmetry).  The  approximation  method  obtains  the 
aircraft  variables  in  terms  of  simply  calculable  functions  of  the  stability 
derivatives.  The  frozen  analysis  of  the  aircraft  equations  suggests  solutions 
of  the  simple  form 


Z 

i 


A.  e 

x 


Aft 


with  A.  and  X.  slowly  varying  in  time.  We  introduce  new  independent 
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variables  (time  scales)  and  T  ^  to  represent  the  amplitude  and  phase 

of  the  modes  and  express  the  aircraft  motion  in  the  form 

E  a.  (T  .  )  e  Tli 
i  1  01 

and  being  determined  appropriately.  The  results  of  the  approximation 

are  compared  with  numerical  integrals  of  the  aircraft  equations  for  the 
third  order  hovering  system  and  the  complete  fourth  order  equations  which 
allow  for  the  transition  from  hover  to  forward  flight.  The  approximation 
is  found  to  be  very  accurate  (to  within  10%  error)  for  the  third  order  system. 
For  the  fourth  order  system  comparable  accuracy  is  obtained  except  near 
the  transition  point.  However,  qualitative  features  of  the  exact  solution  are 
not  lost.  A  uniform  description  of  the  aircraft  motion  is  thus  obtained. 
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CHAPTER  I 


INTRODUCTION 

1.  1  Origin  of  the  Problem 

This  dissertation  concerns  the  analysis  of  systems  described  by 
linear  differential  equations  with  variable  coefficients.  This  is  a 
classical  problem,  one  that  has  attracted  and  occupied  the  interest  of 
mathematicians  and  physicists  since  about  the  middle  of  the  seventeenth 
century.  The  range  of  interest  is  vast,  from  esoteric  fields  of  study  such 
as  topology  and  the  qualitative  theory  of  differential  equations  to  the  very 
practical  task  of  analyzing  actual  systems  arising  in  modern  physics  and 
engineering.  The  fields  of  applied  mathematics  and  physics  are  replete 
with  examples  of  linear  differential  equations  with  variable  coefficients 
embracing  such  diverse  disciplines  as  celestial,  quantum  and  classical 
mechanics,  wave  motion  in  inhomogeneous  media,  rocket  flight  through  the 
atmosphere,  and  so  on. 

The  importance  of  the  study  of  such  systems  has  increased  in  recent 
years  mainly  because  of  three  factors --the  attempt  to  comprehend  the  more 
subtle  phenomena  in  nature,  and  the  advent  of  sophisticated  dynamical 
systems  of  modern  engineering,  and  the  need  to  stabilize  and  control  such 
systems.  In  the  aerospace  sciences,  such  problems  arise  for  example 
in  the  analysis  of  the  dynamics  of  motion  of  a  vertical  take  off  and  landing 
(VTOL)  type  of  aircraft  from  hover  to  forward  flight.  The  aerodynamic 
parameters,  since  they  depend  on  the  flight  condition,  change  continuously 
during  the  transition  and  the  differential  equations  describing  the  motion 
have  nonconstant  coefficients.  Another  example  is  the  motion  of  a  space 
vehicle  negotiating  a  flight  from  or  reentry  into  the  earth's  atmosphere. 

The  variation  of  density  with  altitude  gives  rise  to  variable  coefficients  in 
the  differential  equations. 

Another  way  in  which  such  equations  arise  is  in  the  theory  of 
partial  differential  equations.  In  problems  which  admit  of  the  separation 
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of  variables,  the  ordinary  differential  equations  thus  obtained  will  often 
contain  variable  coefficients.  A  well-known  example  is  that  of  wave  motion 
i  i  inhomogeneous  media.  The  wave  equation  yields,  after  separating  the 
variables,  a  function  periodic  in  time  and  a  space  dependent  function  which 
satisfie  s  an  ordinary  differential  equation  with  variable  coefficients. 

A  similar  problem  also  arises  in  dealing  with  nonlinear  differential 
equations.  The  variational  equations  corresponding  to  a  particular  known 
solution  of  a  nonlinear  differential  equation  often  turn  out  to  be  linear 
equations  with  variable  coefficients.  The  variational  equations  may  be 
useful  in  deducing  stability  information. 

1.2  Historical  Sketch  and  Review  of  Existing  Work 

In  attempting  to  solve  linear  time -varying  equations,  the  analyst  is 
beset  with  considerable  difficulty.  The  history  of  formal  methods  of 
integration  practically  ends  in  the  latter  half  of  the  eighteenth  century. 

The  first  order  linear  differential  equation  (l.d.e.  )  can  be  solved 
exactly  by  means  of  a  quadrature.  The  second  order  l.d.e.  with  arbitrary 
coefficients,  however,  is  quite  another  story.  It  can  be  shown  (Ref.  1), 
though  not  without  difficulty,  that  it  cannot  be  solved  in  general  by  a  finite 
number  of  quadratures  and  elementary  operations.  For  particular 
variations  of  the  coefficients,  however,  the  various  standard  transcendental 
differential  equations  are  obtained.  For  example,  in  the  l.d.e. 

y"  +  uij(  t)  y'  +  uQ(t)  y  =  0  (1.2.1) 

(The  primes  denote  differentiation  with  respect  to  the  independent  variable  t) 
1  n3 

if  Ui  -  -  and  uj  Q  =  1  - -p-  ,  we  obtain  Bessel's  equation  cf  order  n; 
if  u)  i  =  0,  u  =  -t,  the  Airy  equation;  if  Uj  =  0,  u>  =  (a  -  2qcos2t) 

(a,  q  constant)  we  have  the  familiar  Mathieu  equation,  and  so  on.  These 
are  well  known  equations  and  exact  solutions  ire  available. 

In  all  but  the  standard  cases  exact  solutions  are  net  known  and  the 
explicit  forms  of  the  asymptotic  solutions  cannot  easily  be  written  down. 
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One  has  to  resort  to  approximations  under  these  circumstances,  and  the 
literature  on  this  subject  is  very  rich  indeed.  Many  authors  have  contributed 
to  the  state  of  the  art  and  the  theory  is  quite  extensive. 

The  idea  basic  to  most  of  the  schemes  for  approximation  is  to  show 
that  under  certain  conditions  the  given  equation  may  be  represented  by 
another  simpler  equation  which  can  be  solved,  such  that  the  difference 
or  the  error  between  the  two  solutions  is  small.  The  most  direct  approach 
is  that  of  perturbation  theory  introduced  by  Poincare  (Ref.  2)  and  this  can 
be  applied  when  there  is  a  small  parameter  which  exists  in  the  physical 
system.  Then  direct  perturbation  theory  yields  approximations  of  better 
and  better  accuracy.  There  are  many  cases  where  such  a  representation 
becomes  unsatisfactory  over  part  of  the  domain  of  interest.  This 
phenomenon  is  termed  nonuniformity  and  will  be  examined  in  detail  later. 

More  sophisticated  perturbation  methods  have  been  devised  so  as  to 
overcome  this  difficulty  (but  sometimes  only  partially). 

The  available  methods  of  analysis  and  their  application  to  time -varying 
control  systems  have  been  examined  in  Ref.  3.  There  the  methods  are 
categorized  into  classical  theory  of  differential  equations,  matrix  methods, 
methods  of  integral  transforms,  etc.  App_oximations  such  as  the  method 
of  collocation,  Schellkunoff's  wave  perturbation  method,  etc.,  have  been 
discussed  in  detail  and  examples  given. 

Among  the  many  approximations  existing  in  the  literature,  one  which 
has  enjoyed  considerable  fame  is  the  so-called  Liouville -Green  (or  WKBJ) 
approximation,  sometimes  named  with  other  permutations  of  the  letters.  This 
name  refers  to  the  representation  of  the  solution  of  the  differential  equation 

y"  +  ha  u;  (t )  y  =  0  ;  u>  >  0  (1.2.2) 

in  the  form 

y  =  AuT1^4  exp  (  ih  juu1/-  dt  )  +  exp  ( -ih  juu1^8  dt)  (1.2.3) 

where  A  and  B  are  arbitrary  constants.  The  usefulness  of  such  an 
approximation  will  be  examined  later  and  the  conditions  of  validity  established. 
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The  naming  of  this  scheme  of  approximation  has  had  a  diverse  history. 
Apparently  the  use  of  such  approximate  solutions  may  be  traced  to  Carlini  (1817) 
(Ref.  4)  who  considered  a  specific  equation  of  the  Bessel  type.  Liouville  (Ref.  5) 
and  Green  (Ref.  6)  (1837)  derived  the  approximation  for  more  general  equations, 
although  their  derivations  lacked  rigor  and  were  only  valid  in  restricted  regions 
of  the  complex  plane.  More  recently  the  method  has  been  referred  to  by 
physicists  as  the  WKB  method,  after  Wentzel  (1926)  (Ref.  7),  Kramers  (1926) 

(Ref.  8),  and  Brillouin  (1926)  (Ref.  9),  though  the  letter  J  is  often  added  to 
acknowledge  the  contribution  of  H.  Jeffreys  ( 1 923  )(Ref J  0)  in  connecting  the 
approximate  solutions  valid  on  either  side  of  a  transition  point--a  point  on  either 
side  of  wMc’.,  u.(i)  has  opposite  signs;  and  H.  Jeffreys  (Ref.  11)  has  recently 
pointed  out  that  he  himself  had  been  anticipated  by  Gans  (Ref.  12).  Just  to 
complete  the  overwhelming  list  of  the  different  names,  Bailey  (Ref.  13)  has 
chosen  to  call  it  the  "L.  R.  approximation",  after  Liouville  (Ref.  5 )( 1  837)  and 
Rayleigh  (Ref.  14)  (1912).  In  a  recent  review  article  by  B.S.  Jeffreys  (Ref.  15), 
the  name  asymptotic  approximation  method  is  suggested.  Perhaps  among  all 
these  various  names,  that  of  "phase -integral  method"  as  used  by  Heading  (Ref.  16) 
seems  to  be  most  appropriate  as  it  does  not  refer  either  to  its  discoverer  or  to 
its  method  of  derivation.  Though  the  name  WKBJ  seems  to  be  widely  prevalent, 
we  sha.l  follow  Olver  (Ref.  28)  ard  use  the  name  "L-G  approximation"  after 
Liouville  and  Green  who  derived  the  approximation  first. 

It  must  be  noted  in  passing  that  it  is  possible  to  raise  some  criticisms 
about  the  point  of  view  adopted  in  the  above  method  of  approximation. 

However,  approximations  quite  often  serve  a  useful  purpose  in  mathematical 
physics  and  engineering  and  have  led  to  illuminating  results  in  many  cases. 

When  valid  the  above  approximation  scheme  shows,  for  example,  the 
connection  between  classical  and  quantum  mechanics  by  providing  approximate 
solutions  to  Schro  dinger 's  equation.  In  response  to  such  criticisms  the 
author  cannot  do  better  than  to  fall  back  on  Heading  (Ref.  17),  who  says, 

"For  example,  with  no  just  foundation  for  such  remarks,  Smyth  has 
criticized  a  paper  making  use  of  the  method  by  writingi'It  should  be 
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observed  that  the  authors  have  used  a  solution  which  is  a  very  poor 
approximation  to  the  given  problem  as  an  approximate  solution  to  another 
problem.  It  is  certainly  not  to  be  expected  that  the  results  obtained  in 
this  manner  will  have  any  connection  with  the  original  problem.  1 
Introducing  a  new  approach  that  leads  to  a  difference  in  applicability, 

Hines  has  observed  that  his  new  method  yields  an  approximate  evaluation 
of  the  exact  solution,  rather  than  an  exact  evaluation  of  an  approximate 
solution  as  is  found  in  the  WKBJ  method!"  Perhaps  the  wise  remarks  of 
Schellkunoff  (Ref.  18)  concerning  the  approximations  should  be  recalled: 
"There  is  something  in  human  nature  that  makes  one  yearn  for  the  exact 
answer  to  a  given  problem.  In  particular  it  makes  little  difference  whether 
a  given  problem  is  solved  approximately  or  replaced  by  an  approximating 
problem  which  is  then  solved  exactly.  " 

1.3  The  Liouville -Green  (  or  WKBJ)  Approximation 


This  approximation  can  be  derived  by  first  converting  (1.2.2) 
into  the  Ricatti  equation 


by  the  transformation 


1_ 

h 


(1.3.  1) 


where  a  small  parameter  X  has  been  introduced  in  (1.3.1  ).  The  equation 
(1.3  .1  )  can  be  formally  solved  by  expanding  z  as  follows 


z 


X 


t 


dt  r  x,  xk 

k=0  k 


(1.3.2) 


where  t  is  a  constant.  Substituting  this  into  (1.3  .1  )  and  equating  like 
powers  of  X,  we  get  a  set  of  equations 
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=  ±  i/u) 


x 

o 


dx 


k-1 


dt 


k 

-  E  x.x  .  k  =  1,  2,  3,  .  .  . 
j=0  J  k-j 


(1.3.3) 


from  which  the  series  £  x  X  can  be  determined.  This  is,  for  example, 

k=0  k 

the  method  used  by  H.  Jeffreys  (Ref.  10). 

It  is  known  that  this  series  is,  in  general,  not  convergent,  but  is 
only  asymptotic.  From  the  above  we  can  get  the  LG  approximation. 


y(t)  =  (  £a  )~^4  exp[  ±  i  J 


u> 


Va 


dt] 


(1.  3.4) 


if  we  neglect  the  other  terms.  Similar  derivations  of  this  approximation 
are  given  in  many  books  on  quantum  mechanics,  but  few  give  the  precise 
conditions  of  validity. 

The  parameter  X  was  introduced  only  as  a  formal  mathematical  tool 
for  obtaining  the  desired  expansion  for  y.  It  is  seen  that  bot'.  in  the  1.  d.e. 

(1.2.2)  and  the  LG  solution  (1.3.4),  the  parameter  X  and  the  function 

(1) 

“>(t)  appear  only  in  the  combination  ,  and  so  one  can  simply  write 
ua  rfm 

instead  of  ^  .  The  definiteness  of  the  sign  of  must  be  ensured  by 

branch  cuts  in  the  complex  t  plane  radiating  outwards  from  the  zeros 

of  <s,  and  the  integration  must  take  place  along  paths  not  crossing  these  cuts. 

The  above  treatment  is  inadequate  since  neither  estimates  of  the  errors 

of  the  approximation  nor  the  regions  of  validity  in  the  complex  plane  are 

given.  The  solutions  of  equation  (1.2.2)  are  single -valued  in  a  domain 

containing  no  singularities  of  *u(t);  but  because  of  the  fractional  powers  of 

<u(t)  the  LCt  solutions  are  not  single -valued  and  so  it  is  clear  that  the 

solutions  (1.3 .4)  are  valid  only  in  restricted  regions  of  the  complex  t  plane. 

One  of  the  purposes  of  this  work  is  to  demonstrate  a  more  general 

method  of  approximation  which  yields  the  LG  solution  as  a  special  case 


6 


as  well  as  approximate  solutions  valid  where  the  LG  solution  is  not. 
The  study  of  asymptotic  approximations  to  the  solutions  of  l.d.e. 


would  be  incomplete  without  at  least  a  brief  discussion  of  the  Stokes 
phenomenon.  This  name  is  given  to  the  discontinuous  changes  in  the 
arbitrary  constants  that  occur  in  the  asymptotic  solutions  of  certain 
differential  equations. 

It  is  rather  interesting  to  follow  the  beginnings  of  the  observation  of 
this  phenomenon.  In  a  letter  to  a  certain  young  lady,  Sir  George  Gabriel  Stokes 
wrote  in  1857  (Kef.  19)  (and  the  present  author  wishes  to  beg  the  indulgence 
of  the  reader),  "When  the  cat's  away  the  mice  are  at  play.  I  have  been 
doing  what  I  guess  you  won't  let  me  do  when  we  are  married;  sitting  up 
till  3  o'clock  in  the  morning  fighting  hard  against  a  mathematical  difficulty. 

Some  years  ago  I  attacked  an  integral  of  Airy's,  and  after  a  severe  trial 
reduced  it  to  a  readily  calculable  form.  But  there  was  one  difficulty  about 
it  which,  though  I  tried  till  I  almost  made  myself  ill,  I  could  not  get  over, 
and  at  last  I  had  to  give  it  up  and  profess  myself  unable  to  master  it. 

I  took  it  up  again  a  few  days  ago,  and  after  two  or  three  days  hard  fight, 

I  at  last  mastered  it.  " 

The  phenomenon  was  first  observed  in  connection  with  the  Airy 
equation: 

y"  -  zy  =  0 

For  small  |zj  ,  the  general  solution  comprising  of  two  independent  power 
series  solutions,  would  involve  two  fixed  arbitrary  constants.  It  was 
observed  that  if  for  a  certain  range  of  arg  z,  the  general  solution  was 
represented  by  a  certain  linear  combination  of  the  two  asymptotic  solutions, 
then  in  a  neighboring  range  of  arg  z  it  was  by  no  means  necessary  for  the 
same  linear  combination  to  represent  the  same  general  solution.  Stokes  (Ref.  20) 
in  fact  showed  that  the  arbitrary  constants  must  be  changed  discontinuously 
on  crossing  certain  lines  in  order  to  provide  an  asymptotic  representation 
of  a  continuous  function  for  both  ranges  of  arg  z.  These  lines  are  called 
anti -Stoke  a  lines. 
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The  Stokes  phenomenon  occurs  because  an  asymptotic  series  is  not 
unique.  For  example,  the  two  functions 

i  i  \  *(z)  "z 

and  ib  2  (z)  =  — —  +  e 


have  the  same  asymptotic  expansion  for  /  z  j  —  <*>  when  Re  f  z  ]  >  0.  Besides,  the 
asymptotic  expansion  for  (|)  a  (z )  will  change  drastically  from  Re  r  z]  >0  to 
Re,  z'  <  0.  The  asymptotic  form  of  ib  a  (z)  as  arg  z  change s  would  reveal  a 
discontinuity  at  arg  z  =  ^  ^  ,  etc.  These  discontinuities  are  only  apparent, 

however,  and  are  essentially  a  result  of  the  nonuniqueness  of  asymptotic 
expansions  (Ref.  21). 

One  may  also  consider  that  the  Stokes  phenomenon  occurs  because  the 
operations  of  analytic  continuation  and  taking  the  asymptotic  expansion  do  not 
commute.  Tn  other  words,  let  f(t)  have  the  asymptotic  expansion: 
f(t)  =  fQ(t)  +  efj  (t)  +  ... 

Let  f.(t)  be  now  analytically  continued  into  f.(z).  On  the  other  hand,  let  the  function 
i  l 

f(t)  be  analytically  continued  into  f(z).  If  we  now  obtain  the  asymptotic  expansion 
f(z)  =  f^\z)  +  e  f^(z)  +  ... 
we  will  find  that  in  general: 

f.(z)  t  f(0(z) 

This  is  the  Stokes  phenomenon. 

The  Stokes  phenomenon  in  LG  theory  arises  as  follows  (Ref.  22). 


The  approximations  to  the  solutions  of 


y"  +  id  (t)y  =  0 


are  given  by 


l 

y+  (t;  y ,  6  )  =  6  u'  "l/<t  (t)  cos(  j  d,1^3  dt  +  y  ) 


when  u  is  positive  and 


i 

y_  (t;  o,/3  )  =  (-u)(t)  )~lt*  f  a  exp  (  j  dt) 

a 

t  i 

+  ft  exp(  -  J  ( -u)  r  dt  ) 


(1.3.5) 


(1.3.6) 


(1.3.7) 


when  ir  (t)  is  negative.  a,(3,y,f>  are  arbitrary  constants. 
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When  tju  >  0  (1 .  3.  6)  represents  one  solution  of  (1.  3.  5)  with  specific 
constants;  and  when  uj  <  0  (1.  3.  7)  also  needs  specific  constants.  If  uu  (t) 
changes  sign  in  the  interval  of  interest,  then  the  requirement  that  (1. 3.  6)  and 
(1.3.  7)  must  represent  the  same  solution  in  the  entire  interval  thus  correlates 
the  forms  (1.  3.6  )  and  (1.  3.  7)  and  the  correlation  is  determined  by  the 
association  of  the  respective  constants.  As  uj(t)  changes  sign  it  is  seen  that  both 
the  forms  (1. 3.  6)  and  (1.  3.  7)  break  down  for  two  reasons: 

(i)  both  become  infinite  when  u>  (t)  vanishes 

(ii)  the  equation  for  which  (1.  3.  6)  or  (1.3.7)  are  exact  solutions  has 
singularitie s  at  the  zeros  of  a  (t)  and  the  functions  (1 .  3 .  6)  and  (1.3.7) 
are  multivalued  in  general  in  the  vicinity  of  these  singularities, 

i.e.  transition  points.  These  points  are  also  called  "turning  points." 
The  representation  of  single -valued  functions  by  multiple -valued  functions  can 
be  expected  to  be  valid  only  in  a  restricted  region. 

In  fact,  as  Langer  says  (Ref.  ZZ),  merely  because  the  pair  of  solutions 
y (t)  ~  y+  (t;  y ,  f- ) 


y(t)  ~  y  (t;o ,  0  ) 


(1.3.  8) 


valid  respectively  on  either  side  of  a  transition  point  exist,  it  is  a  non  sequitur 
that  the  r,  h.  members  of  (1.3,8)  represent  one  and  the  same  solution  of  (1.3.5), 
The  contrary  is  the  case.  For  every  specified  y  ,  6  there  correspond  specific 
a ,  /3  .  In  order  to  deduce  one  form  of  asymptotic  representation  from  the  other, 
one  depends  on  the  so-called  "Connection  Formulae."  These  can  b?  derived  in 
two  ways.  One  is  by  representing  the  solution  near  the  transition  point  by  the 
Airy  function  and  connecting  this  to  the  asymptotic  solutions  on  either  side. 

The  other  is  by  a  study  of  the  Stokes  phenomenon  and  thus  connecting  the 
asymptotic  solutions  on  either  side  of  a  transition  point.  This  will  be  discussed 
in  a  later  section. 


1.4  Objectives  of  the  Investigation 

From  what  has  been  said  hitherto,  it  is  evident  that  for  the  general 
l.d.e.  of  order  greater  than  two,  the  best  general  result  that  can  be 
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obtained  is  to  get  a  good  approximation  with  a  knowledge  of  the  errors 
committed  in  the  use  of  such  an  approximation.  Unless  the  coefficients  have 
certain  special  forms,  it  is,  in  general,  impossible  to  solve  the  equations 
exactly  in  terms  of  elementary  functions  and  operations.  Once  this  is 
realized,  the  aim  is  to  get  approximations  and  error  estimates.  Mathematicians, 
however,  have  been  for  the  most  part  interested  in  areas  which  afford 
general  conclusions  regarding  the  mathematical  properties.  There  is  an 
extensive  mathematical  literature  on  the  many  aspects  of  linear  differential 
equations  and  one  may  refer  to  the  works  of  Hartman  (Ref.  23),  Feschenko 
et  al  (Ref.  24)  etc  .  which  contain  extensive  bibliographies.-  The  engineer 
and  the  physicist,  on  the  other  hand,  have  been  interested  in  approximations, 
insofar  as  they  describe  the  physical  system  adequately,  in  order  to  glean 
some  quantitative  insight  about  the  system. 

The  LG  approximation  seems  to  fill  this  gap  satisfactorily  for  a 
number  of  applications.  Physicists  have  used  the  method  to  great  advantage 
particularly  in  the  fields  of  quantum  mechanics  and  radio  wave  propagation 
in  the  ionosphere  (Ref.  25).  However,  the  control  systems  engineers  have 
generally  stayed  clear  of  this  rather  powerful  method,  except  for  some 
researchers  such  as  Pipes  (Ref.  26),  who  applied  it  to  analyze  time-varying 
networks.  More  recently  Curtiss  (Ref  27)  has  applied  these  ideas  to  the 
analysis  of  VTOL  transition  dynamics  where  he  has  developed  a  modified 
root -locus  method  to  determine  the  "unsteady"  roots,  as  deviations  from 
the  "quasi-steady"  or  the  variable  "charac. eristic"  roots  of  the  system. 

Using  this  technique  one  is  able  to  draw  sketches  fairly  quickly  at  a 
number  of  points  and  obtain  information  about  the  instantaneous  "damping" 
and  "frequency"  of  the  modes  of  motion.  These  applications  have  broken 
the  ice  in  regard  to  engineering  dynamics  analysis  of  variable  systems  and 
pointed  the  way  to  a  more  complete  treatment  of  the  problem.  However, 
with  reference  to  the  VTOL  transition  dynamics,  it  is  desirable  to  have  a 
uniformly  valid  approximation  throughout  the  transition  from  hover  to 
forward  flight.  This  naturally  leads  to  the  study  of  transition  points  or 
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turning  points.  Error  bounds,  of  course,  make  the  theory  more  complete. 

In  contrast  to  the  LG  approximation  for  the  l.d.e.  of  second  order,  for 
which  there  is  a  substantial  body  of  work  covering  these  areas,  no  such 
complete  theory  exists  for  the  third  or  higher  order  l.d.e.  Since  dynamic 
systems  of  the  vehicular  type  are  generally  of  higher  order  than  the  second, 
it  is  felt  that  the  present  work  might  fill  in  part  this  need. 

In  this  thesis  a  scheme  for  approximations  for  l.d.e.  with  variable 
coefficients  is  developed.  Explicit  formulae  for  the  approximate  solutions 
are  derived  and  this  is  done  by  appropriate  extension  of  the  independent 
variable,  employing  multiple  time  scales  and  proper  "clock'1  functions(which  are 
complex  and  nonlinear!  The  "frequency"  and  "amplitude"  variation  of  the 
solution  are  extracted  separately  and  are  then  combined  to  form  a 
composite  solution.  The  advantage  is  that  one  is  able  to  retain,  to  some 
extent,  the  familiar  ideas  of  stationary  linear  systems  analysis. 

Further,  absolute  error  bounds  are  derived.  It  is  clear  that  these 
are  more  useful  than  the  usual  ©symbols  of  asymptotic  analysis,  which  are 
necessarily  somewhat  vague.  The  question  of  transition  points  is  then 
examined  and  a  technique  is  proposed  to  circumvent  the  accompanying 
difficulties. 

1.5  Arrangement  of  the  Dissertation 

The  results  presented  in  the  dissertation  are  presented  in  the  following 
manner.  Chapter  II  presents  the  theory  of  the  method  of  extension 
and  multiple  time  scales  which  will  form  the  basis  of  the  results  obtained  in 
the  dissertation.  The  method  is  applied  to  simple  examples,  and 
asymptotological  principles  are  presented. 

Chapter  III  contains  the  principal  ideas  of  the  thesis.  Here  the  explicit 
formulae  for  the  approximate  solutions  of  l.d.e.  with  variable  coefficients 
are  derived  using  the  method  ofextension.  It  is  seen  that  for  the  second  order 
l.d.e.  the  Liouville -Green  {or  WKBJ)  solution  is  one  of  the  approximations 
derived.  The  formulae  for  third  and  fourth  order  l.d.e.  are  derived  and  then 
the  theory  is  generalized  to  the  n**1  order  equation. 

1  1 


In  Chapter  IV  absolute  bounds  on  the  errors  of  the  approximations  are 
derived.  For  the  second  order  equation. the  error  bounds  for  the  WKBJ 
approximation  are  derived  in  a  new  and  direct  way  different  from  that 
of  Olver  (Ref .  28)  and  asymptoticity  of  the  approximations  is  demonstrated. 

Chapter  V  contains  the  application  of  the  above  approximation  to 
some  analytical  examples  and  the  analysis  of  the  dynamics  of  VTOL 
aircraft  during  transition  from  hover  to  forward  flight.  First  the  hover 
or  two  degree -of -freedom  case  is  studied,  and  then  the  full  three 
degree -of-freedom  system  it  studied. 

A  brief  sketch  of  the  failure  of  the  LG  approximation  near  transition 
points  is  discussed  in  Chapter  VI.  An  outline  of  some  ©pen  problems  is 
presented,  together  with  a  method  of  shifting  the  transition  point  out  of 
the  physical  domain  of  interest. 

In  this  work  the  word  "canonical"  is  used  to  denote  an  equation  of 

t  Vi  tli 

the  n  order  in  which  the  term  containing  the  (n-1)  derivative  does  not 

appear. 
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CHAPTER  II 


ASYMPTOTIC  REPRESENTATION  AND  UNIFORM  VALIDITY 

2.  1  Nonuniformity  in  Approximations 

In  this  chapter,  the  failure  of  the  conventional  perturbation  approach 
in  certain  regions  of  interest  and  some  well-known  methods  of  dealing  with 
it  are  presented  briefly.  The  purpose  is  to  provide  a  basis  for  the  present 
work  within  the  general  framework  of  the  theory  of  such  approximations. 
Formal  proofs  are  not  presented. 

We  shall  begin  by  illustrating  what  we  mean  by  a  uniformly  valid 
approximation.  Given  a  function  f(t)  of  quite  an  arbitrary  shape  (Fig.  1), 
fQ(t)  is  said  to  be  a  uniformly  valid  approximation  of  f(t)  to  order  c 
(where  e  is  a  "small"  parameter,  i.  e.,  e  <<  1 )  if  and  only  if  for  all  t: 

f  =  f  +  O(e);  (or  f=f  +  o(l)  )  (2.  1.  1) 

o  o 

That  is,  the  error  between  the  function  and  its  approximation  is  uniformly 
small  within  the  domain  of  interest.  A  further  discussion  is  given  in 
Appendix  V. 

Contributions  to  the  theory  of  uniformization  of  asymptotic  expansions 
have  come  from  many  sources  and  it  is  difficult  to  do  justice  to  all  of  them. 
The  work  of  some  authors,  however,  has  been  highlighted  for  purposes  of 
orientation.  In  problems  exhibiting  the  presence  of  a  small  parameter  e, 
approximations  based  on  a  direct  perturbation  expansion  in  powers  of  c 
were  first  introduced  by  Poincare  (Ref.  2)  in  his  researches  on  celestial 
mechanics.  Often  such  a  scheme  leads  to  a  serious  misrepresentation 
of  the  true  function  and  this  phenomenon  is  called  nonuniformity  in  the 
expansion.  For  example,  direct  expansion  about  hover,  of  the  solution  to 
the  transition  equations  of  motion  of  a  VTOL  aircraft  fails  to  yield  the 
correct  long  time  behavior.  Also  the  expansion  of  the  solution  of  the 
Liouville  equation  of  statistical  mechanics,  in  powers  of  the  strength  of  the 
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two -body  interaction,  breaks  down  for  times  of  the  order  of  the  relaxation 
time  to  equilibrium;  and  so  fails  to  give  the  crucial  information  on  how  a 
gas  approaches  equilibrium  (Ref.  29). 

The  precise  nature  of  the  nonuniformity  enables  one  to  classify  in  the 
fashion  of  Sandri  (Ref.  30)  as  follows. 

( 1  )Singular  type 

Nonuniformity  occurs  for  finite  values  of  the  independent  variable. 

(2)  Matching  type 

Nonuniformity  is  manifest  in  that  it  is  not  possible  to  satisfy  the 
initial  or  boundary  conditions.  This  is  usually  because  the  inherent 
simplification  of  perturbation  theory  results  in  lowering  the  order  of  the 
original  equation. 

(3)  Secular  type 

Nonuniformity  occurs  for  large  values  of  the  independent  variable. 
The  classification  is  only  pragmatic  and  it  must  be  noted  that 
sometimes  one  type  of  problem  can  be  transformed  into  the  other. 

The  manner  in  which  the  nonuniformities  arise  is  illustrated  as 
follows  with  simple  examples  without  plunging  into  lengthy  calculations. 

1.  Singular  Perturbation 

Since  only  linear  systems  are  of  concern  here,  the  linear  analog  of 
Lighthill's  well-known  example  (Ref.  31,  32)  (Linear  Lighthill  Model) 
suffices  to  illustrate  the  essential  features  of  the  phenomenon. 

Consider  the  equation 


(t  +  c)  —  +  f  =  0  .  (2.  1.2) 

with  the  condition-  f(l)  -  1  (2.1.3) 

Direct  perturbation  theory  yields: 

f  =  f  +  efx  +  .  .  . 
o 
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(2.1.4) 
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Perturbation  theory  yields 


af  '  +  bf  =0 
o  o 

f  "  +  af,'  +  bf  =  0 
o  1  1 


having  the  solutions 

f  =  k  exp  (  -  { b /  a.)t  )  ,  (a) 

o  o 

(2.  1.10) 

f:  =  kt  exp  (  -(  b  /  a  )  t  )  -  k^t  exp  (  -(b/a  )t  )  (b) 

where  k^,  are  arbitrary  constants.  Clearly  conditions  (2.1.V) 
cannot  be  met.  Furthermore,  the  "correction”  r  eventually  becomes 
larger  than  the  lowest  order  term  f  and  therefore  the  expansion  is  not 
uniform  for  large  t. 

Equation  (2.  1.8)  can,  however,  be  solved  exactly  as 
m  t 

f(t)  =  c  e  °  +  c,  e™1  1  (2.  1.  11) 

o 

where  m^,  mx  are  the  roots  of: 

a  a  b 

irr  +  —  m  +  —  =  0 
e  e 

a  i  ,  ,  a  ,a  4b  i 

m  =  -—+$((-)* - ) 

o  e  e 

m,  =  -  2-  -  i  (  (  i  r  -*±  >* 

e  c  c 

c  ,  c:  can  be  chosen  suitably;  the  solution  is  depicted  in  Fig.  3. 
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3.  Secular  type 


Given  the  equation 


df 

dt 


+  r  f  -  0 


(2.1.12) 


with  f(0)  =  i,  a  direct  perturbation  expansion  yields: 


(2.  1.  13) 


The  approximation  f  =  +  r  fails  for  t  — 


the  exact  solution 


f  =  exp  (  -e  t) 


reveals  the  slow  decay  (Fig.  4). 

Some  techniques  have  been  developed  in  order  to  render  the  approximate 
solutions  uniformly  valid.  These  methods  of  uniformization  can  be  broadly 
classified  as  follows. 

(1)  The  Poincare-Lighthill-Kuo  or  PLK  Method 

(2)  Method  of  Matched  Asymptotic  Expansions  (Inner  and  Outer  Expansions) 

(3)  Method  of  Extension  and  Multiple  Time  Scales 

The  PLK  Method  (Ref.  32)  is  typically  applied  to  singular  perturbation 
problems.  The  method  consists  of  suitably  stretching  the  independent 
variable  and  moving  the  singularity  out  of  the  physical  domain  of  interest. 

This  is  done  in  equation  (2.  1.2),  for  example,  by  expanding  the  independent 
variable  t  also  in  a  series: 

t  =  s  +  c  ta  (s)  +  e  ata  (s)  +  ...  (2.1.2  a) 

The  functions  t:  (s)  ,  tj  (s)  ,  etc.  are  to  be  chosen  so  as  to  eliminate  the 
nonuniform  terms.  This  technique  yields  the  exact  solution  for  the  above 
example.  It  is  worth  noting  that  even  though  in  principle  the  PLK  method 
and  the  time  scales  approach  are  similar,  they  in  fact  differ  considerably 
in  the  mechanical  details  of  the  analysis.  In  Appendix  IV  we  emphasize 


I 


t 


this  point  with  a  simple  example.  One  can,  however,  obtain  connection  between 
the  stretching  and  clock  functions,  as  was  done,  for  example,  by  Sandri  (Ref.  3  0). 
Conceptually,  while  the  multiple  time  scales  approach  obtains  the  description 
of  a  phenomenon  by  following  the  gross  features  on  one  scale  and  fine  ones  on 
another,  the  main  purpose  of  the  PL.K  method  is  to  move  the  singularities 
outside  the  domain  of  interest. 

The  method  of  "inner"  and  "outer"  expansions  consists  in  developing 
separate  expansions  (Ref.  33),  each  valid  within  a  region,  and  matching 
them  at  the  boundaries  of  these  regions.  This  method  has  been  typically 
applied  to  problems  of  the  boundary  layer  type.  The  boundary  layer 
approach  and  that  of  multiple  scales  have  as  a  common  feature  the  existence 
of  separate  scales,  on  each  of  which  the  unknown  function  exhibits  different 
behaviors.  However,  in  the  former,  the  "inner"  and  "outer"  solutions 
have  to  be  matched  at  the  common  boundary  or  in  a  region  of  overlap. 

This  usually  calls  for  criteria  based  on  intuition  and  seems  to  involve  a 
certain  amount  of  art  in  the  process.  The  method  of  multiple  scales,  on  the 
other  hand,  while  it  recovers  the  different  behaviors,  does  not  involve  any 
matching,  but  consists  in  an  extension  of  the  independent  variable  into  a 
space  of  more  than  one  dimension.  Further,  it  is  a  formal  method  and 
may  lead  to  a  systematic  method  of  studying  problems  of  the  boundary  layer 
type  also. 

In  regard  to  (2. 1.8)  it  can  be  shown  that  the  fast  variation  of  the 
solution  for  small  values  of  t  can  be  obtained  in  lowest  order,  by  a 
suitable  extension  of  the  variables,  given  by: 

t  f  T_!  ,  T  q  ,  ...  Tn}  (2.1.14) 

f  =  f  +  efi  +  . . . 
o  1 

T .  (t)  are  defined  by: 

t  n 

T  i  -  *  T  ~  t  ,  T  i  =  €  t  ,  . . .  T  ~  G  t 

*  €  o  *  n 
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The  lowest  order  equation  describes  the  "inner"  solution.  To  obtain  the 
behavior  for  values  of  t  of  order  un:-<y  or  greater,  we  do  not  have  to 
consider  an  additional,  "outer",  expansion  of  the  original  equation. 

The  term  of  the  next  order  in  the  expansion  already  made  obtains  the  "outer" 
solution.  The  time  derivatives  are  now  extended  as: 


d 

dt 


1  a  .  a  .  * 

_  - -  +  - -  +  £  - -  +  .  .  . 

e  a  T.a  a  tq  a  r  x 


(2. 1. 15) 


d2 

d? 


+ 


af _ 

ST-i  ar0 


+  2 


af _ 

3t4  ati 


)  +  . . . 


In  lowest  order  we  have 


9TUa  a  rml 


0 


(2.1.16) 


with  the  solution: 


f  (T  -i»  T  »  •  •)  =  -  -CL  (  T  » 
A  o  a  o 


T1( .. )  exp  (-aT.j  )  + 


MV 


T,,..)  (2.1.17) 


The  equation  in  the  next  order  is: 


a3  f 


at,  aT 


+  a  H  +  M  =  0 


(2.1.  18) 


Substituting  from  (2.  1.  17)  above: 


•  •  • 

2  k  exp  (  -ar  ,  )  -  k  exp  (-aT  ,)  +  ak. 


+  b  [ 


a 


exp(  -ar^)  +  1^  ] 


0 
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Regrouping  terms: 


(k  -  —  k  )  exp(  -  at,  )  +  (ak.  +  bk.  )  =  0  (2.1.  19) 

o  a  o 

We  now  equate  the  terms  in  each  parenthesis  to  zero,  which  gives: 

kQ  =  (Tx  ,  .  .  .  )  exp(  +  -  T  )  (a) 

(2.  1.20) 

b 

ki  =  La  (Tj  ,  .  .  .  )  exp(  -  -  T  )  (b) 

Upon  restricting  along  Tml  =  —  ,  r  q=  t  we  have: 

f  =  c;  exp(  -  ~  )  exp(  +  ^  t)  +  ca  exp(  -  ^  t)  (2.  1.21) 

This  process  can  go  on  to  obtain  higher  approximations.  It  can  be  seen  that 
(2.1.21)  describes  the  correct  behavior  of  the  exact  solution  (2.  1 .  11)  to 
leading  order. 

For  purposes  of  the  present  study,  the  method  of  extension  and  multiple  scales 
is  more  pertinent;  the  primary  interest  here  will  be  as  it  applies  to  multiple  time 
scales,  though  in  its  general  concept  (Ref.  29)  the  method  includes  the  other  schemes 
also.  A  further  discussion  is  given  in  the  next  section.  An  example  of  singular 
perturbation  will,  however,  be  considered  and  it  will  be  solved  by  the  time  scales 
treatment. 

The  main  aim  is  to  show  that  the  failure  of  the  direct  perturbation  expansion 
has  as  its  raison  d'etre  an  inappropriate  scale  on  which  the  unknown  function  is 
observed.  The  natural  scales  or  "clocks"  which  afford  a  uniform  description  of  the 
phenomenon  are  determined  by  knowing  the  precise  nature  of  the  breakdown  of 
the  direct  expansion. 

2.2  The  Concept  of  Extension 

The  method  of  extension  was  recently  introduced  as  a  mathematical  technique 
designed  to  exploit  as  much  as  possible  the  presence  of  a  small  parameter  if  one  is 
available  in  a  problem.  The  aim  is  to  render  approximations  of  the  perturbation 
type  uniformly  valid.  The  origin  of  the  concept  can  be  traced  to  the  works  of  Bogoliubov, 
Krylov,  and  Mitropolsky,  who  allowed  all  the  constants  of  the  lowest  order  perturbation 
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theory  to  be  slowly  varying  functions.  The  original  ideas  of  Bogoliubov  and 
Krylov  (Ref.  34)  were  extended  and  modified  in  the  recent  work  (translated  into 
English)  of  Bogoliubov  and  Mitropoldty (R ef.  35),  which  provides  a  broad 
theoretical  framework  for  the  method  of  averaging.  The  method  of  time  scales 
in  its  early  form  was  applied  to  certain  nonlinear  differential  equations  by 
Cole  and  Kevorkian  (Ref.  36);  Frieman  (Ref.  37)  and  Sandri  (Ref.  29)  used  it  in 
the  theory  of  irreversible  nrocesses.  Sandri  (Ref.  38)  has  also  discussed  a 
general  technique  of  un'(  rmization  of  asymptotic  expansions  and  has  shown 
some  of  the  well-known  unixurmization  procedures  to  emerge  as  special  cases 
of  such  a  general  technique,  by  introducing  a  complete  reparametrization  of 
the  lowest  order  term  of  the  perturbation  expansion.  The  development  of  the 
theory  given  in  Ref.  38  is  rather  abstract  and  relies  on  the  composition  of 
mappings.  Precise  conditions  of  validity  of  any  special  form  of  the  method  are 
not  established  and  it  is  here  that  the  present  work  seeks  to  fill  a  gap. 

The  fundamental  idea  is  to  extend  the  domain  of  the  independent  variable 
using  suitable  "clocks"  determined  by  knowing  the  precise  nature  of  the 
nonuniformities  arising  in  direct  perturbation  theory.  It  should  be  noted  that 
the  variables  in  general  are  not  restricted  to  be  real.  The  "clocks"  are  so 
chosen  that  the  new  terms  that  arise  due  to  extension,  called  "counterterms", 
eliminate  the  nonuniformities  of  direct  perturbation  theory  so  that  in  the 
extended  domain  uniform  approximations  to  the  unknown  function  can  be  obtained. 

The  concept  becomes  more  transparent  by  a  re-examination  of  (2.1.  12). 


The  solution  f  is  represented  as 


f  =  1  -  e  t  + 


e  2  ta 


c3t3 


(2.2.1) 


which  is  a  convergent  series  and  can  be  summed  to  the  exact  solution 

f  =  exp(-et).  In  general  the  perturbation  series  is  not  summable  and  one  has  to 

,  ,  th  .  . 

resort  to  the  k  order  approximation: 

23  ,  k  k 

f  ~  1  -  et  +  -  ...  +  (-1)  (2.2.2) 

Clearly  this  fails  for  t  ~  ,  since  all  terms  will  attain  the  same  order  of 

€ 

magnitude.  The  fact  is  that  in  representing  a  function  by  a  series,  we  want  the 
leading  term  to  give  maximum  information  and  hence  we  look  for  an  asymptotic 
expansion  rather  than  a  convergent  one.  These  have  the  property  that  successive 
terms  decrease  in  magnitude  up  to  a  point,  beyond  which  they  may  start 


growing.  The  series  is  terminated  at  this  point  and  provides  a  useful  means  of 
computing  the  functidn  (Ref.  39,40).  With  this  in  mind  we  see  that  the  first  few 
terms  of  the  series  (2.2.2)  do  not  represent  the  true  solution  adequately  for 
long  times. 

A  clearer  physical  picture  comes  into  view  when  we  consider  the  function 


f(t)  =  exp(  -  - —  ) 
1 * 


(2.2.3) 


from  a  different  standpoint,  t  is  a  fixed  constant  with  dimensions  of  time,  and 
f  represents  a  physically  observable  quantity  such  as  displacement  from  a 
reference  position,  or  temperature  difference  between  two  insulated  bodies. 

An  observer  who  measures  f  and  records  it  using  a  clock  whose  unit  of  time  is 
t^  will  have  to  wait  a  long  time  (the  longer  for  smaller  e  )  before  he  can 
observe  a  perceptible  change  in  f.  Instead,  if  our  observer  were  to  use  the 
slow  variable  Tx  =et,  or  a  "super"  clock  which  measures  t  in  giant  units  of 

f! 

the  phenomenon  is  seen  much  better,  for  then  (2.2.3)  is 

f(t)  =  exp  (-  p1  ) 

* 

which  is  indeed  the  exact  representation  of  f.  Thus  the  method  of  extension 
enables  us  to  perform  readings  on  appropriate  scales  by  employing  a  sufficient 
number  of  independent  observers. 

Consider  a  three-dimensional  space  (Fig.  5)  with  orthogonal  axes 
T  ,  Tj  ,  and  f .  Readings  on  "fast"  and  "slow"  clocks  are  represented 
respectively  by  points  along  Tq  and  T  j  coordinates  and  f_  is  defined  to  be 
the  function 


MT  .  T  x  )  =  e  exp(  -T !  ) 


(2.2.4) 


where  c  is  a  constant.  Graphically,  f(T  ,  T*  )  is  represented  by  a  cylindrical 

surface  in  Fig.  5  which  is  constant  in  Tq  ,  but  decays  exponentially  in  T4  . 

To  relate  f  (  T  ,  T .  )  to  f(t)  ,  let  T=  j-  , 

-  o  ** 

From  (2.  2.  3): 

f[  T  ]  =  exp  C  -ct] 


Choosing  c  =  l,  (2.2.4)  gives 
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f  (  T,  e  T)  =  exp(-fT) 


(2.2.5) 


i.e.  ,  f(T,  E  T)  =  £  f  T] 

and  {_  (T  ,  )  is  said  to  be  an  extension  of  f(t). 

Based  on  these  considerations,  Sandri  (Ref.  29)  defines  extension 
formally  in  the  following  manner. 

Definition.  Given  a  function  f(t)v&ere  t  is  in  general  an  n  dimensional 

vector,  and  a  function  fjTj  ,  Ta  ,  ...  of  the  N  independent  variables 

T,  ,  T,  ,  .  .  .  T  (each  of  which  is  an  n  dimensional  vector),  f  is  said  to 
*  N  — 

be  an  extension  of  f  if  and  only  if  there  exists  a  set  of  N  x  n  equations 

T  =  T  (t)  ,  n  =  1,  2,  . . .  N 
n  n 

which  when  inserted  into  f  give: 

t_  (Tj  (t),  Ta  (t),  .  .  .  TN(t)  )  =  f[t] 


The  space  of  N-tuplets  T  -  [  \  ,  Ta  ..  .  }  is  called  the 

extension  of  the  domain  {. T i  and  the  equations  T  =  T  (t)  are  called 

n  rt 

the  "  trajectories"  in  the  extended  domain.  In  dealing  with  differential 
equations,  the  derivatives  and  indeed  thr  'ntire  differential  expression 
itself  can  be  treated  as  a  function  and  can  be  suitably  extended.  Thus, 
given 

T«  — +  *f  =  ° 

one  extension  of  cp  can  be  written  as 


where 


+ 


dT 


Cf 


0 


Tl 


i 


with: 


T  =  t,  Tk  =  et 


In  general  T  =  T  (e^t). 
n  n 

It  is  evident  that  there  are  infinitely  many  extensions  which  correspond 
to  a  given  function.  Two  degrees  of  freedom  are  available:  choice  of  the 
trajectories  and  choice  of  the  extension  itself.  This  freedom  is  utilized 
in  obtaining  an  f  with  simpler  and  smoother  dependence  on  the  parameter 
than  that  offered  by  f,  and  requiring  the  approximate  solutions  to  be 
uniformly  valid  in  the  domain  of  interest.  Figures  6,  7,  and  8  illustrate 
the  concept  of  extension. 

2.3  Application  to  Simple  Examples 

(a)  Equations  with  Constant  Coefficients 

The  theory  discussed  in  the  last  section  will  now  be  applied  to 
simple  examples.  First  l.d.e.  with  constant  coefficients  are  discussed, 
beginning  with  the  first  order  equation;  and  the  method  is  then  shown  to  work 
for  two  special  types  of  equations  with  variable  coefficients.  The  aim  is  to 
extract  the  leading  behavior  of  the  solutions  and  this  is  done  by  an  extension 
of  the  domain  of  the  independent  variable  alone.  Throughout  the  rest  of 
this  work  primes  denote  differentiation  with  respect  to  (w.r.t)  the  independent 
variable. 

(1)  Slow  Exponential  Decay.  Consider  the  first  order  1.  d.  e.  (2.1.12) 
which  is: 

+  cy  =  0  ;  y(0)  =  1  ;  0<G<<1  (2.1.12) 

at 

The  variable  y  and  t  are  extended  as  follows 

y  e=>-  L 

t  <=>  {  Tq  ,  T  i  ,  T  a  ••• 

with  T  =  t,  T  i  =e  t  ,  T  a^T  3  , .  .  =  T  =0.  Then: 
o  n 
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d  ?  *  dTj  *  dT  2  r_  dTn  a  * 

dt  3T  +  3Ti  dt  +  STa  dt  "  '  *  AT  dt  “  aF  +  €  dF, 

o  no 


(2.3.1) 


On  applying  (2.  3.  1)  to  (2.  1.  12)  and  equating  like  powers  of  e.  we  have: 


-y. 

5T 

o 


=  0 


(2.3.2) 


ati 


+  L 


=  o 


(2.3.3) 


From  (2.3.2)  and  (2.3.3) 


X.  (T0>  Tj  )  =  A(Tj  )  =  exp(-Tl  ) 

which  is  the  exact  solution  of  (2.1.  12)  when  is  restricted  along  T  =t; 

T j  =et.  Fig.  9  a  shows  a  schematic  of  the  root  configuration. 

(2)  Second  Order  Equation.  Consider  the  equation  (Fig.  9  b) 

y"+(a  +  e)y’  +  eay  =  0  (2.3.4) 

where  a  is  a  constant  of  order  unity.  Direct  perturbation  theory  results 
in  secular  nonuniformity  as  follows 


etc. ,  giving 


y  =  y0  +  eyi+--- 

y  "  +  ay  1  “  0 
o  o 


yi  "  +  ayj  =  -  (y^  +  ayQ) 


c 

o 

y  - - exp(  -at)  +  Cj 

o  & 

C2 

yi  =  -Cjt - exp(  -at)  +  c3  ,  etc. 

A 
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where  c^,  Cj  ,  c3  ,  c3  are  constants. 

On  the  other  hand,  extension  yields: 


y(t) 


y(T0.  ) 


=>-  f  t  ,  rx  } ;  r  =t.  rx  =  et 


at 


\  +  a  2JL  =  0 


ar0 


(a) 


ilx-  +a  IX  +  ix_  + 


aTo3T^ 


ar,  dr 


ay  =  0 


(b)  (2.3.5) 


+  AX  =  o  (c) 

arx 

Integration  gives: 

y(T  .  Tx  )  =  -A  exp(  -aTo)  +  B(TX  )  (2.3.6) 

Substituting  in  (2.3.5  b): 

A*  exp(  -aTQ)  +  a(B'  +  B)  =  0 

Since  T  and  Tx  are  independent,  A'  =  0  and  B'  +  B=0.  This  gives 

A  =  pure  constant  (2.3.7) 

B  =  Cexp(-r1  ) 

From  (2.  3.  6 )  and  (2.  3.  7)  restricting  V(T  Q>  Tx  )  along  the  trajectories 
T  =t,  r  i  =  c  t,  we  have  the  general  solution  of  equation  (2.3.4)  as 


(2.3.8) 


y(t)  =  cq  exp{-  at)  +  cx  exp(-  et) 

c  ,  Ci  are  constants,  which  is  the  exact  solution, 
o 

Similarly  in  the  second  order  oscillatory  case  of  (i)  low  frequency, 
large  damping;  and  (ii)  high  frequency,  low  damping,  the  choice  of 
simple  time  scales  separates  the  damped  from  the  oscillatory  motions. 

One  can  then  extract  these  separately  and  combine  them  to  give  the  exact 
answer. 

(3)  Third  and  Fourth  Order  Equations  .  Using  the  same  approach  as 
above,  third  and  fourth  order  l.d.e.  with  constant  coefficients  can  be  solved 
exactly  by  a  judicious  choice  of  time  scales.  For  example,  the  following  two 
cases  of  the  third  order  equation  can  be  considered. 

(i)  oscillatory  mode  with  low  frequency  and  low  damping; 
heavily  damped  non-oscillatory  mode  (Fig.  9  c(i)  ) 

(ii)  oscillatory  mode  with  high  frequency  and  large  damping; 
lightly  damped  non-oscillatory  mode  (Fig.  9  c(ii)  ) 

As  an  example  of  (i)  consider  the  l.d.e. 

y"'  +  (a  +  2e)y"  +  2e(a  +  c)y'  +  2a*ay  =  0  (2.  3.  9) 

where  a,  e  are  constants;  a~l ,  0  <  e  «1. 

Direct  perturbation  theory  fails  because  of  the  appearance  of  secular 
terms  as  shown  below. 

y  +  ay  "  =  0 
o  o 

yx  +  ay  =  -  2(  y  "  +  ay  ') 
o  o  7o 

c 

giving  y  (t)  =  -y  exp(  -  at)  +  cx  t  +  ca 

O  a 

yx  +  ayx  =  -  2acx 
Integration  leads  to  secular  terms. 
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The  choice  of  simple  time  scales  as  in  (2.  3.  1)  yields,  in  a 
straightforward  manner 

y=6£lL  exp(-aT  )  +  B(Ti  )  T  +  T( Tj  )  (2.3.10) 

a  o  o 

A'  exp(-aT  )  +  2a(B'  +  B)  =  0 
o 

giving  A  =  pure  constant 

B  =  Dj  exp(  -Tj  ) 

C  =  Da  exp(  (-1  +  i)  Tx  )  +  D3  exp(  (-1  -  i)  Tx  ) 

Upon  restriction  along  T  =t,  Tj  =  et,  we  can  write 

y (t )  =  Cx  exp(-at)  +  Cs  exp(  ( -1+i)  ct)  +  C3exp(  (-1-i)  ct) 

which  is  the  exact  general  solution  of  (2.  3.  9). 

The  procedure  is  the  same  for  fourth  and  higher  order  systems. 

When  the  motion  has  well-separated  modes  a  proper  choice  of  time  scales 
yields  the  correct  answer  in  a  straightforward  way.  A  typical  fourth  order 
example  is  that  of  the  airplane  longitudinal  equations  of  motion  which 
exhibit  two  oscillatory  modes,  one  of  high  frequency,  heavy  damping,  and 
the  other  of  low  frequency,  low  damping,  being  recognized  respectively  as 
the  short  period  and  phugoid  modes. 

Independent  work  in  this  connection  has  been  recently  reported  by 
H.  Ashley  (Ref.  52).  He  considers  the  constant  coefficient  l.d.e. 
describing  the  aircraft  motion  ,  and  obtains  approximate  solutions,  order 
by  order.  He  also  achieves  a  rough  separation  of  the  performance  and 
dynamic  response  problems.  For  both  these  questions,  he  employs  linear 
time  scales,  in  the  fashion  of  Kevorkian  (Ref.  36).  Our  approach  differs 
from  Ashley's  in  that  we  are  able  to  recover  exact  solutions  of  linear 
equations  with  constant  coefficients.  This  is  done  by  choosing  a  proper 
power  of  the  small  parameter  e  as  the  expansion  parameter.  Further  remarks 
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on  the  comparison  of  Ashley's  work  and  ours  are  made  in  Chapter  V, 
Section  5.3. 


The  method  is  illustrated  by  the  following  example  of  the  fourth 
order  system: 


y""  +  2(a+e)y +  (b  +  4a c  +  2^  )y"  +  2e(b  +  2ae)y'  +  2be*y  =  0  (2.  3.  11) 


where  a,b  are  constants  of  order  unity  and  0  <  e  «  1.  Direct  perturbation 
theory  obtains 

y ""  +  2ay  +  by  "  =  0  (2. 3.  12) 

o  o  'o 

yi  ""  +  2ayx  +  byx  "  =  -2(y^'"  +  2ayo"  +  by^)  (2.  3.13) 


etc.  Integrating 


yJO  =  %  exP(  ^jt)  + 
o  mj  mg 


exp(mgt)  =  Cgt  +  C3 


(2.3.14) 


where  CqI  Ci  ,  Ca  ,  C3 


are  constants  and  mj  ,  mg 


are  the  roots  of 


ma  +  2  am  +  b  =  0 


(2.3.15) 


Substituting  in  (2.3.13): 

y1  ""  +  2ayx  +  by"  =  constant 

Integration  clearly  leads  to  secular  terms. 

Extencmg  t  as  in  (2.  3.  1)  yields: 


y(VTi)  = 

A(Tj)  emi  T° 

+  B(TX  ) 

em=To  +  C(TX  )  To  + 

D(Tj) 

(2.3.16) 

2 

+  3a  &3y 

,  ba3y  ,  2,  s2y 

+  b 

n  t?  % 

i  J  d  j  _ 

ST  STj 
o  1 

■) 

t-° 

K 

ST  STX  ST  * 

o  o 

+  bs “ 

o 

-  u  \C .  J 

6 


a4y 


+  6a 


.,3 

5_L 


ST  ST 

o 


*7 

1 


+  6 


ii _ 

K 3Ti 


+  8a 


_ 

&To  S  T» 


2o 


(2.3.  18) 


+  2 


5  2'y 

37  s 


2b  — 


3Ti 


4a  ^ 

a  r 


+  2by  =  0 


Substituting  (2.3.  16)  in  (2,3.17)  and  observing  the  linear  independence  of 
the  exponentials,  A  and  B  are  deduced  to  be  pure  constants  and 
C  =  E  exp(-T!  ).  Substituting  this  in  (2.3.  18)  and  simplifying, 


E  s  0  and  f>  =  Oi  ^ 1  +  ps  e0® 


(2.3.20) 


where  nx  ,  rig  are  the  roots  of  n*  +  2n  +  2  =  0.  The  solution  therefore  is: 


y(To'Tl 


,  m,  r 

)  =  Cj  e  1  °  +  c: 


e"”  T°  +  0, 


,'Hr, 


+  c4e 


not', 


(2.3.21) 


where  ml  ,  m3  and  nj  ,  ng  satisfy 


and 


ma  +  2am  +  b  =  0 


(a) 

(2.3.22) 


na  +  2n  +  2  =  0  (b) 

respectively  (Fig.  9  d).  The  restriction t.Tj  =  $t  obtains  the  exact  solution. 

A  similar  approach  can  be  used  for  higher  order  equations  also.  The 
important  point  to  note  is  the  existence  of  separate  time  scales  as  evidenced  by 
the  presence  of  a  small  parameter  e  .  The  precise  power  of  e  that  appears 
in  a  time  scale  can  be  obtained  by  applying  Kruskal's  principle  of  maximal 
balance  discussed  in  the  last  section  of  this  chapter.  Though  l.d.e.  with 
constant  coefficients  are  not  difficult  to  solve,  the  examples  above  were 
presented  mainly  for  purposes  of  preserving  some  order  in  the  development 
of  the  method,  rather  than  for  pedantry.  Furthermore,  in  high  order  systems 
the  extraction  of  the  different  behaviors  individually  is  useful  in  providing  a 
different  point  of  view  and  may  obviate  to  some  extent  the  labor  of  factoring 
high  order  characteristic  polynomials. 
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(b)  Equations  with  Variable  Coefficients 


(1)  Singular  Perturbation.  The  method  will  now  be  applied  to  a 
problem  of  the  singular  perturbation  type.  The  linear  Lighthill  model  is 
reconsidered  in  light  of  the  theory  of  multiple  time  scales.  The  governing 


equation  is: 


(  t+e  )  —  +  f  =  0  ;  f(l)  =  1 


(2.1.2) 


We  have  already  seen  that  direct  perturbation  is  singularly  nonuniform. 
The  independent  variable  is  extended  as 

t  *=>[  r  .f,  )  with  T  =  c  +  t  or  T=1 
o  1  o  o 

r  j  =  e  k(t) 

where  c  is  a  constant  and  k(t)  is  as  yet  an  undetermined  clock  function. 
Case  (i)  c  =  0(1 ) 

The  extended  equations  are  f  to  order  e  : 


(T  -  c )  +  f  =  o 

o  3  T  -o 

o 


(2.3.23) 


<Vc>  y  *  ■ 

O 


1°  +  k  (  T  -c ) 

Tn  °  STX 


(2.3.24) 


Integrating  (2.3.23)  : 


f  =  A(.T x  ) 

— o  T  -  c 

o 


If  fQ  +  c  fx  is  to  be  an  approximation  to  f  uniformly  then  the  ratio  must 
be  uniformly  bounded.  On  integrating  (2,  3. 24)  this  ratio  can  be  written  as : 


■K1 


— )>)dTot  A{T'] 


(2.  3.25) 


The  counterterm  k  must  be  chosen  so  as  to  cancel  the  nonuniformity  arising 
in  direct  perturbation  theory.  Setting  the  integrand  to  zero  is  sufficient  to 


31 


ensure  the  boundedness  of  ~  uniformly  in  T  .  This  proves  to  be  convenient 

f  o  r 

as  it  enables  us  to  determine  A(TX  )  and  k(T  ).  However,  it  must  be  noted 

o 

that  this  is  only  a  particular  choice  and  the  freedom  in  the  choice  of  the  clock 
function  can  be  exploited  in  other  ways  also.  We  require  uniformity  of  the 
extended  function  along  the  trajectories  T^ft)  and  Tx  (e  ,  t)  .  But  one  may 
for  example  demand  uniformity  in  the  T  ,  Tj  plane.  In  this  case  A(Tx  ) 
should  not  vanish  at  T  x  =  0.  The  behavior  of  B(Tl  )  must  be  determined  by 
going  to  the  next  order. 

Further,  we  may  note  that  if  only  the  independent  variable  in  (2.1.2) 
had  been  extended,  the  above  condition  would,  of  necessity,  have  to  be 
satisfied . 

Thus  we  can  write 


A(TX  ) 

<vT? 

o 


kA' 


(2.3.26) 


i. 


e.  : 


Hence: 


X-  (T>)=EFT-Ci»  -  =  constant 

o 

A  -  D  exp(cx  T  l  ) 


_1 _ 

(T  -  c  ) 
o 


(2.3.27) 


After  restricting  Tq  =  t+c  ,  Tx 


c  k(t)  we  can  write: 


f 

— o 


(T 


D 

t 


e 


-c/t 


(2.3.28) 


This  is  an  improvement  on  direct  perturbation  theory  in  that  f^(t) 
is  finite  at  the  origin;  however,  it  is  not  very  useful  as  f^t)  i*  forced  to 
go  to  zero  at  t  =  0. 

We  shall,  therefore,  consider  the  next  case. 

Case  (ii)  c*0(c);  c=«a,  a  =0(1) 
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The  extended  equations  now  are: 


*  f 

T  +  f  =  0 

o  o  T  o 

o 


(2.3.29) 


r  +f,  = 

o  ST  1 


(1  .  *fn  ;  df 

1  -a  — 0  +  r  k  — — o 
ST  o  at, 

o  1 


(2.3.30) 


As  before: 


f  (r.  r,  l. 
o  o  T 


The  uniformity  ratio  is  given  by: 


i  .f 

f 

O  o 


/  A’  L 
(  A"  ^ 


TTT  )  dt  +  ~  (T  j  ) 
T  o  A  1 
o 


(2.3.  31) 


The  second  term  in  the  integrand  gives  rise  to  singular  nonuniformity  in 
straight  perturbation  theory.  In  order  to  eliminate  this  we  may  put  a=l 
and  k=0.  Exact  solution  is  now  obtained  and  is  given  by: 


,  .  .  D 

f  (t)  = -  ;  D  =  constant 

o  t  +  f 


(2.3.32) 


On  the  other  hand,  we  may  equate  the  integrand  in  (2.3,31)  to  zero  and  solve  for  k. 
Then , 

-(-‘---a)TAr> 1  ■  iviT.) 


and 


A1  ,  ,  (1 -a)  ,  . 

T~  Tl  )  =  7TT 3""  (T  )  =  cl  =  constant 
A  *  kT  o  1 
o 


(2.3. 33) 


giving 


A  =  D  exp(cl  Tj  ) 
a  -1 


k  = 


c,  T 

*  n 


After  restriction  along  Tq=  t+c,  Tj  =  k{t )  : 


C  t.  ,  \,  ,  D  ,  ,  e  (a  -1  )  ,  D  ,  ,€(3-1), 

o(To'T‘  >|t  1  777“  >  =  'tTTI1  exp(TT7T 


(2.3.34) 
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On  substituting  (2.3.34)  into  (2.1.2)  the  value  of  a  is  fixed  to  be  equal  to  1 
and  the  exact  solution  (2.3.32)  is  obtained  (Fig.  2). 

It  seems  desirable  to  develop  a  criterion  of  uniformity  in  terms  of 
conditions  on  the  time  scales  that  would  enable  one  to  proceed  systematically. 
One  may  therefore  consider  the  following  criterion: 

Clock  Uniformity  Criterion  (CUC).  The  time  scales  T  (t) 

- L -  o 

and  Tx  (t)  must  be  chosen  such  that 


i  run' 

'  T'  (t) 

o  i 


i.e.  ,  the  slope  of  the  T  a  (T  )  curve  must  be  0(r  )  uniformly  in  t 
when  the  parameter  e  is  used  to  separate  the  time  scales. 

In  the  light  of  this  criterion,  we  may  note  that  the  clock  uniformity 
ratio  (CUR)  fo-  case  (i)  of  this  example  is  given  by 


rj  (t)  _  e 

|r;(t)  -7T? 

and  therefore  the  CUC  cannot  be  satisfied  for  all  t.  However,  for  case  (ii) 
the  CUR  is  given  by 


Tj  (t)  J  _  c  (a-)  ) 

T '  (t)  I  “  (t  +  r  a)3 

o  | 

If  we  now  demand  that  CUC  be  satisfied  for  all  t  we  must  have  a  =  1. 
This  leads  to  the  exact  solution  (2.  3.32)  without  having  to  substitute 
(2.  3.  34)  back  into  (2.1.2). 

It  has  thus  been  demonstrated  that  a  singular  perturbation  problem 
can  sometimes  be  solved  by  a  proper  choice  of  time  scales. 
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(2)  Simple  Dynamic  Model.  The  following  example  is  given  here 
because  it  is  the  forerunner  of  a  phenomenon  of  deep  significance  in  the 
approximation  of  solutions  of  higher  order  equations.  For  instance,  in 
second  order  l.d.e.  the  change  of  sign  of  the  coefficient  is  associated  with 
the  failure  of  the  approximation  and  the  related  analysis  of  the  Stokes 
phenomenon.  The  present  example  exhibits  the  breakdown  of  the  "frozen" 
and  the  perturbation  approximations.  The  example  may  be  considered 
as  the  simplest  model  of  a  flight  vehicle  of  the  VTOL  type,  being  characterized 
by  initial  apparent  instability;  the  system  is  stable,  however,  for  long 
times.  Again  the  equation  is  solved  by  the  method  of  extension. 

Consider  the  equation: 


df 

dt 


/  1  ~  C  t 

'  1  +  e  t 


)  f  = 


0 


(2.3.35) 


f(0)  r  1;  0  <  e  <<  1 

Fig.  10  illustrates  the  variation  of  the  characteristic  root  and  the  solution. 

The  simple  "frozen"  approximation  is  a  growing  exponential  and  does  not 
match  the  true  solution  anywhere  except  near  t  =  0  and  gives  incorrect 
stability  information.  Another  approximation,  which  is  a  slightly  more 
relined  scheme  of  "freezing"  the  system,is  to  treat  the  coefficient 
essentially  as  a  constant  is  far  as  the  solution  is  concerr  d,  but  to  vary 
on  a  slower  ct  time  scale,  and  can  be  viewed  as  a  simple  application  of 
the  time  scales  method.  The  approximation 

r (t )  =  exp  {  (  ~ ~  )  t  } 

I  +r  t 

thus  obtained  gives  the  correct  initial  behavior  and  stability  information,  but  is 
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quite  wrong  in  representing  the  true  solution  in  other  respects.  For  example, 

the  correct  asymptotic  behavior  for  large  t  is  not  described;  besides,  the 

maximum  value  which  occurs  at  t  as  0.4/  e  for  the  approximation  F  (t)  ,  is 

given  as  exp(0.  17/  e  )  whereas  the  true  maximum  is  exp(0.4/  c  )  and  occurs 

at  t  =  l/c  ;  further,  f  IT  as  exp(  0.2 /e)  .  Direct  perturbation 

itiex  max 

expansion  on  the  other  hand  is  secular  and  yields 

f(t)  =  et  {  1  -  et  +  . . .  ) 

1 

and  fails  for  t  ~  —  . 

f 

We  shall  now  see  that  a  proper  choice  of  time  scales  results  in  a 
uniformly  valid  solution. 

The  variables  are  extended  as  follows: 

*  =*-  (  VT> 1 

with  T  =1;  T  =  t  +  constant 

o  o 

Tx  =“  k(t )  (2.3.36) 

€ 

Now: 

*  U  *- 

dt  3T  e  ST, 

o  * 

Equation  (2. 3.  35)  can  be  written  as  : 

(-  +t)^f  -  (  -  -t)  f  =  0  (2.3.37) 

g  dt  c 

This  suggests  that  the  constant  in  T  in  (2.  3.36)  is  0(— ).  Therefore, 
c 

let  T  =  t  +  —  ;  c  =  0(1).  The  extended  equations  are: 


The  choice  of  c 
equation  (b) 


d-c)  k|=  =0  (a) 

d  T  i 

(1-c)  -  U+c)f  +r  k  ^  =0  (b)  (2.3.38) 

STq  -  o  STi 

T  (  +  f  )  =  °  (c) 

o  3  T  — 
o 

from  (a)  and  f  =  A (T,  )  exp(-T  )  from  (c)  yields  in 

o 


2 

r  £ 


(r  )  =  c 
o 


l 


=  constant 


whence : 


A  =  D  exp(c,  Tt  );  k(T  )  =  —  /n  T 
4  O  1  C|  o 


Restricting  f  (T  ,  T,  )  along  T  =  t  +  —  ,  tj  =  —  k  and  f(0)  =1,  we  obtain 

—  o  *  o  e  e 


2  ..  ,  -t  ...  2/c 


f(t)  -  exp{  -t  +  —  {,n  (1+Gt)]=  e  (1+et) 


which  he  exact  solution.  The  asymptotic  behavior  of  the  function  can  be 
wr  on  ..  ^ : 


tt 


f(t) 


f(t) 

t“*OD 


1 


f  2  ,  1  1  ,  1  si 

expt  T  (  ~  '  27*1*  +57rP 


f(t) 

tj  t-0 


exp(  t  -  c  t“  +  —  c  t 


2  _3.3  e3t* 


+  .  .  .  ) 


It  has  thus  been  demonstrated  that  in  dealing  with  equations  having  variable 
coefficients,  the  generality  of  a  nonlinear  clock  function  is  mandatory.  The 
clock  itself  can  be  a  highly  nonlinear  function  even  in  simple  problems. 

2.4  Asymptotology 

This  chapter  concludes  with  a  brief  look  at  one  aspect  of  asymptotic 
analysis  which  has  hitherto  been  known  as  an  art,  at  best  as  a  quasi -science. 
Most  people  who  have  worked  with  asymptotic  phenomena  have  acquired 
implicit  knowledge  useful  with  different  problems  but  not  general  enough  to  be 
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explicitly  formulated. 


In  a  highly  instructive  lecture  M.  Kruskal  enunciated  seven  principles 
governing  the  philosophy  of  approach  in  asymptotic  analysis  (Ref.  41).  This 
section  (beginning  with  the  title  which  was  first  used  by  him)  is  a  brief  review 
of  these  principles  motivated  by  their  usefulness  in  later  sections. 

Kruskal  defines  asymptotology  as  the  art  of  dealing  with  applied 
mathematical  systems  in  limiting  cases;  alternatively  it  is  the  art  of 
describing  the  behavior  of  a  specified  solution  (or  family  of  solutions)  of  a 
system  in  a  limiting  case.  The  principles  are  enumerated  below;  however, 
the  one  most  important  for  the  present  work  is  the  sixth--the  principle  of 
maximal  balance  (  or  minimal  simplification  ). 

1.  Principle  of  Simplification. 

Asymptotological  analysis  tends  to  simplify  the  system  considered, 
thus  facilitating  the  generation  of  approximate  solutions.  Simplification 
occurs  for  example  in  perturbation  theory;  another  way  this  can  occur  is  in 
the  separation  of  autonomous  subsystems.  The  system  f(x,y)  =  0;  g(x)  =  0 
has  the  autonomous  subsystem  g(x)  =  0. 

2.  Principle  of  Recursion. 

The  dominant  terms  only  are  retained  and  solved  for  and  the  other 
terms  are  treated  as  known.  Iteration  enables  one  to  obtain  an  asymptotic 
representation  of  the  unknown  function  irrespective  of  the  forms  of  the  terms 
appearing.  This  principle  is  also  useful  in  deducing  general  properties  through 
mathematical  induction. 

3.  Principle  of  Interpretation. 

This  advises  us  to  suitably  formulate  the  problem  so  that  the  limiting 
case  is  meaningful.  Overdeterminism  as  occurring  in  matching  problems 
of  the  boundary  layer  type,  results  in  simplification  at  the  cost  of  losing 
importart  information. 

4.  Principle  of  Wild  Behavior. 

This  states  that  apparent  overdeterminism  occurs  if  in  the  limit  (at  least) 


38 


some  solutions  have  peculiar  (e.g.  singular)  behavior.  For  instance,  if 

f(t)  =  exp  (  -  —  )  ,  f  -*-<*>  when  t  vanishes  faster  than  c.  More  general 
e 

forms  must  be  used  in  the  asymptotic  representation  (e.g.  inverse  powers 
or  logarithmic  terms).  Underdeterminism  on  the  other  hand  results  in 
nonuniqueness  of  solutions. 

5.  Principle  of  Annihilation. 

An  annihilator  is  an  operator  which  results  in  zero  when  applied  to 
a  mathematical  entity.  It  is  used  to  eliminate  such  terms  as  may  lead  to  under  - 
determinism  in  the  limit. 

6.  Principle  of  Maximal  Balance. 

This  dictates  the  choice  of  the  terms  to  be  neglected  (leading  to 
simplification)  from  among  the  competing  terms  when  a  comparison  of  the 
relative  asymptotic  magnitude!  is  made.  This  is  based  on  sound  sense  because 
neglecting  the  minimum  number  of  terms  retains  maximum  information. 

When  there  is  more  than  one  maximal  set  of  terms,  each  set  describes  one 
asymptotic  behavior.  Simply  stated,  the  principle  requires  that  the  ordering 
be  so  chosen  that  the  maximum  number  of  terms  is  retained.  For  example, 
in  the  asymptotic  analysis  of  the  roots  of  the  cubic  equation 


3rax3  +  x3  -  ex  -  4=  0 


(2.4.  1) 


in  the  limit  e”*0,  one  may  choose  a  general  representation  x  =  ae  V  and 
determine  v«  The  terms  can  be  ordered  as: 


3\<  +2  Z\; 

c  :  c 


v  +1  o 
c  :  c 


(2.4.2) 


V  is  chosen  such  that  the  maximum  number  of  terms  (which  are  dominant)  is 
retained  after  neglecting  the  terms  which  are  small.  Among  the  different 
choices  of  v,  it  is  found  that  the  maximal  ordering  is  given  by  v=0  or  -2. 

A  graphical  technique,  which  can  be  traced  to  Newton  but  used  in  this  context 
by  Kruskal,  can  be  exploited  to  determine  the  value  of  v  for  maximal 
ordering  (Fig.  11a).  Each  term  of  (2.  4.  1 )  is  represented  as  a  point  on  a 
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graph  with  the  power  of  x  as  the  abscissa  and  the  power  of  r  as  the  ordinate. 
The  relation  between  the  exponents  of  x  and  r  is  given  by  a  line  and  for 
small  e  the  values  of  v  for  minimal  simplification  are  given  by  the 
lower  convex  support  lines  of  the  set  of  graphed  points;  i.e.  lines  passing 
through  at  least  two  points  such  that  there  are  no  graphed  points  below  them. 

The  present  example,  shown  in  Fig.  11a,  is  discussed  in  detail  by  Kruskal 
in  Ref.  41.  For  large  e  ,  on  the  other  hand,  the  upper  convex  support 
describes  maximal  balance.  These  ideas  are  further  discussed  in  the  next 
chapter. 

7.  Principle  of  Mathematical  Nonsense 

This  is  the  simple  idea  that  if  during  an  asymptotological  analysis 
an  absurd  conclusion  is  reached,  the  analysis  has  not  been  done  correctly 
or  carried  far  enough. 

During  the  ensuing  chapters,  more  than  one  of  the  above  principles 
will  be  invoked  and  this  is  the  reason  for  their  inclusion  here.  To  be  sure, 
it  is  more  desirable  to  arrive  at  conclusions  in  a  systematic  and  logical 
manner  after  proper  asymptotological  analysis --as  Kruskal  says--"like 
remarkable  coincidences  in  a  well -constructed  mystery  story,  "  (Ref.  41). 

2.5  Summary  of  the  Chapter 

Nonuniformities  of  perturbation  theory  and  some  well-known  uniformizing 
methods  are  discussed.  The  difference  between  the  PLK  method  and  the 
multiple  time  scales  approach  is  emphasized. 

The  method  of  time  scales  is  applied  to  simple  examples.  Constant 
coefficient  l.d.e.  are  solved  exactly  by  this  method,  the  fourth  order  equation 
being  representative  of  aircraft  motion.  A  singular  perturbation  problem 
and  a  simple  dynamic  model  are  also  studied. 

Asymptotological  principles  are  enumerated. 
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CHAPTER  III 


DEVELOPMENT  OF  THE  APPROXIMATION  SCHEME 


This  chapter  seeks  to  obtain  approximations  to  linear  differential 
equations  with  variable  coefficients  in  various  limiting  cases.  The 
coefficients  are  assumed  to  be  slowly  varying  and  the  precise  definition  of 
slowness  is  discussed  in  the  next  chapter.  The  equations  are  parameterized 
by  introducing  an  ►:  and  the  system  is  studied  in  the  limits  of  small  and 
large  e  .  In  each  case  the  choice  of  suitable  "clocks"  results  in  the  extraction 
of  the  leading  behavior  of  the  solutions;  composite  solutions  are  obtained  by 
combining  the  behaviors  on  different  time  scales.  One  of  the  major  tasks 
of  asymptotic  analysis  is  the  determination  of  the  "natural"  variables  in 
which  the  given  problem  can  be  treated  as  a  perturbation  problem.  The  present 
approach  is  intended  to  relax  the  requirements  on  knowledge  given  in  advance 
and  ad  hoc  assumptions  and  to  provide  a  systematic  way  to  deal  with  equations 
as  they  are  given.  Throughout  this  chapter  the  domain  of  the  independent 
variable  alone  is  extended  into  several  dimensions  and  nonlinear  clock  functions 
are  employed.  The  coefficients  of  the  differential  equations  can  in  general 
depend  on  r  and  t;  however,  in  this  analysis  they  are  assumed  to  depend  only 
on  t.  The  discussion  begins  with  the  first  order  equation  and  is  continued  to 
higher  order  equations. 

3.  1  First  Order  Equation 

Consider  the  equation: 


v'  +  c  ^Q(t)y  =  0  ;  y { 0 )  =  1 


(3.1.1) 


Direct  perturbation  expansion  y  =  y^  +  r  y^  +  ...  exhibits  both  secular  ard 
singular  nonuniformity  depending  on  a  (t).  The  nonuniformity  ratio  often 
indicates  the  nature  of  the  breakdown.  It  is  seen  that 


XL  = 


y 


o 


(3.1.2) 
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Clearly,  for  example^  when 

u)  ^  =  t°  ;  n  >  -1  perturbation  theory  is  secular 


exact 


=  t  ;  n  <  -1  "  "  singular 

=  eint;  n  real  "  "  neither  secular  nor  singular 

In  order  to  uniformize  the  perturbation  expansion,  the  following 
extension  of  the  independent  variable  is  made: 


int 

=  e  ;  n  real 


f  T  ,  Tx  } 
o  * 


ro=t;  tj  =f.k(t) 


(3.1.3) 


The  coefficient  is  taken  to  vary  on  the  T  ^  scale  (to  pick  up  the  varying 
nature  of  the  coefficient)  .  Thus 

i-^Tr  *  cMTc» 

O 

i*  =0 

?r°  , 

MTo>^  +“o<To)>r  =  0 


giving 


y(ro.  rx  )  =  a(t x  ) 


k  (T  )  A'(TX  )  +»„(T  )  A(Tj  )  =  0 
o  *  o  o  * 


A  * 

—  (T x  )  =  -  -J*  (Tq)  =  s  =  constant 

the  1.  h.  s.  and  r.h.  s.  being  respectively  functions  of  Tx  and  T  only.  Hence: 


y(To»Tx  )  =  A(rx  )  =  c  exP  («Ti  ) 

k(T  )  =  -  —  f <o  (T  )  dt 
O  8  J  O  O  O 


(3.1.4) 


Now  restricting  the  solution  along  the  trajectories  Tq=  t  and  rx  =  rk(t)  the 

exact  solution  y(t)  =  c  exp  (-r  I  ®  (t)  dt)  is  recovered.  Henceforth  s  can  be 

J  o 

set  equal  to  unity  without  loss  of  generality. 

We  may  observe  that  in  this  extension  the  CUR  is  given  as: 


T  i  (t) 

T '  (t) 
o 


eu>(t) 


Clearly  the  CUC  cannot  be  met  as  u>(t)  may  have  singularities  in  the  domain 
of  interest.  This,  however,  does  not  really  matter  as  this  difficulty  can  be 
removed  by  making  TQ(t)  also  a  nonlinear  function.  For  example. 
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if  a) 


— >-  —  and  u>  is  bounded  for  large  t,  we  may  choose  T  (t)  =  t  +  (,n  t 

t^  t  o 


and 


T{ 

C  0) 

T  ' 
o 

1  +  1/t 

0 


The  choice  of  a  nonlinear  T  (t)  thus  enables  us  to  enforce  the  clock 

o 

uniformity  criterion. 

It  is  interesting  to  consider  a  pictorial  representation  of  the  function 
(Fig.  7,8).  For  small  e  .  T  q  represents  the  fast  time  scale  and  Tj  the 
slow  time  scale.  The  extended  function  surface  is  essentially  a  constant 
along  t  ,  but  decays  exponentially  on  the  Tj  axis  and  can  therefore  be 
described  naturally  (and  uniformly)  as  dependent  on  Tx  .  The  natural 
clocks  are  depicted  as  trajectories  along  which  the  extended  function  is 
restricted  and  are  solely  determined  by  the  coefficient.  For  instance, 
when  u)q  is  unity  the  trajectory  is  a  straight  line  through  the  origin  at  an 
angle  whose  tangent  is  c  ;  a  linear  u>Q(t)  leads  to  a  parabola  for  the 
trajectory.  In  general  the  magnitude  of  the  parameter  governs  the  proximity 
of  the  trajectories  to  the  t  or  axis. 


3.2  Second  Order  Equation 


In  the  previous  case  it  was  evident  that  whatever  the  nature  of  the 
coefficient  u>0(f)  and  the  magnitude  of  the  parameter  e.  one  could  always 
determine  a  natural  clock  to  describe  the  function  uniformly.  For  the  second 
and  higher  order  equations  this  is  no  longer  possible,  and  one  has  to  be 
content  with  approximations  valid  in  different  regions  of  the  domain  of 
interest.  First  the  canonical  form  of  the  equation  and  then  the  noncanonical 
form  are  discussed. 


(a)  Canonical  Form 

Consider  the  equation: 


y"  +  eu)(t)y  =  0 


(3.2.1) 
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Again  direct  perturbation  can  be  shown  to  fail  for  various  t  (t),  the  nature 
of  nonuniformity  depending  on  uj(t)  (Table  I).  We  wish  to  make  use  of  the 
simplification  afforded  by  perturbation  theory  and  improve  upon  it.  The 
extension  sought  is: 


Ti  ) 


T  =  t; 

o 


T,  = 


-  m  u 

=  e  k, 


(t) 


(3.2.2) 


The  extended  derivatives  are  given  in  Table  U..  Equation  (3.  ?.  1 )  can  be 
written  as : 


3  T 


4- 


m  / 

e  ( 


a_L_ 

ST, 


+  21^ 


sa  y 


,  2m  •  a 
)  +  C  (kta 


+c«|T„ly  =0 


(3.2.3) 


The  terms  containing  kj  are  the  counterterms  introduced  by  extension; 
we  shall  determine  the  clock  functions  in  various  limiting  cases  in  the 
following  way. 

The  various  terms  can  be  written  as: 

o  m  .  2 m  .  ,  n,  ,  . .  .  A . 

c  (  )  +  e  (  )  +  €  (  )  +  e  (  )  =  0  (3. 2. 4) 


n  is  a  given  constant  and  m  is  to  be  determined. 

The  quantities  not  containing  e  are  implied  t  a  be  of  order  unity.  The 
failure  of  the  approximation  is  indicated  when  the  above  ordering  breaks  down. 

In  order  to  use  the  graphical  technique  the  coefficient  of  m  in  the 
exponent  of  e  in  (3.2.  1)  is  plotted  along  the  abscissa  and  the  constant  term 
in  the  exponent  of  c  along  the  ordinate.  Thus  each  term  in  (3.2.4)  denotes  a 
graphed  point.  In  the  fashion  of  Kruskal,  for  small  c  ,  the  lower  convex 
support  line  of  the  set  of  graphed  points  gives  m  for  maximal  balance. 

On  the  other  hand  the  upper  convex  support  line  yields  maximal  balance  for 
large  e  .  In  the  present  context  we  consider  e  to  be  small  and  thus  m=0 
for  maximal  balance.  But  this  choice  is  rejected  as  it  is  not  useful  and 
corresponds  to  straight  perturbation  theory;  i.  e.  no  extjnsion  has  been  made. 
What  is  desired  is  a  compromise  between  completeness  and  simplicity--  a 
system  as  complete  as  possible  and  still  simple  enough  to  solve.  With  this  in 
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mind  we  pick  m  =  l,  which  corresponds  to  "submaximal  balance,  "  i.e. 
maximal  balance  in  second  rank  of  terms. 

We  may  therefore  venture  to  state  below  the  principle  of  submaximal 
balance.  If  the  maximal  balance  of  terms  results  in  an  equation  that  is  either 
too  difficult  to  solve  or  yields  too  little  information  then  submaximal  balance 
obtains  maximum  information  consistent  with  simplification.  Again  this  is 
mere  prudence;  short  of  solving  the  complete  system  (which  may  be  impossible) 
we  take  the  next  best  course  of  action  and  s^lve  a  system  which  is  not  the 
most  complete  but  one  next  to  it  in  the  order  of  importance.  That  is,  the 
system  chosen  contains  more  information  than  all  the  others  except  the 
maximally  complicated  one,  i.  e.  a  system  maximal  in  second  rank  of  terms. 

It  appears  that  it  is  precisely  the  method  of  extension  which  permits 
uniformly  valid  approximations  which  would  otherwise  not  be  possible.  In  the 
above  example  T  =  t  must  be  a  time  scale  since  u‘  is  a  general  function  of  t. 
But  since  the  parameter  is  present  it  is  likely  that  there  is  another  time 
scale.  Unless  u  (t)  is  a  known  function  of  t  we  cannot,  in  general, 
transform  e  out  of  the  equation.  Maximal  balance  is  therefore  not  possible 
without  extension  and  hence  submaximal  balance  is  resorted  to  after  a 
suitable  extension  of  the  variables. 

In  terms  of  the  graphical  technique,  submaximal  balance  for  small  r 
would  correspond  to  the  lower  convex  support  line  of  the  set  of  graphed  points 
except  for  one  point  which  may  be  beneath  the  support  line  (Fig.  11  b). 

For  large  c  this  corresponds  to  the  upper  convex  support  line  except  for 
one  point  which  may  lie  above  the  support  line. 

With  reference  to  (3.2.  3)  we  see  that: 

c  small  m  =  0  obtains  maximal  balance 

m  =  1  obtains  submaximal  balance 

c  large  m=^f  obtains  maximal  balance 

m=  0  or  1  obtains  submaximal  balance 
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(i)  Short  Time-Analysis.  From  the  above,  for  small  e  the  choice 
m  =  1  leads  to  the  following  set  of  partial  differential  equations  replacing 
equation  (3.2.  1): 


14- 


(a) 


k  +  21^  L_L_  +  u,(T  )y  =  0 

i  ati  ^  aT0STi  ° 


(b)  (3.2.5) 


i^a  f-**  =  o 

^  STx  3 


(c) 


Integrating  (3.2.5  a)  we  have: 


y(To  •  Tl)=A(r1)To  +  B(Tj  )  (3.2.6  a) 

Now  A(T x  )  Tq  and  B(Tx  )  are  linearly  independent  v/.r.t  Tq  and  can  be 

used  to  generate  separately  the  corrections  to  the  lowest  order  result, 

giving  rise  to  two  clock  functions.  Substitution  of  (3.2.6  a)  in  (3.2.5  b)  leads  to: 


k*  i  i  A'T  +  2  itjj  A'  =  -u)  (T  )  AT 
4  *  o  o  o 

The  choice  of  A  =  e*^1  solves  this  equation  as: 


(3.2.6  b) 


(a) 

(3.2.7) 


and  similarly: 

r 

skxa  =  - 

J 

The  parameter  s  can  again  be  set  equal  to  unity  without  loss  of  generality. 


1 


w  (r  )  dT 

o  o 


dr 


=  -r 


L 

•j 


dr 


K 


dr 


(0 
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On  restriction  the  two  independent  approximations  are  obtained  as 


%  sj  Cj  t  exp  f  e(  j-J*  t2  wdt  -  J  tu)dt)  } 
y  2  «  c2  exp  (  e  |  [  a  dt  dt  ) 


(3.2.8) 


In  each  case  the  equation  (3.2.5  c^  which  is  left  unsatisfied,  defines  the  error. 


a  3  3  y 


The  approximation  breaks  down  when  the  neglected  term  3  - la 

j  Ti 

becomes  of  order  —  . 

e 

Higher  order  corrections  are  obtained  by  introducing  more  time  scales. 
These  are  obtained  as  multiplicative  corrections  instead  of  additive  ones  as  in 
direct  perturbation  theory,  and  will  therefore  reflect  improvement  in  some 
types  of  problems.  Thus  if  t  is  extended  as 


{t  «Tj,Tj,...,T  } 

o  n 


with  Tq=  t;  T1=ck1(t);  Ta  =  Ga  kg  (t);  . .  .  ;  Tn=  e  k  (t) 


d  3  I  3  ,  j  !  3  . 

- —  ■=*-  —  +  e  k.  - —  +  c  Kg  —  +  . 

dt  a  T  ^  3TX  dTa 


n  •  a 

.  +  e  k  - — 

n  3T 


da  ^a2 


;  a1 


dta  a 


—TT  +  e  (  k.  —  +  2k,  - -  )  + 

v  o7':  ^  JT  ST, 


e  3  (  ka  -  +  2kg 


i.  aa 


r —  +  2kg  - -  +  k.  a - 

3Ta  3T 3Tj  ^ ara 


+  ...  +  C  (k  +  k  3  | - a  +2 T  k.  k  .  -  )  +  c  *(  )  +  . 

n  ar  n/2  a T  2  i  =  0  i  n-i  3T.  at  .  '  1 

n  i  =  n/2  *  n-i 


The  term  k 


n/2  aT 


=  0  when  n  is  odd. 


The  extended  equations  can  be  solved  order  by  order  as  follows: 
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3  3  v 

r-r-4  -  0;  y(T  ...T  )  =  A(Tx  ,  Ta 
a  Tn  o  n 


TnJ  To  +  B(Tl  ,Ta . Tn) 


C>  fr,  |r?Tl  +  » <1-  0>y  =  o 


Again  using  At  and  B  to  generate  two  independent  approximations: 


kn  r 


3  A 


3  A 


3TX  " 


-  U)T  A(T,  ,  .  .  .  ,  T  ) 
o  n 


Now  choose 

A(T  j  iTj  i  •  •  «  •  T  )  —  A(T  a  » T  3  ,  •  •  <  i  T  )  e 
n  i  n 


(3. 


Therefore: 


kll  Tq  +  2kl  i  =  -  wTo 


Integration  giv^s  the  clock  function  obtained  earlier  (3.2.  7  b). 

The  slower  clock  kg  is  obtained  from  the  equation  of  next  order: 


k31  2X.  +2k, 

1  3ra  a 


T0aT* 


.  i  a 


Substituting  from  (3.2.  10)  and  (3,2.  7  a)  one  obtains 


1 

’  r8 

o 


V  (k^^a  dTo 


(a) 


i.  e. : 


=  r -  f  T  3  {  ki  i  )3  dT  -  f  T ■  (  kj  x  )3  dT 
T  O  o  o  *  o 


o  J 


(b)  (3. 


with: 


A(Tx  ,  •  •  •  T^)  =  Ax  (Ta  >  •  »  •  t  T^  )  e  *  =  Aa  (T  3  »  • « »  »T^  )  e  1  e  (c) 


The  general  result  for  n  >  1  can  be  written  recursively  as  follows. 


(3.2.9) 


2.  10) 


2.  11) 
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(3.2. 12) 


and  is  subject  to  the  conditions  (3.2.  13).  Thus  clock  functions  can  be 
determined  to  define  slower  *ime  scales  and  improve  the  accuracy  of  the 
approximations. 
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The  domain  of  *t  for  which  the  short  time  approximation  is  valid  can 
be  determined  as  follows.  In  obtaining  Equation  (3.2.  8)  the  third  term  of 
Equation  (3.2.3)  was  neglected  in  favor  of  the  second  and  fourth  terms. 
Clearly  failure  of  the  approximation  occurs  when  this  condition  is  violated; 
i.  e.  when: 


The  condition  for  failure  is  given  by 


<k)a  =  — 
e 


substituting  from  (3.2.?)  and  simplifying,  we  obtain: 


/e  (  Jta  (i)  dt)  =  Jta 


D  iff  e  r  entiat  ing 


/etau)=$  ujt3  +  2t  w  ^ 


l.  e. 


1  -3/2.  2  y 

t  in  u)  +  —o)  =  ✓  e 


l.  e. 


d  .  “J-  .  2  y 

~(u)  )  -  —  (  U)  )  =  -  ✓< 


Using  t"a  as  an  integrating  factor,  the  above  condition  becomes: 

1_  i  «“  )  =  -4£- 

dt  't5^  '  t8 


Integrating 
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The  approximation  will  fail  when 


e 


(3.2.14) 


The  criterion  can  be  derived  in  a  different  way  also.  The  approximating 
functions  A(T  i  )T  and  BlTj  )  are  linearly  independent  w.  r.  t.  Tq  . 

Upon  restriction  along  T  q  =  t,  Tj  =  e  k(0*  this  property  may  not  be 
satisfied  throughout  the  domain.  The  approximations  can  therefore  be 
expected  to  fail  in  a  region  where  the  constancy  of  the  Wronakian  (Ref.  59) 
is  destroyed.  From  (3.2.  8)  the  Wronskian  can  be  written  as 

w(%  /}&  )  =[  -1+e  (^-Jt3u)dt  -tjujdt)]  exp(e  (£-  j  t3  u>dt  -  (  twdt  -  JJuodt1'  )) 

i.  e.  to  lowest  order  in  c,  W(yj  /ya  )  i®  a.  constant.  Hence  failure  is 
indicated  when  either  the  exponent  is  of  order  unity  or 

^-Jt3u>dt  -tjuidt  ~ 


i.  e. : 


-2c 


uu  dt  dt 


t 


Differentiating: 

-  2 e t  j u> dt  ~  1  or  -2c^|u)dt  ~  ~ 
Differentiating  again: 


Thus  the  approximation  (3.2.  8)  fails  near  a  value  of  t  for  which  vt8  ~  — 

c 

as  obtained  earlier.  Substituting  this  shows  that  the  exponent  in  the 
exponential  function  of  W(yx  ,  ”ya  )  is  of  order  unity. 
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For  the  Airy  equation,  for  example,  u)=t,  and  the  appr oximations 

,3  .3 

£  t  0  t 

t  exp(  --jy  )  and  exp(  -  — — )  obtained  f* om  (3.2.8)  break  down  according 
to  the  above  criterion,  when  t  ~  e  "l^3  .  Any  attempt  to  improve  upon  the 
approximation  by  going  to  higher  order  in  e  is  foiled  as  the  es  approximation 
also  fails  when  t^e  The  reason  for  this  is  clear.  The  method  then  tells 

us  that  it  is  not  possible  to  effect  any  improvement  by  using  slower  clocks. 

On  considering  the  Airy  equation 

y"  +  e  ty  =  0 

it  is  seen  that  when  t  ~  e”1^  the  parameter  e  is  completely,  removed  from 
the  equation  thus  indicating  a  region  in  which  the  equation  must  be  solved 
exactly. 

However,  for  a  different  us(t)  the  criterion  can  be  utilized  to 

advantage.  For  instance,  if  u)(t)  =  p  ,  the  above  criterion  says  that  the 

approximation  will  fail  when  e~l.  For  small  e  breakdown  is  not  indicated 

and  improvements  can  be  effected  by  going  to  higher  order  terms  in  e  . 

This  is  indeed  the  case,  for  the  exact  solution  for  e  >  j  is  oscillatory 

1  4 

but  it  is  not  so  for  eS 

4 

( ii)  Long  Time  Analysis.  The  interest  is  now  shifted  to  the  long  time 
behavior  of  the  solutions  of  Equation  (3,2.1).  It  is  shown  that  the  LG 
approximation  can  be  derived  easily  by  a  proper  choice  of  time  scales.  The 
reason  for  deriving  this  well-known  approximation  is  not  merely  pedagogical 
but  to  provide  a  systematic  method  for  higher  order  equations  and  to  emphasize 
a  clear  physical  picture  of  the  phenomenon.  In  reference  to  Equation  (3.2.  1) 
the  analysis  is  carried  out  in  the  limit  of  large  e  .  If  tu  (t )  is  an  unbounded 
monotonic  function  for  large  t,  the  correspondance  between  e  and  t  is  clear. 

In  any  case,  the  growth  of  the  magnitude  of  u>(t)  must  be  properly  associated 
with  the  limit  of  large  e  .  It  is  seen  that  in  this  limit,  the  choice  m  =$ 
obtains  maximal  ordering  for  the  terms  in  Equation  (3.2.3). 

The  extension  evolves  as  follows,  wittf  2m  =  1  and  denoting  by  X  , 
and  leads  to  the  following  partial  differential  equations,  which  replace  (3.2. 1). 
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x»  :  k»  l~la  +  a)  (1*o  )y 


;•  2x-  +2k  2 
bfx  a 


(b)  (3.2. 


X°  :  r“pi  =  0  (c) 

3o 

Equation  (a)  above  can  be  treated  as  a  constant  coefficient  one  w.  r.t. 
Tx  and  the  solution  is  given  as 


y (Tq  ,Tx  )  =  a(T  q)  exp(r  x  ) 


(3.2.  16  a) 


whence,  the  clock  k  satisfies  the  equation 


(k)3  +  oj(T  )  =  0 


(3.2.  16  b) 


obtained  by  substituting  (3. 2.  16  a)  in  (3. 2,  15  b).  Substitution  into  (3. 2.15  b) 
yields 

d  .  .  1  d 

a?  It"">  =  T  3F  (t"") 

O  O 


a(T  )  =  u)-1^4  (T  ) 
o  o 

Restriction  along  Tq  =  t,  Tx  =  Xk(t)  yields  the  composite  solution 


"y  (t)  =  Cj  id'1/4  exp(iXjw^dt)  +  ca  uj  _1A  exp(  -iX Juj  ^  dt) 


(3.2.  17) 


which  may  be  recognized  as  the  Liouville -Green  solution. 

This  approximation  will,  however,  break  down  when  the  neglected 

a8  y 

term,  viz.  —  %  becomes  of  order  X.  Substituting  (3. 2. 16  a)  in  (3.  2. 15  c): 
o  T 
O 


=  a  eTl  =  (oi  *1/4  (T  )  )"  eT 
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Clearly  the  failure  of  the  approximation  (3.2.  17)  i>  indicated  near  th<s  zero* 
of  d)  (t ). 

(b)  Second  Order  Equation;  Noncanonical  Form 

The  analysis  proceeds  in  the  same  fashion  as  before  and  different 
approximations  can  be  obtained.  In  thi'j  case,  however,  we  may  slightly  alter 
the  point  of  view  and  determine  a  class  of  equations  which  is  maximally 
informative  with  respect  to  a  given  extension*  We  consider  therefore  the 
equation 

m  .  n  _ 

y  +  e  ®i  y  +  s  r^y  =  0 

with  the  extension 

T0  =  t  ;  Tt  =  ck(t) 

and  wonder  what  values  of  m  and  n  would  correspond  to  the  maximal  or 
submaximal  balance  of  the  terms,  together  with  simplification.  This  is 
considered  in  the  limits  of  small  and  large  e  •  The  extended  equation  is; 


+  e  <k  IS-  + 

OTj 

(1) 

(2) 

(3) 

+  *  “  («o  x  )  + 

dTo 

em+1Uik  )  +  eD(  <roy)  =  0 

(3.2.18) 

(4) 

(5)  (6) 

The  terms  are  ordered  as  e  °  :  e  1  :  e*  :  e  m  :  c  m+*  :  e.  D  ::  1,  2,  3,  4,  5,  6. 
Using  the  graphical  technique  (Fig.  11c)  the  relations  for  maximal  ordering 
are  (i)  m  =  n  and  (ii)  m  +  1  =  n.  This  is  obtained  by  balancing  the  exponents 
of  c  taking  two  terms  at  a  time.  Each  balancing  defines  a  curve  in  the  m,  n 
plane.  The  number  of  terms  balanced  for  each  (m^  ,  )  is  given  by  the 

number  of  curves  passing  through  (m^  ,  nj  )  plus  1.  Maximum  balance 
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It  may  be  noted  that  by  means  of  the  transformation 


y(t)  =  z(t)  exp  (-  £  «it) 
the  noncanonical  equation 

y"  +H-1  (t)y'  +  flatty  =  0 


can  be  transformed  into  the  canonical  form 

a  n 

*” + 1  A.  -  -4-  -  -f- 1  *  ■  0 


(3.2.22) 


(3.2.23) 


(3.2.24) 


We  may  therefore  choose  to  study  approximations  for  the  noncanonical 

equation  directly  or  the  canonical  form  after  the  above  transformation. 

The  clock  functions  given  by  (3.2.7)  and  (3.2.20)  are  different  in  the  two 

cases.  However,  the  approximations  (3.2.8)  and  (3.2.21)  are 

unaffected  to  leading  order  after  taking  into  account  the  transformation 

(3.2.22).  The  difference  in  the  clock  functions  (3.2.7)  and  (3.2.20) 

in  the  case  when  il,  =  eui.  and  -0-  =  c  w  exactly  corresponds  to 

o  o 

the  noncanonical -canonical  transformation  (3.2.22). 

(ii)  Long  Time  Approximation:  (e  large).  The  alternative 
balancing  of  m+l=n;  m  =  1,  n  =  2  leads  to  a  different  approximation 
for  large  values  of  e  •  The  equation  now  is: 

y"  +  eu)j  y1  +  ?a  a>oy  =  0  (3.2.25) 

This  ordering  can  be  obtained  in  the  following  way  also.  Consider  the 
equation: 


+  u  o(Xr  )y  =  0 


(3.2.26) 
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The  coefficienta  are  slowly  varying.  If  X  <<  1.  With  the  transformation 
XT  =  t,  the  equation  can  be  written  as; 


dr 


^f"  +  e  (t)  +  e''ttiA(t)y  *  0  ;  e  =  >  >  1 


The  extension  T  =  t; 

0 

T  i  =  €  k(t) 

c a  :  ica 

t*  * 

If, 

e1  :  k 

il  +  2k 
ar. 

ili —  + 

8ToaT' 

iiL  = 

3T 


=  0 


(a) 


(b)  {3.2.27} 


ar 


Xr  o 


(c! 


The  coefficients  of  (a)  can  be  treated  as  constants  w.  r.t.  Tx  and  the  solution 
can  be  written  as 

tm 

y(To.r1  )  =  a  (T q)  e  1 

where: 

ka  +  ojx  (k)  +  u)q  =  0  (3.2.28) 

The  amplitude  variation  ir  obtained  from  (b)  by  substituting  and 
integrating  as: 

Qf(To)  =  (2k  +  wx  )Va  Wher*  o7o  (  fn  V  J  =  (41c  +  2UJJ  )  (3.2.27) 

• 

When  idj  is  a  constant,  or(T  )  =  (2k  +  u'i  )  *  .  However,  even  when  a'j  is  not  a 

0  l 

^ 

constant,  consider  the  function  a(TQ)  =(2k  +  u>l  )  and  the  approximation: 


y  (T  »TX  )  =  <*(T  )  exp(T i  ) 


(3.2.30) 


This  depends  on  whether  y(TQ)  is  a  slowly  varying  function.  In  any  case, 
(3.2.22)  can  be  used  as  an  approximation  if  the  error  estimates  are  known.  Such 
an  analysis  is  made  in  the  next  chapter.  This  approximation  affords  a  method  of 
shifting  the  point  at  which  the  standard  LG  approximation  breaks  down,  and 
will  be  discussed  in  later  chapters. 
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3. 3 'Third  Order  Equation 

Consider  the  equation: 

y1"  +  u>ay"  +  (!>!  y'  +  tw^y  =  0  (3.3.1) 

The  transformation  y  =  z  exp{-  J'oua  dt)  converts  the  above  equation  into  the 
canonical  form: 

z'"  +  z'  +  H  z  =  0  (3.3.2  a) 

o 

Let  us  therefore  consider  the  canonical  equation 

y,M  +  e  m  u)j  y1  +  e  n  y  =  0  (3.3.2  b) 

and  the  extension  r  =  t;  Ti  =  e  k(t).  We  determine  a  class  of  equations 

o 

for  which  this  extension  obtains  maximal  or  submaximal  balance  of  terms, 
in  the  light  of  the  discussion  preceding  (3.2. 14).  The  choices  are  found  to 
be  (i)  m  =  n  or  (ii)  m  +  1  =  n.  Let  (i)  m  =  n  =  1  and  (ii)  m  =  2,  n  =  3. 

(i)  Short  Time  Approximation:  (c  <<  1) 

The  equation  is  now  ordered  as  : 


y'"  +  e  u)i  y'  +  e  u>  y  =  0 
o 


(3.3.2  c) 


Using  the  extension  (3. 1.3): 


°  a3  y  n 

1  'TT*  -  0 

o 


(a) 


e  k  ^-  +  3k * —  +  3k  ~ \ +  WJ  =  o  (b) 
*TX  aTjdTj,  8Toaarx  *3To  O7 


(3.3.3) 


»:  3kk  +  3k8  pL-  3  +  .xk|l-  =0 

ati a  sr  at! 3  1  a rx 


if  ^  =o 


ar, 


(C) 

(d) 


For  small  e  the  pertinent  equations  are  (a)  and  (b)  above.  Solving 
them  in  this  order: 


y(T q.  rx  )  =  a(tx  )  r  oa  +  b(tx  )  tq  +  C(rx  ) 
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Since  terms  on  the  r.h.  8.  are  linearly  independent  w.  r.t.  T  (Wronskian  t  0), 

o 

each  can  be  used  to  generate  a  clock  function.  Substituting  y=A(Tj  )T  * 
in  equation  (3.3.3b): 


kj  A'  T  3  +6k.  A'r  +  6k.  A*  =  -  A(2o>.  T  +  •*»  T  3  ) 
o  K  *ooo 


Choosing  A  =  exp{T1  )  leads  to  the  second  order  equation  for  the  clock: 


T.4ik!  +  4T  akl  + 


.  -g—  +  ot  y -  » -  +  6(k-  )  =  -  (2u>.  T  +  U)  T  3  ) 

o  dr  o  dr  ^  1  o  o  o  . 


(3.3.4) 


Though  this  is  a  variable -coefficient  equation*  it  can  be  readily  solved 

being  recognized  as  the  inhomogeneous  equidimensional  or  Euler -Cauchy 

2 

equation.  The  transformation  T  =  e  reduces  Equation  (3.3.4)  to  a 
constant  coefficient  equation  in  z.  Alternatively  it  can  be  written  as 


da$  A  6  d«I>  6 

d +F  +  * 

o  o  o  o 

where  <6  =  1^  and  f  =  -  (  +  uj  )  . 

T  o 

o 


The  solutions  of  the  corresponding  homogeneous  equation  are  obtained  as 


(3.3.5) 


$  =  T  where 
o 


m(m-l)  +  6m  +6  =  0 


i.  e.  nr^  =  --3  ;  mg  =  -2 


and:  =  r  "*  ;  =  T 

*  o  3  o 
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The  particular  solution  is  given  by: 


$>  =  -$  l  (T  ) 


ds  +  $a  (t  )  f 

i.«.)  J 


%a  (s)f(s) 
W($ 


li liliiiL  d8 

W(*x  ,  $a  ) 


W($x.$a  )  = 


_3 


-3r 


-2r 


=  T 


-6 


♦  =  K  =  T, 


-3 


f 


r  a  (2a),  +  w  r  )dr  -  r 

o  o  o  o  o 


-2 


Ta  (2u).  +  T  r  )dT 
o  4  o  o  o 


Integrating  by  parts  and  noting  that  T  =  t 

ki  =  f  j  t3  ga  dt  “  fia  j  t3  ga  dt  -  J  tga  dt 


(3.3.6  a) 


where: 


ga  (t)  =  2U)I  +  u>  t 

-3-3  i 

For  example,  when  uux  =  0;  uj^=t  ,  g3  =  t  and  kj  =  constant  -jflnt. 

The  approximate  solution  is  given  as  : 

y  (t)  =  cta  exp(-  y  In  t)  =  c  t^S  ^  (3.3.6  b) 

The  exact  solution  for  this  example  can  be  obtained  as  y=  tm  where 
m  satisfies  the  equation 

m(m-l)(m-2)  +  e  =  0 

i.  e.  :  m3  -  3ma  +  2m  +  c  =  0 
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For  small  t  ,  expanding  m  =  +  em;  +  ...  and  taking  m0=  2  the 

correction  to  order  e  is  obtained  as: 


3m  3  m,  -  6m  m,  +  2m,  +1=0  or  rru  =  -& 
o  1  o  1  1 


Hence,  for  small  e  ,  the  approximate  solution  to  order  e  is 
(2  -e/a  ) 

y  =  c  t 


(3.3.6  c) 


which  is  indeed  predicted  by  the  approximation  via  time  scales  (3.  3.  6  b). 

The  clock  functions  corresponding  to  the  other  two  solutions  ot 
(3.3.2  c)  are  similarly  obtained. 


2-t  ^  (<i)  j  +  u)  Qt  )t  +  ^  (oc  x  +u;  Qt)  t  dt  ~  k  J  (<*>  \  +  tw  ^t)  dt  dt 

k3  =  -  JJ|„od'3 


and  k 


Using  these  clocks,  the  approximations  to  order  e  can  be  written  as 


yjo  =  exP(  "  f 

o  o 


where  ga  =  2^  +  uj  t;  =  uj^  +  uj  t. 

-3 

For  the  case  when  uul  =  0,  uo  =t  , 


5  t  j~J t3  ga  dt  -  |  t3  ga  dt  -fc  J tga  dt)  ) 

(a) 

l  4jt=gl  dt  +  jtgidt-  ||  gl  dt3  )  ) 

(3.3.7)  (b) 

ifk-> 

(c) 

_  A  8  -  e/a  ) .  (1  +e  )  .(-  e£  ) 

ya  =  c*  t  ;  Yl  =  cx  t  ;  Y0=Cot 

That  these  are  indeed  the  correct  approximations  to  order  e  can  be  seen  by 
making  the  expansion  m  =  mQ  +  e  m :  +  .  .  .  and  evaluating  mt  from  the 
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2 


relation  3m^2  mt  -  6m^m1  +  2mj  +1=0  and  then  successively  taking 

m  =  2,  1  and  0. 
o 

The  failure  of  the  approximation  can  be  studied  as  before  by  comparing 
terms  from  (c)  and  (b)  in  the  Equation  (3.3.3);  failure  occurs  when: 


(»„y) 


(  — 


1 


k 


lx_ 


1 

e 


Taking  kg 


u)  dta  ,  condition  for  breakdown  is 
o 


c 


u)  dt3 
o 


(3.3.8  a) 


When  =  C,  it  is  better  to  obtain  the  criterion  differently  as 


taking  the  terms  from  (3.3.3  b)  and  (3.3.3  c)  respectively; 
i.  e.  : 

3e  (jj\dt3  )(  Ju)odt  )  ~  u>o  (3.3.8  b) 

The  approximation  breaks  down  near  the  value  of  t  for  which  it  is  satisfied. 

When  tu  =  4  ,  we  see  from  (3.3.b  b)  that  the  approximation  fails  when: 
o  t 

3c  ~  t3  '  i,e‘  When  ~  1 

For  small  c  the  time  scales  approach  gives, in  this  case, the  correct 

approximation  to  the  leading  order  in  c  ,  and  fails  as  e  increases  in 

2 

magnitude  towards  unity.  The  value  e  =  may  be  verified  to  exactly 

correspono  to  the  occurrence  of  multiple  roots  for  m  in  looking  for  a  solution 

m  2 

of  the  type  y  =  t  .  The  solution  is  therefore  oscillatory  when  e>j7'j  *nd 


62 


nonoscillatory  when  e"  £  jyr^-  .  Thus  for  small  c  ,  the  method  yields 
useful  approximations. 

Higher  order  corrections  can  be  obtained  by  employing  slower  clocks 
as  in  the  second  order  equation  and  the  validity  of  the  approximation  scheme  can  be 
studied  as  before. 


(ii)  Long  Time  Approximation  (  e  large  ) 

The  alternative  balancing  in  Equation  (3.3.2  b)  with  m  =  2,  n  =  3, 
leads  to  the  equation 


y'"  +  e  u.'i  (t)y 1  +  c  ou  o<t)y  =  0 


which  v/ill  be  studied  for  large  values  of  c.  Using  the  extension 


t  T  »  T x  }  ;  T  =  t;  rx  =  ek(t) 
o  o 


the  extended  equations  for  large  c  ,  are  in  the  order 


k3  9 


— +  IDJ  k  +  w  Y  =  0 

I  1  0*1  o 


:  3kk  ill,  +  3k3  — * a  +  Wl 

arl a  ar  arx  3  1  ar 

o  4  o 


=  0 


:*!■*-+  3k  M 

atx 


aroaT: 


+  3k  — \ =  0 

^  m  •  ^  m 


3  To  3  T 1 


i-4- 
at J 


=  o 


As  before  these  can  be  solved  as 


y(T  ,  Tl  )  =  a(T  )  eT  1 


where: 


k3  +  ui  i  (k)  +  u)  =0 
1  o 


(3.3.9  a) 


(3.3.9  b) 


(a) 

(b) 

(3.3.  10) 

(c) 

(d) 


(3.3.  11) 

(3.3.  12) 
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Substitution  of  this  into  (3.3.10  b)  results  in 
V  (T  ) 

°^T  o}  =  •  - - * 

(3k3  +  u*1r 

where: 

Un  y)=^ -  =  *  4^T-  (4n  (3k3  +  «,  )  )  (3.3.13) 

dTo  2(3 k3  +  )  hTo 

restriction  along  T  =  t;  Tj  =  e  k(t)  leads  to  the  approximation  for  y. 

The  noncanonical  third  order  equation  can  be  similarly  studied, 
yielding 


(To)  = 


v3  (r  ) 

_ nc  ° 

r  3k3  +  2uja  (k)  + 


where 


(M  V 3  )  =  £  ~  (tn  (3k3  +  2uja  (k)  +  u>x  )  ) 
T  o  nc  o 

2ul* 3  (k)  4-  uo a _ 

2(3 ic3  +  2uja  (k)  +  .0x  ) 


Again  the  extended  function  y(T  ,  Tj  ),  must  be  restricted  along  T  q=  t; 
Tj  =  c  k(t)  to  obtain  the  approximation. 

The  third  order  noncanonical  equation  is  equivalent  to  the 
canonical  equation  of  the  form 


y"'  +  (e  a  o» x  +0  (e  )  )y'  +  (e  3  oi  q  +  0(c  3  )  ) y  =  0 

for  large  e  .  It  can  be  seen  that  to  leading  order,  the  frequency  of  the 

solution  i6  determined  by  u>l  ,  u)^;  however,  the  amplitude  is  affected  by 

ir,  ,  ir  and  also  by  the  terms  0(e  )  with  y*  andO(ea  )  with  y. 

*  o 
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3.4  The  Linear  Equation  of  O’der  n 


The  l.d.e.  of  fourth  and  higher  orders  can  be  analyzed  in  a  similar 
manner.  In  each  case,  after  choosing  the  two  time  scale  extension,  the  proper 
balance  of  terms  for  short  or  long  time  (or  alternatively  for  small  or  large  values 
of  the  parameter  e  )  can  be  determined  by  the  principle  of  maximal  or  submaximal 
balance.  In  the  fourth  order  case  (for  short  times)  for  example,  the  clock 
functions  are  found  to  satisfy  an  inhomogeneous  third  order  l.d.e.  with  variable 
coefficients,  but  again  of  a  particular  type,  namely  the  Euler-Cauchy  equation. 
Without  going  into  the  details  of  this,  however,  in  this  section  the  method  of  time 
scales  will  be  applied  to  the  general  n^  order  l.d.  e.  and  some  general  results 
obtained.  The  extension  of  the  n^  order  derivative  operator  is  derived  in 
Appendix  I.  Consider  the  equation 


Sjl  +  17  _ 

dtn  n-l  dtn 


dn_1z  n 

T  +  ,  ,  ,  +  ^  Z  =  0 

-  o 


This  can  be  transformed  into  the  canonical  form 
dn 


q  y  . 
dtn  +  « 


Jn_a 

d  y 

n-a  dF=^ 


where:  z  =  y  exp( 


,  rn 

n  J  n- 


+  .  . .  +  u)  y  =  0 
o 

dt) 


(3.4.1) 

(3.4.2) 

(3.4.3) 


(i)  Short  Time  Approximation:  (small  e  ) 

t 

Consider  the  parameterized  equation 


,n 


,n-a 


+  r( +  ••• +  %y)  =  ° 


(3.4.4) 


in  the  limit  as  c-0,  Direct  perturbation  theory  can  be  shown  to  be  nonuniform, 

depending  on  the  forms  of  the  u>  ,(t);  i=0,  1,  ( ,  n-2.  The  extension 

t  «=$>“(  T  , T i  ]  ,  T  =t;  Ti  =  f  k(t)  leads  to  a  set  of  (n+ 1 )  partial  differential 
o  o 

equations.  For  small  e  only  the  lowest  and  first  order  equations  are  retained  in 
this  analysis.  Using  the  result  derived  in  Appendix  I  the  extended  equations  can 
be  written  order  by  order  as: 


■  ■ 0 
o 


(a) 


Cl 


/  \  n_1  /  X  -vr  +  1  *  n'8 

k<n)  P-  +  z  ( n)  k(n'r)  3_1.  .  -u, 

brl  r  =  1  r  8T>i  n-9  ar 


(3.4.5) 


rfa  +•  •  *+a;oy 


(b) 
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plus  higher  order  terms. 

Integration  gives  : 

(3.4.6) 

y(T  ,r  j  )  =  A  . (r  i  )r  n"  +a  Atx  )t  n"  +  . . .  +  Ax  (^  )t  +  A  (Ti  ) 
o  n  - 1  o  n  -2  *  o  o  o  * 


At  this  point  we  make  use  of  the  linear  independence  of  A.(T  x  )T  ,  i=o,  1, 
and  generate  n  clock  functions,  each  corresponding  to  one  solution  of 
(3. 4.  4)  starting  with  -  ich  A.(T  j,  )T  *.  Considering  A^  ^(Tx  )T  *  , 
substitution  into  (3.4.5  b)  gives: 

,n,  ,(n-r), 

d  k  ,  ,  n-1  d  k  r  „  . 

dT-^  A'»-1(T‘  )To  +rZL1  <  r  1  ,7Tn—  A  n-1  ^  (To  > 

o  dT  o 


n-1 


n-2 


"(u,n-2  +  “*  +WoY  ]  (3.4.7) 


Choice  of  an  exponential  T  1  dependence  of  A^  ^  will  result  in  the  following 
equation  for  the  clock; 

_  n-1  dn  1  (k  . )  n 


°  dT"* 

o 


k  )  n  _  ,n-2  * 


dn"3(k  ) 

+  ("Kn-l)(n-d)ro“‘3  g^-irrf--  +...  +(n!)(kn.,) 

o 


=  -(u;  ,T  +  u)  ,  T  3  +  .  . .  +  ui  r  n_1  ) 
n-2  o  n-3  o  o  o 


This  can  be  written  as 
,n-l 


n-1 


^  n-1 


dT 


n-1 


+  a 


n 


-2  o 


n-2 

2  °  a  1 

7  n  - ^  +  . .  .  +  a  ip  =  f(T  )  (3.  4.  8) 

,  n-2  o  n-1  o 

dT 


where  ©.  =  k.  .  This  is  recognized  as  the  inhomogeneous  Euler -Cauchy  or 
equidimensional  equation.  The  solution  of  the  corresponding  homogeneous 
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m 


equation  can  be  expressed  as  ^  f  where  m  satisfies  the  algebraic 
equation 

m(m  -1  )(m  -2 )  ...  (m-n+2)  +  m(m -1  )(m  -2 )  ...  (m-n+3)(  ^  )(n-l ) 


+  . . .  +  n!  =0 


(3.4.9) 


which  can  be  written  as 


n  - 1  n  -2 

m  +a  ,  m  +  . .  .  +  a,  m  +  a  =  o 
n-2  1  o 


(3.4.  10) 


having  (n-1)  roots  which  will  be  assumed  to  be  distinct.  The  homogeneous 
solution  is  given  by 


m 


....  ma  n-1 

cp  J=  k  J=  c,  T  +  . .  .  +  c  r 
n,- 1  n-r  *  o  n-1  o 

nom  horn 


(3.4.  11) 


which  will  be  denoted  respectively  by  ,  oo#,  .  .  .  0  .  The  particular 

n-1 

solution  can  be  written  as  (Ref.  ) 


V,  f  , 

wT)  dT»  '  , 

O  J 


JL  dT 

W(To)  To 


n 


+  .  .  .  +  (  -1 )  cp 


n-1 


fv  f  dr 

n-1  o 

w(r  ) 

o 


where  W=  W(^1  ,  cpa  »  •  •  •  >  cA  , )  is  the  Wronskian  and  V.  is  the  determinant  of 

th 

the  matrix  obtained  by  replacing  the  i  column  of  the  (n-1)  square  matrix 


Cpl 

Ta 

*n-l 

, _ 

~ . 

cp  'l 

cp  '8 

’  n-1 

Cp  "l 

~  (n-2) 

cpi 


~  (n-2) 

cpa 


CP 


(n-2) 

n-1 
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by  the  column: 

0 

0 

f 


One  more  integration  of  cp^  ^  gives  Thus  obtaining  the  proper  time 

scales  and  after  restriction  along  TQ=t,  T  j  =  e  k^  ^(t),  one  approximate 
solution  of  (3.4.4)  to  order  e  can  be  written  as: 


,  .  n-1 

y(t)  =  c  t  exp(  e 
n  -  i 


%-i at) 


The  other  independent  approximations  are  similarly  obtained  by  determining 

the  clock  functions  cp  i  .  o  .  and  combining  each  with 

n -2  n-3  n-*  o 

t  ,  t  ,  . . . ,  t,  t°  as  in  (3.4.6).  Therefore,  the  approximate  general 
solution  of  (3.4.4)  to  order  e  (for  small  c  )  can  be  written  as 


y(t)=  c  tn_1exp(e  [cp  ,dt)+c  tn"2exp(c  f cp  ,dt) 
n-1  J  n-1  n-c  J  n  -»• 

+  . . .  +  Cjt  exp(e  Jcpi  dt)  +  cq  exp(e  K  dt'  (3.4.13) 

where  c .  are  constants  and  cp.(i  =  0,  1,  ....  n -1 )  are  clock  functions. 
i  i 

For  any  given  n,  the  breakdown  of  this  approximation  may  be 
investigated  as  shown  in  the  last  section  by  studying  when  the  terms  which 
are  neglected,  viz.  of  order  ea  ,  become  of  the  same  order  as  terms  of  order  e. 

Higher  order  corrections  can  again  be  obtained  by  employing  slower 
clocks,  of  order  ca.e3.  etc. 

(ii)  Long  Time  Approximation:  («  large). 

The  ordering  of  the  parameterized  l.d.e.  of  order  n  can  be  obtained 
by  requiring  all  the  coefficients  to  be  slowly  varying.  For  example,  the 
equation 
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dn 

dr 


,n-l 


n  +U)n  ~  +•••  +  U>1  (XT  )  ^-  +  (1)  (Xt)y  = 

n  n-l  n-1  dT  o 

df 


(X  <<1) 


can  be  transformed  into 


(3.4.14) 


d 

dt 


+  e  n_i'^  "~n-*l~  +  •  •  •  +  e  n  1ul'  (t)  +  C  Q(t)y  =  0  (3.  4.  15) 

dt 


where  t  =  XT  and  e  *  r  >>  1.  This  equation  can  now  be  studied  in  the  limit 
of  large  e  .  To  use  time  scales,  the  domain  of  the  independent  variable  is 
extended  as  beforer 


t  T  .Tj  }  ;  T  =  t,  Tx  =  c  k(t) 

o  o 


(3.1.3) 


Equating  powers  of  c  will  lead  to  a  ret  of  (n+1)  partial  differential 
equations.  For  the  limit  of  large  %  only  equations  retaining  terms  of 
order  e"  and  e**  *  need  be  considered.  Using  the  results  derived  in 
Appendix  I  these  equations  can  be  written  as: 

<k)n|^n+‘"n.1(k)n'lL~J_1  +...  +u)x(k)|^  +^oy  =  °  (3.4.16) 

S  T  i 

n(i°n'1  fe-T-V  +<"-»>»„.! W-2  +(n-*kn>,n-3^ 

o  O  dT  dTx 


+ . . .  +  2o)  q  k  y  —  +  uj  .  —  y 


bTdT: 


ST 


n 


-2 


I  n(n_1>  fUn_2ir‘  1  n'ly  ■  (n-l)(n-2)  •  n-3-  5  y 

+  ~ — <k>  k  r,-r.i  + - 5 - “,n-i<k)  kFf7^-a 


arl 


+  ..  .  +  u>8  k  =0  (3.4.  17) 

Equation  (3.4.  16)  has  coefficients  independent  of  and  we  seek  solutions 


n-3 
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of  the  form  y{T  ,  T,  )  =  a  (f  )  0(T,  )  =  a  (r  )  exp(T,  ).  The  clock  function 
0**0  o  * 

satisfies  the  equation 

(k)+u>  (k)  +...+u>x(k)  +  ii)  ~0  (3.4.18) 

n-i  o 

the  roots  of  which  are  taken  to  be  distinct  for  this  analysis.  When  any  of 
the  roots  coalesce,  the  approximation  (3.4.21)  fails  because  the  amplitude 
factor  becomes  unbounded.  In  this  case  a  nonelementary  function  such  as  the 
Airy  function  is  necessary  in  order  to  represent  the  true  solution.  To  relate 
the  approximations  via  elementary  functionron  either  side  of  a  point  where 
the  roots  coalesce,  one  is  faced  with  a  nontrivial  connection  problem.  For 
distinct  roots,  the  explicit  amplitude  variation  is  obtained  from  (3.4.17)  as: 

,  n(n-l )  w;  ,n-2  (n-l)(n-2)  4n-3  ,  ,  „.t 

d  <  — 2 —  *  +  - 2 -  “n-1  {k*  +  •  •  •  +  ®a  )  (k) 

—  (In  a)  = - - - ; - = - ; - 

To  (n  (k)n  >(n-l)i«n  ^  (k)D  +...  +  2oj  a  (k)  +  u>  j  ) 

=  '^llk)  (k)n'Z  +  ...i«,)}Il  (3.4.19) 

This  can  be  written  in  compact  form  as  follows.  Given  the  l.d.e.  (3.4.15), 
consider  the  corresponding  "characteristic  expression": 

F(x,t)  =  x°  +  u)  xn  +  . . .  +  m  x  +  u)  (3.4.20), 


F  =  0  is  defined  to  be  the  "characteristic  equation.  "  The  explicit  amplitude 
variation  can  be  determined  from  the  relation: 


5t  (tn  a)  =  -*  h  l°  1 


a  F  .  dx 
ax  dto 


Integration  gives 


where  y (T  )  satisfies  the  equation: 
o 


d_ 
d  T 

o 


in  <  Tx1 


* 


Combining  the  Tq  and  T*  behaviors  and  imposing  the  restrictions  To=t 
and  Tj  =  e  k(t)  the  result  may  be  stated  as  follows. 

Given  an  l.d.e.  (3.4. 15)  the  approximation  via  time  scales  for  large 
e  is  obtained  as 


n 

E 

i=l 


(a) 


where  y  .  is  given  by 


Sr  WB  V  *  ‘  ‘  Hr’* »  ib)  (3.4.2D 

i 

where  F  is  given  by  (3.4.20)  and  x.  are  the  roots  of  F  =  0. 

Clearly  when  the  u>.  ,  (i=l,  2, . . .  ,  n-1)  approach  constants  as  t-*<», 

y  approaches  a  constant.  Also  if  w ^  is  the  only  varying  coefficient  y  is  a 
pure  constant.  In  this  case  the  approximation  is  given  by 

y(t)  =  E  c  (-|-^— )"^  «xP(e  f  x:dt)  (<0 

i=l  id*.  J  1 

where  c.  are  arbitrary  constants.  This  can  also  be  interpreted  to  mean  that 

we  impose  the  condition  that  <11,  ,  a>«  .  .  .  tv  vary  more  slowly  than  uu  . 

1  n  -l  o 

It  can  easily  be  verified  that  the  formula  (3.4,21)  recovers,  for  second 
order  equations  the  results  cf  the  standard  Liouville -Green  (WKBJ)  theory  as 
well  as  those  obtained  by  Curtiss  (Ref.  27). 

3.5  Summary  of  the  Chapter 

The  main  theme  in  this  chapter  has  been  the  demonstration  that  when 
the  direct  perturbation  theory  breaks  down,  a  natural  time  scale  can  be  found 
on  which  the  solution  can  be  described  uniformly.  The  method  of  time  scales 
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is  shown  to  systematically  lead  to  the  determination  of  the  natural  "clocks" 
of  the  problem. 

Linear  differential  equations  are  studied  in  the  light  of  multiple 
time  scales,  beginning  with  the  first  order  equation.  Reparameterizing 
the  equations  enables  separation  of  the  different  behaviors  and  a  clear 
physical  picture  is  shown  to  emerge.  In  each  case  failure  of  the  classical 
perturbation  approach  is  examined  and  the  time  scales  method  is  shown  to 
improve  on  the  approximations. 

After  deriving  the  exact  solution  of  the  first  order  equation,  the 
second  and  higher  order  equations  are  studied  in  the  two  limits  as  a 

parameter  e.  becomes  small  or  large.  The  ordering  in  each  case  is  justified  by  max¬ 
imal  or  submaximal  balance  of  the  terms.  Higher  order  corrections  are 
derived  for  the  second  order  equation  and  a  criterion  of  validity  of  the 
approximation  is  derived.  In  the  large  e  limit  the  standard  Liouville -Green 
(WKBJ)  approximation  is  derived  via  time  scales.  The  noncanonical  formulation 
is  shown  to  lead  to  a  different  approximation;  the  significance  of  this  new 
approximation  becomes  apparent  in  later  chapters  where  it  is  shown  to  be 
useful  in  the  transition  point  analysis. 

The  theory  is  then  generalized  to  the  order  equation.  For  the 
small  c  behavior  the  clock  functions  are  shown  to  satisfy  the  Euler-Cauchy 
equation  and  hence  can  be  determined  exactly.  For  the  large  c  limit,  a 
compact  formula  is  derived  which  enables  one  to  write  down  the  approximation 
by  inspection.  It  is  observed  that  the  time  scales  are,  in  general, 
nonlinear  functions  in  addition  to  being. complex  quantities. 
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CHAPTER  IV 


ERROR  ANALYSIS 

The  aim  of  this  chapter  is  to  examine  the  validity  of  the  approximations 
derived  in  the  last  chapter  and  to  obtain  bounds  on  the  errors  incurred. 

The  first  order  equation  presents  no  problem,  since  the  solution  obtained 
from  the  proper  choice  of  time  scales  is  exact.  For  higher  order  equations 
one  has  necessarily  to  resort  to  approximations  and  an  estimate  of  the 
errors  is  sine  qua  non  for  the  analysis.  Though  there  is  an  extensive 
literature  establishing  precise  conditions  for  the  existence  of  asymptotic 
solutions,  strict  upper  bounds  for  the  errors  have  not  in  general  been 
formulated.  Blumenthal  (given  in  Ref.  28)  did  obtain  uuch  bounds  for  the 
second  order  equation  as  early  as  1912,  but  his  results  have  not  been  generally 
known.  More  recently  Olver  (Ref.  28)  showed  that  it  is  possible  to  deduce  from 
the  existence  proofs  sharp  upper  bounds  for  the  errors  instead  of 
the  O-symbols .  He  derived  these  error  bounds  for  the  Liouville -Green 
approximations  (or  WKBJ  functions)  and  their  derivatives,  and  showed  that 
these  are  indeed  both  realistic  and  easy  to  evaluate.  But  it  must  be 
remembered  that  these  bounds  are  valid  only  in  certain  regions  of  the  complex 
plane  which  are  free  from  transition  points  and  hence  allow  the  use  of  one 
and  the  same  form  of  asymptotic  expansion. 

In  this  section  a  few  preliminary  results  are  quoted  and  these  are 
used  in  subsequent  sections.  For  the  second  order  equation  a  new 
approximation  theorem  is  proved  based  on  Olver's  results,  and  this  is  seen 
to  deal  with  the  noncanonical  equation  directly.  This  will  be  further  exploited 
in  the  next  chapter  for  transition  point  analysis.  Since  the  derivation  of 
Olver's  results  is  quite  involved,  the  extension  of  his  technique  to  higher 
order  equations  is  not  readily  apparent.  Error  bounds  of  a  similar  type  are 
derived  in  this  chapter  using  a  direct  approach  starting  with  the  exact  solution 
and  this  is  extended  to  higher  order  equations.  The  discussion  is  made 
considerably  simpler  than  Olver's,  at  the  cost  of  imposing  restrictions  on 
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the  coefficients  of  the  differential  equations. 

4.  1  Some  Basic  Definitions  and  Useful  Lemmas 

Some  useful  results  will  now  be  stated  in  order  that  they  may  be 
used  later. 

Nonoscillatory  equations  and  oscillation  criteria  (Ref.  23). 

A  homogeneous  second  order  1.  d.  e. ,  with  real  coefficients  defined 
on  an  interval  J  is  said  to  be  oscillatory  on  J  if  one  (and/  or  every) 
real  solution  (j£  0)  has  infinitely  many  zeros  on  J.  Conversely,  when  every 
solution  (/  0)  has  at  most  a  finite  number  of  zeros  on  J,  it  is  said  to  be 
nonoscillatory  on  J.  Further,  if  in  addition  every  solution  (/ 0)  has  at  most 
one  zero  on  J,  the  equation  is  said  to  be  disconjugate  on  J. 

The  oscillation  theorems  of  Sturm  can  be  stated  in  many  ways.  For  the 
present  purpose  the  comparison  theorem  can  be  stated  as  follows  (Ref.  42), 

Sturm's  Comparison  Theorem.  Let  f(x)  and  g(x)  be  nontrivial 
solutions  of  the  1.  d.  e. 

u"  +  p(x)u  =  0  and  v"  +  q(x)v  =  0 

respectively,  where  p(x)  >  q(x).  Then  f(x)  vanishes  at  least  once  between 
any  two  zeros  of  g(x),  unless  p(x)  =  q(x)  and  f  is  a  constant  multiple  of  g. 

Corollary:  If  q(x)  s  0,  then  no  nontrivial  solution  of  the  1.  d.e. 
u"  +  q(x)u  =  0  can  have  more  than  one  zero,  i.  e.  "q(x)  £  0  on  J"  is 
sufficient  for  the  1. d.e.  to  be  disconjugate  on  J. 

The  proof  is  by  contradiction.  By  the  Sturm  comparison  theorem,  the 
solution  v  =  1  of  the  1. d.e.  v  "=  0  would  have  to  vanish  at  least  once 
between  any  two  zeros  of  any  nontrivial  solution  of  the  1. d.e. 
u"  +  q(x)u  =  0. 

One  of  the  very  useful  results  is  a  lemma  essentially  due  to  Gronwall  (Ref.  46). 

Gronwall’s  Lemma.  Let  X  (t)  be  a  real  continuous  function  and  jj  (t) 
a  non -negative  continuous  function  on  the  interval  (a,b).  If  a  continuous 
function  y(t)  has  the  property  that 
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.t 

y(t)  i  X(t)  +  J  u(s)y(a)ds 
a 

for  a  ;»  t  £  b  the  a  on  the  same  interval : 

t  t 

y(t)^X(t)+  J  X  ( s )  u  ( 8 )  exp(  j*  u(r)dT)ds 
a  a 

In  particular  if  X  (t )  =  X,  a  constant : 

t 

y(t)  *  X  exp(  J  u  (s)  ds) 
a 


Next  we  state  two  results  of  F.W.  J.  Oiver  (Ref.  28)  in  connection  with  second 
order  1.  d.  e. 

Theorem  1  (Oiver).  Let  u  be  a  positive  parameter,  and  f(u, x)  be  a 
continuous  real  or  complex  function  of  x  in  the  interval  a  £  x  £  b.  Then  in 
this  interval  the  differential  equation 


da  w 


f  u3  +  f(u,x)}  w 


has  solutions  wt  (u,x),  wa  (u,  x),  such  that 


wl  (u.x)  =  e^l+ej  (u,  x)  ),  wa  (u,  x)  -  e  '“(l  +  Ca  (u,x)  ) 

^■wl  (u.x)  =  u  e^l+Znj  (u.x)  ),  ~  w8  (u,  x)  =  -u  e  "“*(  1 a  (u,x)  ) 

where  |ej  (u,x)J  ,  |T11  (u,x)|  £  exp  {  ^  "  1 


ra(u,x)  ,  na(u.x)|i  exp{  )  -1 


b 

and 

r  i  1 

Fj  (u.x)  =  j  f(u,t)  jdt, 

a 

Fa(u.x)=J  |  f (u.  t )  |  dt 

The  interval  (a, b)  may  be  infinite  provided  that  the  integrals  converge. 
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Theorem  2  (Qlver).  With  the  conditions  of  Theorem  2  the 
differential  equation 

=  (  -u3  +  f(u,x)  ]  w 

iux 

has  solutions  w:  (u,x),  w3  (u,  x)  such  that  wx  (u,x)  =  e  +  £  \  (u,x), 

/  i  iux  .  ,  d  ,  .  iux  ,  . 

w2  (u,x)  =  e  +  e2  (u,x),  —  w:  (u,x)  =  lue  +  ur^  (u,x), 

d  ,  .  iux 

—  wa  (u,  x)  =  lue  +  uria  (u,x)  where 

J  cx  (u,x)  |  ,  j'n  (u,x)  J  ,  |  e  a  (u,  x)  ,  |r)3  (u,x)|  s  exp{  )  -1 

and  x 

F(u,  x)  =  J  J  f(u,t)J  dt 
c 

c  being  an  arbitrary  point  such  that  a  £  c  £  b.  The  interval  (a,b)  and  the 
value  of  c  may  be  infinite  provided  that  the  integral  converges. 

The  following  lemma  on  integral  equations  also  provej useful  (Ref.  43). 

The  Fredholm  integral  equation  is  written  as 

b 

■» 

cp(x)-A  k(x,  s)  cp  (s)  ds  =  f(x) 

a 

If  the  kernel  k(x,  s)  is  identically  zero  whrn|  s>x  (which  is  true  of  causal 

dynamic  systems),  the  integrand  is  zero  when  x  <  s  £  b  and  the  integral  becomes 
x 

J  k(x,  s )  cp  (s)  ds.  This  leads  to  the  Volterra  equation 

a 

x 

Cp  (x)  -  A  J  k(x,  s)  cp  (s)  ds  =  f(x) 
a 

Lemma.  (Ref.  43).  If  the  "free  term"  f(x)  in  the  Volterra  equation  is 
absolutely  integrable  and  the  kernel  is  bounded,  then  successive  approximations 
for  this  equation  converge  for  all  values  of  A. 

The  sequence  of  successive  approximations  is  given  by 
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X 

cd  (x)  =  f(x)  +  X  J  k(x,  a)  cp(s)ds 
a 

®  X  X  X 

+  £  X  f  k{x,  a,  )  )  k(s.  ,  sa  )  .  .  .  I  k(8  ,  a  )  f(8  )  ds  ...  da 

J  J  Jm-lmmm 

m=s  a  a  a 

If  Jk(x,s)|  £  andjf(x)|  £  Ma  ,  then  it  can  be  proved  (Ref.  39) 

by  induction  that  the  modulua  of  the  general  term  in  the  aerie8  for  rp(x) 
does  not  exceed 

IXfV  m  Ma  (x-a)m  /  (m!)  s  jX|m  m  Ma  (b-a)m  j  (m!) 

The  aeriea  convergea  uniformly  for  all  valuea  of  X. 

4.2  Approximation  Theorems  for  Second  Order  Equations 

We  will  now  prove  approximation  theorems  for  the  second  order 
non-canonical  l.d.e.  using  the  results  obtained  by  Olver. 

Consider  the  equation 

y"  +  eu)1y'  +  e?ujy  =  0  (4.2.1) 

o 

valid  in  an  interval  (a  £  t  £  b).  In  the  light  of  the  time  scales  treatment 
for  large  e  ,  the  characteristic  equation  is  : 

xa  +  u) x  (t)x  +  u  Q(t)  =  0  (4.2.2) 

The  time  scales  approximation  fails  when  the  independent  variable,  t,  has 
a  value  for  which  the  roots  coalesce,  i.  e.  when  the  discriminant  vanishes. 
Excluding  this  we  have  two  cases  for  distinct  roots,  viz.  when  the  roots 
are  real  or  complex  conjugates.  Each  case  will  be  discussed  separately. 

( 1 )  Case  of  Real  and  Distinct  Roots;  D(e,t)Huu,3  -4ir  >0 

1  o 

In  this  case  the  approximation  via  time  scales  is  given  as 
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**(*)  =  C\  { 
and 

Vb  *(t)  =  ca  f 


/  f  dt 

,xp(JT^~ 


2D 


T^- -  }  exp(  -  |- Ju^  dt)  exp  (  j  J  dt)  (a) 


exp{ 


°  -  )  exp(-|-  Ju)x  dt)  exp(  -  j  J^dt)  (b) 


(4.2.3) 


D1/4 


where  Cj  ,  ca  are  arbitrary  constants.  Let  us,  however,  consider  the 
following  approximations  which  are  valid  when  0  (Ref.  Eq.  (3.4.21  c)  ). 


Yi  (0  =-^ —  exp(  -  j  f  ui!  dt)  exp(  +  ~  fD^dt) 
Di/*  *  J  £  j 

ya(t)=~'4  e*P(“|-Ju>  dt)  exp(-  j^D^dt) 


(a) 

(4.2.4) 

(b) 


New  variables  £  and  W  are  introduced  and  defined  by  the  relation: 


y  =  mW 


where  m  and  n  are  as  yet  undetermined  functions  of  t  and  e  .  We  choose  them 
such  that  £  has  a  one-to-one  correspondance  with  t.  Let  a  and  0 
be  the  values  of  £  corresponding  respectively  to  a,b.  The  differential 
equation  (4.2.  1)  is  now  transformed  into 

mna  - — ^C-*  ^  +  (e  rnn  +  m'n  +  (mn)'  )  —■  W(e  ,  P  ) 
dt»  1  d? 

+  (e  3  u)  m  +  euilm'+m")W(c,£)  =  0  (4. 2. 5) 

o 


The  primes  denote  differentiation  w.  r.t.  t.  We  now  seek  mapping  functions 
rr.  and  n  such  that  in  (4.2.5) 

(i)  the  coefficient  of  the  first  derivative  w.r.t.  £  vanishes 

(ii)  the  coefficient  of  the  second  derivative  and  that  of  the  ea 
term  multiplying  W(e  ,  £  )  are  the  same. 
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The  first  condition  above  gives  us  the  differential  equation: 

(mn)'  +  m'n  +  cuj^  mn  =  0  (4.2.6) 

This  can  be  readily  integrated,  using  the  integrating  factor  m.  Thus^ 

™  f  (m2  n)}  =  -  e  u  x 
and  hence 

ma  n  =  c  exp(  -  p  J  ioj  dt)  (4.  2.  7) 

The  second  condition  gives 

m3  n  =  to  m  +  coefft  of  e  in  (a1,  m')  +  coefft  of  ea  in  (m")  (4.2.8) 

°  h 

(a)  If  m  is  independent  of  c  ,  (4.2.  8)  gives  n  =  cu  2  .  But  from 

i  o 

(4.2.7)  m3  =  io  exp(  -e  dt),  which  results  in  a  contradiction  unless 

°  J  -1/4  1 

u),  ■  0.  In  this  case  m  =  <*>  '  ,  n  =  co  ,  and  this  leads  to  the  standard 

o  o 

LG  result, 

(b)  If  m  is  also  allowed  to  depend  on  e  ,  the  choice  of 

exp(  -  f-  Ju)i  dt) 

m  - - 

D1/4 

leads  to 


which  satisfy  both  the  conditions  (4.  2.  7)  and  (4.  2.  8).  The  transformed 
equation  for  W  now  is: 

^  =  f  C  3  +  2C  -  4D"^4  (D  V  )  W 

d£ 3  D 

This  can  be  written  in  the  form 

—  =  (  ea  +  «c.t)  }  W 
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where: 


He  ,t)  = 


-  4D"^4  (D'^4)" 


(4.2.  10) 


We  now  suppose  that  f  is  a  continuous  function  of  £.  This  is  true  if 
D  is  twice  differentiable  and  does  not  vanish  within  the  interval.  Applying 
Olver's  first  theorem  a  solution  (z  ,  £,)  exists  such  that 

Wx  (e  ,f)  =  e€*(l  +  Ex  ).  ^  Wx  (e,£)  =  e  e£  f  (1  +  2ni  ) 


where 


|  Ex  |  ,  |nx|  *  exp  (  j^)  -  1;  Fx  =  J  |  f(e.A)  j 


-1/4  D7 

In  the  original  variables,  yx  =  D  and  d£  =  — ^ —  dt  ,  and  similarly  for 

the  second  solution  ya  . 

The  case  when  the  characteristic  roots  are  complex  conjugates 
(discriminant  ujj  3  -  4t ^  <  0)  can  be  similarly  treated  using  Olver's 
second  theorem.  Therefore,  the  following  two  approximation  theorems  have 
been  established  as  extensions  of  Olver's  results  for  the  standard  LG 
approximation. 

Theorem  3.  The  differential  equation 


y"  +  fuij  (e  ,t)y' +  ea  tu0(c  ,t)y  =  0 


(4.2.1) 


has  solutions  yx  and  ya  such  that 


y*  (e  ,t)  £  yx  (r,t)  exp  (  —  J  J  f (c  , t ) j  dt  ) 

a 

~  1  b 
ya(f.t)s  ya(t,t)exp(  +  —  J  |f(e,t)/dt 


(4.2.11) 


~  ^  F  f 

where:  p  is  a  positive  parameter;  y.  =(- — )  exp(  e  x.dt);  i=l,  2 

l  o  x  J  l 

and  are  the  distinct  real  roots  of  the  characteristic  equation  F  =  0,  with 


F  s  x3+  on  (e.t)x  +  w 
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and 


f(e  ,  t)  =  C 


■e  U'. 

D 


-  4D 


m-V 


{D  '  )"  ] 


?!  =[ 


-  2d“^4(D’^4  )"3 


D  being  the  discriminant  of  F,  is  positive  and  is  assumed  to  be  twice 
differentiable  in  the  interval.  The  interval  (a,b)  may  be  infinite  provided 
that  the  integrals  converge. 

When  the  characteristic  roots  are  complex,  the  following  theorem 
expresses  the  error  of  the  time  scales  approximation. 

Theorem  4.  The  differential  equation  (4.2,  1)  has  conjugate  solutions 


y  and  y*  such  that 

y  =  y  +  E 

(4,2.  13) 

where  : 

y  =  (  ^  exp(e  J  x  dt) 

(4.2.14) 

exp(-|-  [a)X  dt)  t 

|E|*  -l/4  (exp(  JJ  |f(e,t)j  dt  J  )  -1 

(4.2.15) 

a  <  c  <  b 


D,  F  and  f  are  as  given  in  Theorem  1  ,  D  being  understood  to  be  the  absolute 

value  of  the  discriminant;  and  x  is  the  complex  root  of  the  equation  F  =  0. 

We  see  from  theorems  3  and  4  that  the  error  of  the  approximations 

as  stated  is  0(1)  as  c  -ay.  However,  if  we  impose  the  condition  that  u>x 

varies  more  slowly  than  w  ,  then  the  error  is  o(l)  as  e  becomes 

o 

increasingly  large.  Thus  as  e  -<*> , 

(JL'X  =  IDj  (  — ) 

*  r 

.  1  - 

and  u),  =  —  uj, 

1  e 

In  this  case,  for  any  fixed  t,  the  error  of  the  approximation  becomes 
vanishingly  small  as  . 

Each  of  these  theorems  is  valid  on  one  side  of  a  transition  point 
(where  D  =  0)  since  they  are  corollaries  of  Olver's  results.  Nevertheless 
the  difference  in  the  form  of  the  approximate  solutions  can  be  used  to 
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advantage  when  dealing  with  the  case  of  multiple  characteristic  roots. 
A  brief  look  at  the  two  approximations  may  be  in  order  here. 

In  regard  to  (4.2.  1)  the  usual  LG  approximation  is  given  as: 


Yi  (0 


exp(-^-  dt) 


<%- 


(U  • 


dt  ) 


+  ca 


exp(  -ie 


>*dt  )  } 


However,  by  treating  the  noncanonical  form  of  the  equation  directly  by 
time  scales,  one  obtains,  under  certain  restrictions,  the  approximation: 


)^dt)  +  cs  exp(-ie  j(u)c 


} 


The  errors  of  the  approximation  have  been  given  in  each  case.  Yet  another 
approximation  can  be  written  as: 


fr(t)  = 


(w 


o 


{ 


c.  exp 


)^dt) 


+  c3 


)^dt)  } 


The  error  of  this  approximation  cannot  be  readily  written  by  application  of 
the  Olver  theorems.  It  can,  however,  be  estimated  by  a  different  method, 
which  will  be  the  subject  of  the  next  discussion. 
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4.3  Derivation  of  Erfor  Bounds 

In  order  to  assess  the  merits  of  the  time  scales  approximation  the 
estimates  of  the  errors  will  be  examined  using  a  well-known  method 
of  successive  approximations.  This  consists  in  treating  an  initial 
value  problem  and  writing  an  integral  equation  which  is  satisfied  by  the 
unknown  function.  Iteration  then  gives  a  sequence  of  successive  approximations. 
Bounds  of  the  Olver  type  will  be  rederived,  though  they  are  not  as  sharp 
and  the  conditions  are  more  restrictive  .  The  emphasis,  however,  is 
on  the  directness  and  simplicity  of  the  method  and  applicability  to 
equations  of  second  and  higher  order. 

Second  Order  Equation 

(i)  Self-Adjoint  Form.  Consider  again  the  equation 

y"  +  c3u;y  =  0  (4.3.1) 


in  an  interval  (a,b).  The  approximation  via  extension  is 

y  (t)  =  a)  ^  exp(±  iejui^dt) 
and  is  found  to  satisfy  exactly  the  equation 

y  "  +  (c  3  uj  +  f(t)  )  y  =0 

.  j3 

where:  f(t)  =  -a^4  -jp  (w  ) 

Equation  (4.3.  1 )  is  therefore  written  as 

y"  +  (c  3  a'  +  f)  y  =  fv 


(4.3.2) 


U) 

(4.3.3) 

(b) 


(4.3.4) 


with  f  as  given  by  (4.3.3  b). 
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If  f  is  small  when  compared  to  ea  a',  then  (4.3.  1)  and  (4.3.3  a) 
will  be  nearly  the  same  and  y  can  be  expected  to  be  a  good  approximation 
to  y.  In  order  to  estimate  the  difference  y  -  y" ,  (4.3.4)  is  regarded  as  an 
inhomogeneous  equation  with  the  r.h.  s.  as  a  known  forcing  function. 

The  method  of  variation  of  parameters  enables  one  to  write 


y(t)  =  y(t)+J  h(t,  s)  f(s)  y(s)  ds  (4.3.5) 

€ 


where  £  is  some  fixed  point  in  the  interval;  y  (t)  is  the  inhomogeneous 

solution  and  h(t,s)  is  the  Green's  function  (or  the  time -varying  impulse 

response)  respectively  of  (4.3.4).  For  a  fixed  s  in  (a,b),  h(t,s)  is  a 

function  of  t,  and  satisfies  (4.  3.  3  a)  together  with  the  initial 

conditions  h(sf  s )  =  0;  —  (sfs)  =  W(y  a  ,yl  )  =  1.  In  order  to  find 

3  t 

approximations  to  y  characterized  by  conditions  at  an  interior  point  of 
(a,b),  or  by  its  behavior  as  t~*a  o  -  t-*b,  we  shall  follow  the  argument 
presented  by  Erdelyi  (Ref.  44)  and  thus  £  is  chosen  to  be  the  point  in 
question  and  y  (t)  to  be  that  solution  of  (4.3.3)  which  is  characterized  by 
the  same  conditions  as  y(t). 

(4.3.5)  can  be  written  as  a  Volterra  integral  equation 

t 

y(t)  =  y  (t)  +  |  k(t,  s)  y(s  )  ds 
c 


where  k(t,  s)  =  h(t,  s)  f(s). 

From  the  theory  of  integral  equations,  the  sequence  of  successive 
approximations  converges  uniformly  when  y  and  k(t,s)  are  bounded  (Ref.  39,43). 
Let  the  equation  (4.3.5)  be  rewritten  as 


y(t)  =  AyA  (t)  +  Bya  (t)  +  y  x  (t)  J 


ya  (s)  f ( s )  y (s ) 
W(ya  ,yx  ) 


-ya  (t)  j 


y  i  (s)  f ( s  1  y(s) 

w(?a  ) 


(4.3.6) 
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where  y  j  ,  y2  are  the  two  independent  solutions  of  (4.3.3)  and  W(ya  ,  yx  ) 
is  the  Wronskian  defined  by: 


Ya  Yi 

Ya  Yi 


(and  is  equal  to  1);  A,B  are  constants.  Let  yx  (t)  and  ya  (t)  be  defined 
respectively  by  the  pair  of  conditions  A  =  l,  B  =  0  and  A  =  0,  B  =  l.  Also  let 
the  subscripts  1,2  respectively  correspond  to  the  +  and  -  signs  in 
the  exponent  in  (4.3.2).  The  definiteness  of  the  fractional  powers  of 
oj(t)  is  assured  by  agreeing  to  consider  positive  quantities  only. 

Case  (1)  Nonoscillatory  Case,  uj  (t)  £  0  by  Sturm's 
comparison  theorem  and  the  characteristic  roots  are  real.  For  distinct 
roots  we  require  uu(t)  <  0  in  (a,b).  The  conditions  under  which  fy  is  of 
constant  sign  in  (a,b)  can  be  determined  by  using  Sturm's  theorems  (See 
AppendixITT  for  a  particular  case). 

The  equation  for  yx  can  be  written  as: 

t  t 

Yi  =  y\  (0  ^  1  +  J  Va  (s)  f(s)  Yi  (s)  dsj  -  ya  (t)  J  yx  (s)  f(s)  yi,  (a)  ds 

£  £ 


From  the  above  discussion  and  observing  that  yx  and  ya  are  positive  functions, 
we  can  write  the  inequality 


t 

Yl  (t)  «  X  (t)  [  1  +  J  (s)  f(s)  yx  (s)  ds  1 

£ 


and 


_ Vj _ 

t 

1  +  1  Ys  (*)  f<s)  yx  (s)  ds 

£ 


s  y i  (t) 


Since  the  integral  is  positive.  Multiplying  both  sides  by  yaf  and 
integrating  between  £  and  t: 
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~  *  1  +  j*  Ya  (»>  *(»)  Yi  (s)  da  *  exp  (  J  yx  (s)  ya  (a)  f(s)  da) 
Yi  (t)  £  '  £  ' 


Thus  we  have  proved  the  following  lemma. 

Lemma.  If  in  the  interval  (a,b) 


y  5  y+y  J  n(»)  y(s)  ds 


(4.3.8) 


then  on  (a,  b), 

t 

y(t)  s  y  (t)  exp(  J  uyds) 


where  y,  ^  ,  and  y  are  positive. 

Substituting  for  y1  ,  ya  ,  and  f  we  get  the  final  reault: 

Yi  (t)  *  y i  (0  e*P  C  J  w"^4  I  (u>‘l/V|  da  ) 

i  1 


(4.3.9  a) 


The  equation  for  ya  (t)  can  be  written  as: 

t 


ya  (f)  =  Ya  (t)  +  y X  (t)  J  y a  (s)  f(s)  ya  (s)  da  -  y a  (t)  f  yx  (a)  f(a)  y8  (a)  da 

£  £ 


Considering  yx  ,  ya  and  fyj  as  positive,  we  can  then  write  the  inequality: 

ya  (*)  s  y*  (t)  ^  1  +  J  y i  («)  f(s)  ya  (»)  da  J 

t 


By  an  analysia  similar  to  the  caae  of  yx  we  arrive  at  the  reault: 


ya  (*)  *  7a  (0  ®xp(  J  |  (tu  )"  J  da  ) 


(4.3.9  b) 


Equation  (4.3.9)  (a,b)  with  £  =  a,b  respectively  are  precisely  those 
obtained  by  Olvcr  for  e  =1.  Sharper  bounds  are  derived  later.  The  various 
cases  can  be  discussed  aa  follows. 
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t  t 

y(t)  =  Ayj  +  By a  +  yx  (t)  J  y a  fy  ds  -  ya  (t)  J  yx  fy  ds 

€  € 

when  fy  >  0. 

(i)  f>  0,  y>  0 


(4. 


Y\  (0  s  Y\  (t)  exp  (  J  yA  7a  f  ds  ) 

i 


B=0;  A=1 


Ya  (t)  *  ya  (t)  exp  (  J  yx  ya  f  ds  ) 

t 


A  =  0;  B=1 


(ii)  f  <  0  .  y<  0 


Yi  *  Y\  t A  +  J  y«  fYi  d«3 
i 

|yi|  *  |yi|  [a  +  I  |y«/|  ‘Hyi/  ds3 


IrJ  1/al  M 


A  +  J  |y»|  |f|  |yi|  ds 


5  IyxI  |ya|  Ul 


Hence  |  yx  (t)|  £  |  Vi  (t)  |  exp(  J  |*y  i|  |ya|  |fl  ds) 
and  |y3  (t)|  ^  [ya  (t)|  «acp(  J  |yxy8J|f|ds) 


when  fy  <  0. 


(iii)  f  >  0,  y  <  0.  If  aj{t)  is  such  that 
conclusions  on  errors  of  case  (ii)  hold. 

(iv)  f<  0,  y>0;  f  =  -g  and  g  >  0 

t 

y  =  Ayx  +  By a  -yl  /  ys  g  y  ds  +  ya 

£ 


yx  <0*  yx  <0*  ya  <°.  the 


;  yx  g  y  d. 
* 
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yi  =  yi 


C 

L  1  +  J  y2  g  Yi  ds*  *  ya  J  Yi  g  Yi  ds 

t  t 


Thus  yx  <;  Yi  exp  (  J  yx  ya  g  ds  ) 


Similarly  ya  £  ya  exp  (  J  yx  ya  g  ds) 

£ 


Hence  we  can  write,  under  the  restrictions  given  above,  when  f>  0: 


l 

|yi  (0|*  fyi  (0|  exp  (  j  JyiYajf  ds) 


j  ya  (t)|  *  I  Ya  (t) )  exp  {  J1  jyi'ya  i  |  ds) 


(4.3.10) 


when  f<  0 


\Y\  (Oj  5  |  Yi  (0  |  exP  (  J  jytYa  f  |  ds) 

t 

|y»  (t)|  5  |  Ya  (t)j  exp  (  j  |  yx  ya  f  j  ds) 


Uniform  bounds  can  be  obtained  as: 


o 

j  Ytra(0|  *  I'yi.a  (t)  |  exp  (  J  |  y  x  y  a  i  j  ds  ) 


(4.3.  11) 


t  b  b 

Denoting  Fx  =  J  |  Yj  Ya  f  J  ds  ;  Fa  =  J*  j  Y\  Ya  f  j  ds;  F  =  J|  yx  ya  f  (  dt 


and  using  the  inequality  (n+1)!  >  2  ,  we  can  obtain  bounds  which  are  less 
sharp:  if  F.  <  2  f 
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i  ***  /  2  F  •  _ 

[  Yi  I  s|  Yil  - L  i  =  1,2 

2  -F 

i 


or  |yi,a|  *  J  y I.a  I 


2F 

2-F 


if  F  <  2 


(4.3.  12) 


Case  (2)  Oscillatory  Case,  uo  >  0  and  characteristic  roots 
are  imaginary.  As  before: 


y(t)  =  Ayj  +  By 3  +  J  L  Y  i  (t)y  a  (s)  -  %(t)  ^  (s)  3  f(s^  y(s)  ds 

Let  us  write  y1  and  ^a  as  'yi  =  _1/4  8  in  fl'(t)  and 

~  /  t  i 

ya  =  w  cob  a(t)  where  a(  t)  =e  J  u)*dt.  And  therefore: 


y  (t)  =  Yi (t )  +  J  C  sin  aft)  cosa(s)  -  coa  Q(t)  sin  ff(s)]  u>"^4(t)  u>  "^(s)  f(s)  y(s)  ds 
1  « 

t 

=  Yi  (t)  +  J  sin(  o(t)  -  o(s)  )  uj  ~^4  (t)  u;  _^4  (s)  f(s)yl  (s)  ds 

£ 


Yl  s  7i  +  j  |  u«  _1/4(t)  |  |uo  *  (s)J  J  f(s)j  jy^sjJ 


ds 


Let  oi  be  bounded  in  (a,b)by  M.  Hence  J  y t  j  £  M 


{ yi  J  ^  j  y  i 


+  m3  /  |fl|yJ  ds 

e 


By  Gronwall's  lemma 


k  k 

[Yi  (t^lYi  (t)|+  Ma  |  j  yx  (s )  f  |  exp(  J  Ma  J  f  |  dT  )  ds 


Y\  (^Myi  (0|  -  Ma  |  ~  L  exp(  |  M3  f  J  dT  )-i  ds 
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i.  e. : 


Yi  (t)l  -  M*r  e*p(  j1  Ma  |  f  |  dr  )] 


Iyi  (t)|+  M*C  exp(  J  M3  |fj  dr)  -1  ] 

£ 

Similarly: 

|ya  (OHya  (t)|+  M“r  exp(  J  Ma  J  f  j  dr  )  -P 


(4.3.13) 


Sharpening  of  Bounds;  Oscillatory  Case.  These  bounds  can  be 
sharpened  as  follows.  Applying  the  Liouville  transformation 


y  =  tt’  ~^*z, 


£  =  J'tw  dt 


,  transforms  (4.  3.1)  into 


z»  +  (ca  +  fa  )z  =  0 


(4.3.14) 


where: 


f,  •  «-**  <»  I 


(4.3.15) 


Solution  of  (4.  3.14)  can  be  written  as 


:(£  )  =  A*j  (£  )  +  Bza  (£  )  +  j  (zx  (£)  *3  (•)  -  *a  (O  *i  (*)^)fi  (■)*(»)  ds 


where  ,  za  are  the  solutions  of: 


z"  +  e  a  *  =  0 


With  the  same  restrictions  on  initial  conditions  as  before,  and  choosing 


~  *  M.  ~  . 

Zj  =  —  cos  c£,  za  =  sin 


so  that  the  Wronskian  W(zj,”za  )  =  1,  we  can  write 
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~  1  *■* 

2a  (S  )  =  2  a  (£  )  +  —  J  (cos  sines  -  sin  c|  cos  e  $)  fj  (s)  z.  (a)  ds 

O'  ' 

!  f 

=  2aU)+-  J  (-sin  (?  s  -  f  £)  )  fx  (s)  za  (s)  ds 

a  * 

i.  e. : 

i  ^ 

z3  (i  )  *  za  ($  )  +  -  j  jfj  (s)||*a  (s)jds 
Further: 

i-  «»/  5  |  ({  )J  +  |  |f1(.)||«a(.)|  d. 

* 

By  Gronwall's  lemma 

<£)H  |s»«>| +  7  [  l*° (■>/  k <*>/  «*  j 

i.  e. : 

|2a  U)|  *  |za  («)  |  +  J  “  ~  exP(  ~  J  |fx(r)|dr)da 

since  za  (s)  is  bounded  and  this  property  is  used  only  in  the  integral. 

Thus: 

(za  (€  )  |  ^  |  za  (€  )  |  -  [  exp(  j  J  |  fj  (T  )  |  dr  )  ]  ^  _  ^ 

i.e.:  |ca  (£)j  £  |sa  (£)|  +  exp  (  1  J  {  fj  (T  )  I  dT  )  -1 

Q  1 

Now  transforming  back  to  the  original  variables,  and  noting  that: 
d£  =  ®^dt 

JyalOj*  u»  {  sin  (e  J®^dt)  +  exp{  ~  J  |  u>~^4(ui )"  j  dt)  -1  J 
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This  can  be  written  in  the  form 


y 2  (t) 


(4.3.  16) 


where 


a  -l/4  |  dtJ  )  .1 


and  c  is  in  (a,  b).  This  is  Olver's  result  for  the  envelope  of  the  solutions. 

The  other  solution  can  be  studied  in  the  same  manner  and  similar  bounds 
can  be  obtained.  This  also  embraces  a  theorem  of  Wintner  (Ref.  45) 
which  is  derived  for  the  case  in  which  e  =  1  and  the  interval  is  infinite. 
Asymptoticity  of  the  solution  for  e  -  “  is  thus  demonstrated. 

Sharpening  of  the  Bounds:  Nonoscillatory  Case.  We  will  now  consider 
the  equation  : 

y"  -  e  2  a  y  =  0  (4.3.  17) 

As  before  the  Liouville  transformation  leads  to  the  equation 


z"  -  (ca  -  f)z  =  0  (4.3.  18) 

with  f;  being  given  by  (4.  3.  15).  Again  treating  (4.  3.  18)  as  an  inhomogeneous 
equation,  the  solution  is  written  as 

£  ^  * 

z  =  Azfx  +Bz”a  +  zx  (£  )  |  ^  fj  z  ds  -Zs(C)  J  Zj  (s)  fj  z  ds  (4.3.19) 

where  zx  and  z3  are  the  solutions  of  the  equation 
z "  -  e  3  z  =  0 


chosen  such  that  the  Wronskian: 

W(  zx  ,  za  )  =  1 
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Thus  z 


(4.  3.20) 


-e  £ 


e  £ 

e 

IT” 


From  (4.  3.  19)  we  can  write: 


Z.  =  Zj  +  Zj 


J1  Zg  Zj  ds 


I 


Zj  Zj  ds 


(4 


For  the  conditions  under  which 


ft  zt  s  0 


we  can  write  the  inequality: 


i 

£  Zj  +  Zj  I*  z3  f  Zj  d s 
a 


Now  lemma  (4.  3.  8)  leads  to  the  result 


(4 


zi  <;  z:  exp(  J  zx  za  f  ds) 


for  f  >  0.  On  substituting  from  (4.  3.2  0)  and  transforming  back  to  the 
original  variables,  the  following  result  is  obtair  d. 


y:  (t)  =  uu  -V4  exp  (  -  R  Jx^dt)  f  1  +  Ej  } 

j  t 

where  |  E1  |  £  exp  (  —  j  a’-1/4  (  a-  _1^4  )"  ds )  -1 


(a) 


(b) 


(4, 


Similarly  the  following  result  is  obtained  for  the  other  solution. 


where 


ya  (t)  =  iu  exp(  e  Jw^dt)  f  1  +  Ea  } 

1  b 

(Ea|s  exP(  2  /  u; 'V'4  (u,'^4  )"  ds)  -1 


(a) 


(b) 


(4. 
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.3.20) 


.3.21) 


3.21) 


3.22) 


\ 


These  bounds  are  precisely  those  obtained  by  Olver  (Ref.  28); 
here  they  have  been  obtained  by  a  direct  procedure,  but  under  the  conditions 
f  >  0  and  z  >  0.  For  any  given  a  (t)  the  sign  of  f(t)  can  be  checked.  The 
sign  of  the  solution  must  be  determined  using  Sturmian  theory. 

(ii)  Second  Order  Equation:  Noncanonical  Form.  The  equation 


y"  +  c  x i  y'  +  e  u^y  =  0 


(3.2.  19) 


can  be  treated  in  a  similar  manner.  The  time  scales  appi  ^xipiations 
y  i  f  a  =  D  '^4  exp  (  c  J  (  -  ^  ±  )  dt 


are  the  exact  solutions  of  the  equations 


y  41  +  e  v  i  y 1  +  (c  3  uj  q  +  f }  y  =0 


•where 


f(c  ,t)  =  (D_l/4)" 


and  B  is  the  discriminant  =  a  j3  -  4u)^.  The  original  equation  is  therefore 
written  in  the  form: 


y"  +  e  uj  i  y1  +  (e  3  u:  q  +  f)y  =  fy 


The  general  solution  of  this  is  given  by: 


y(t)  =  Ay;  +  Bys  +  yx 


f  ya  fy  ds 

J  W(ys  ,y*  ) 


-  y a 


/ 


Yi*) — 

W(y8  .ft 


) 


The  equation  (3.2.  19)  is  not  in  self-adjoint  form  and  hence  the  Wronskian 
W(ya  ,  ft  )  is  not  a  constant.  By  Abel's  formula  (Ref.  60) 
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W(ys  .  Yi  )  =  exp(-e  J  «:  x  dt); 


and  therefore 


y(t)  =  Ayx  +  Bys  +  yl  J  ya  g  y  ds  -  ya  J  y i  g  y  d« 


where  g(c  ,  t ) 


f(c  ,t) 

W(ya  7^  ) 


(4.  3, 


By  a  similar  line  of  reasoning  as  in  the  self-adjoint  case  we  arrive  at  the 
bounds  (4.  3 .  1  0/  provided  f  is  replaced  by  g(  e,t).  It  can  be  verified 
that  this  leads  to  the  results  obtained  by  application  ofOlver'a  lemma  and 
stated  in  Theorems  1  and  2  (4.2.  11;  a,b). 

4.4  Third  Order  Equation: 

We  now  apply  the  ideas  used  in  the  last  section  to  higher  order 
equations  and  derive  error  bounds  for  the  approximation  obtained  by  the 
use  of  multiple  time  scales.  The  analysis  is  more  difficult  since 
oscillation  criteria  such  as  the  elegant  theorems  of  Sturm  are  not  readily 
available  for  higher  order  equations.  As  Hartman  says,  "The  difficulty 
arises  from  the  fact  that  the  theorems  of  Sturm  do  not  have  complete 
analogues"  in  higher  order  systems  (Ref.  23,  p.  384).  Nevertheless  the 
asymptotic  solutions  are  oscillatory  if  the  characteristic  ’-oots  are 
complex  and  monotonic  (in  the  sense  of  having  at  most  one  zero)  if  the 
roots  are  real  (Ref.  46).  Conditions  on  the  coefficients  can  be  determined 
such  that  the  solutions  are  disconjugate  (Ref.  23,  p.  384). 

We  shall  start  with  the  third  order  equation  (3. 3.  9  a),  i.  e.  : 


y"'+c:au)1y,  +  f3uby  =  0  (3.3,9 

o 

in  the  interval  (a,b). 

The  approximation  via  the  time  scales  method  is  obtained  from  (3.4.21  c)  as 

t 

y(t)  =  D(t)  exp(e  J  x  dt)  (a) 

S 


where  £  =  a 
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where  D(t)  -  (ix"  + 


(b) 


and  x  is  a  characteristic  root.  Now  y  satisfies  exactly  the  equation 


(4.4.  1) 


whe  re 


y'"  +  e“  ux  y1  +  (r  3  x  ^  +  f(e  ,  t)  )  y  =  0 


2  1,  ,  „  3x'D'  3xD"  D'" 

f(e't)=  t  (X  +—  +—  >  -  5“ 


(c) 


(d) 


Equation  (3.  3.  9  a)  is  written  as 


y"1  +  c  2  u-  i  y1  +  ( c  3  u.  ^  t  f)y  =  fy  (4.4.2) 

with  f  as  given  by  (4.4.  1  b).  Again  when  f(c  ,t)  is  small  in  comparison 
with  c  3  x  ,  solutions  of  (3 .  3 .  9  a )  and  (4.4.1  c )  will  nearly  be  the  same. 

The  integral  equation  corresponding  to  (4.4.2)  is  jiven  as 

t 

y(t)  =  Ayx  +  By3  +  Cy3  +  y*  J  (y'a  y3  -  y'3  y3  )  fy  ds 

S 

t  t 

-  ya  f  (y’i  y3  -y3yt  )  fy  ds  -  y3  J  (y'a  ft  -y'i  ya  >  fy  ds 

(4.4.3) 

with  -W(y1  ,  ya  ,  y3  )  =  1 .  The  discussion  follows  the  one  for  second  order 
equations  and  y(t)  is  chosen  to  be  that  solution  of  (4.4.1  c)  which  is 
characterized  by  the  same  conditions  as  y(t)  at  any  point  £  in  (a,b). 

As  before  let  the  conditions  A-l,  B  =  C  =  0;  B=l.  A  =  C  =  o  and  C  =  l.  A  =  B-0 
correspond  respectively  to  the  subscripts  1,  2,  3.  Further,  we  notice 
that  for  the  canonical  equation  the  sum  of  the  characteristic  roots  will 
be  zero. 

N'onosr illatory  Solutions.  In  this  case  the  characteristic  roots  are  real, 
Case  (1)  fy  >  0.  Let  the  roots  be  ordered  such  that: 

x3  >  xx  >  x3  (4.4.4) 


96 


Hence: 


Yi  (0  =  y:  -  1  +  J  (y 2  y3  -  y'3  ya  )  fyi  ds!* 

i 

t  t 

-  y2  f  (y'i  y3  -y'3  yx  )  fYi  ds  -  y3  J  (y'a  yi  -  y'i  ya  )  fyx  ds 
«  € 

Let  us  consider  now  that  the  conditions  are  such  that  yj  >  0  in  the 
interval  (a,b).  Also  in  y  which  is  the  approximation  for  y  the  explicit 
amplitude  variation  expressed  by  D  is  understood  to  be  positive. 

Since  y  (t )  =  D(t)  exp(  e  r  x  dt)  where  D(t)  =  (3X2  +  uu  x  ^  , 

€ 

y'<t)  =  (  £  +  ex> v 


y'a  y3  -  y'3  y3  =  f  ( 


)  +  e  (Xg  -  x3  )}  y3  y3 


Now  ( 


D'3  .  d  ,  ^  d 

of 1 1  s- tn  Da  -  d r  ^ 


L  ~  /  n  (3x%  +uj1)-~-in  (3xaa  +  aJ:  )  j 


=  4: 


-  (3xa  2  la,)1  (3x53  +  ^  )’  3 


(3x  2  +a'j)  (3x  3  +  u,j) 


Now  if  xa  >  xa  >  x3  and  the  roots  do  not  coalesce  the  root  variations 

D 1  D 1 

may  be  assumed  to  be  of  the  same  order  of  magnitude.  Thur  ( — - - 2-)  2  0; 

_  d3 

and  in  any  case,  for  large  e,  %'y3  -  y3'  ya  >  0. 

In  the  special  case  when 


'i  3  (t) 


«s  2  (t) 
o 


=  constant, 


D(t)  has  a  particularly  simple  form  given  by  D(t)  =  ^  (t)  =  (t)  (a) 

and  y's  y3  -  y3  ya  =  e  (xjj  -  x3  )ya  y3  =  c  (xg  -x3  )  ujf1  (t)  exp(  c  f  (Xj,  +x3  )  ds)  (b) 


97 


Simila  rly 


Y  l  Ya  -  y  3  Yi  =  e  (xj  -  x3  )yx  y3 


(c)  (4.4.5) 


y  a  Yi  -  y  1  Ya  =  e  (x3  -  x:  )y,  y3  (d) 

A/ 

The  exponentials  in  y  have  real  argument  and  the  following  inequalities 
can  therefore  be  written,  for  large  e  . 

rJ  /J 

y  a  Ya  -  y  3  Ya  >  0 

y’x  Ya  -  Y'aYi  >  0  (4.4.6) 

y'a  Yi  -  y'i  Ya  >  0 

t 

Hence:  yx  (t)  £  7  (t)  C  1  +  J  (y's  y3  -  y'a  ya  )  fyt  ds] 

€ 

By  applying  the  lemma  (4.3.8)  this  immediately  leads  to: 

t 

Yi  (t)  S  9i  (t)  exp(  J  Yx  (y'a  y3  -  y'3  7a  )  C  ds)  (4.  4.  7  a) 

€ 

uu  3 

But  when  — =  k  =  constant 

U) 

■> 

Yl  (Y  a  Ya  -  Y  3  Ya  )  =  e  (xa  -  x3  )  y:  y3  y3 

=  <•  (xa  -x3  )  D3  (t) 

since  Xj  +  Xg  +  x3  =  0  .  Therefore 

Yl  (Y'a  Ya  “  y'a  Ya  )  =  eoi1'3/3(t)  (x,  -  x3  ) 
since  :  D  =  x  j  ^  ( t ) 
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The  approximation  can  be  written  as 


y i  =  y i  (i  +  Ei  ) 


v  ..4,  7  b) 


whe  re: 


I  Ex  I  s  exp  [  e  J  ui  -^a  (s)(xa  -  x3  )  I  f(e  ,  s)  j  ds  }  -  1 


Now  for  A=C  =£0,  B  =  1 ,  from  (4.4.3)  we  can  w rite  the  inequality : 


(‘)s  ya  (t)  +  yi  (t )  J  (  vi  y3  -  y a  )  fya  d8 


Applying  Gronwall's  lemma  to  the  ratio  ya  (t)/y^  (t),  we  obtain: 

(4.4.8  a) 

t  t 

y3  (t)  S  ya  (t)  +  yi  J  ds  ya  (9a  y3  -  y3'  ya  )  i  expf  J  yx  (£'  y3  -  ya  )  fdT  J 

f  s 

When  y1  ,  ya  ,  y3  are  bounded  by  Lj  ,  Mj  ,  N,  respectively,  a  simplified 
bound  can  be  obtained  as 


ya  (0  £  9a  (0  +  yi  C  exp(  J  y:  (ya'  y3  -  y3'  ya  )  f  ds)  -  lj  (4.4.  8  b) 

which  when  as3  ~  x  3  ,  can  be  written  in  the  form 

’  i  0  > 


ya  (t)  =  ya  (t)  +  Ea 


where: 


EaN  >i  (t)|exp(  r  j  x'i  '^a(xa  -  x3)  j  f(e  ,  s)  J  ds)  -  1 1  (4.4.8  c) 


Similarly,  choice  of  A=B  =  0,  C  =  1  leads  to  the  bound  on  y3  (t): 


v 

y3  (t)  *  y3  (t)  +  yx  |  (%,'  y3  -  y3  ys  )  fy3  ds 

t  t 

y3  (0  s  y3  (0  +  Yi  (t)  J  (>&'y3  -y3'y8  )y3  *  exp(  j  yt  (£'  y3  -yj'ya  )fdT  )  ds  (a) 

i  » 
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A  simplified  bound  is  obtai  .ed  as 


y3  (t)  *  y3  (t)  +  yx  (t\jexp(  j  y*  (%'y3  -y3'y a)fda)~lj  (b)  (4.4.9) 
or  alternatively,  when  ay3  =  k  a)Q2  ,  k  being  a  constant 

y3  (t)  =  y3  (t)  +  e3 

where: 

|e3|  <  yx  (t)  [  exp(  C  ^  o:1'^2(xs  -  x3  )  |  f  ( e  ,  s )  |  ds )  -  1  ]  (c) 

Case  (2 )  f  y  <  0.  With  reference  to  (4.4.3)  let  the  integrals 
be  denoted  by  IT  ,  Ia  ,  I3  in  that  order.  When  fy  is  negative  the  signs  of  all 
the  integral  terms  are  changed.  Therefore,  if  these  terms  are  to  have  the 
same  signs  as  in  the  previous  case  the  ordering  of  the  roots  ha^  to  be 
changed.  In  place  of  (4.4.5)  the  following  inequalities  are  needed. 

vb  ya  -  y3  ya  <  °;  %  y3  -  y3  yi  <  o  (4.4. 10) 

The  roots  are  therefore  ordered  such  that  x3  >  >  x3  ;  i.e.  the  roots 

xa  ,  x3  are  interchanged  w.r.t.  case  (1). 

The  analysis  is  carried  out  in  the  same  manner  as  for  case  (1)  and 
the  results  (4.  4.7),  (4.  4.  8),  (4.  4.  9)  and  (4. 4'.  1  0)  the  following  approximation 
theorem  has  been  obtained. 

Theo-om  5.  In  a  given  interval  (a,b)  the  differential  equation 
y'"  +  r  2  tr  i  (c  ,  t)y'  +  e  3  w  ( e  ,  t)y  =  0  (4.4.11) 

possesses  the  solutions  yx  ,  y3  ,  y3  such  that 

yi  =  Yi  (1  +  El  );  y3  =  ya  +  Ea  ;  y3  =  y3  +  E3  (4. 4.  12) 


100 


where 


(4.4.  13) 


y.(t)  =  D.  e:  (r  X.)  ;  i=  1 , 2,  3 
1  1  a  1 

D.(t)  =  (3x3  +  a)!  )"?  (4.4.  14) 

and  x.(e  ,t)  are  the  roots  cf  the  characteristic  equation 


x°  +UL',(e,t)x+u'  ( e  ,  t )  =  0 
o 

l  ^  |  t 

|  Ex  |  £  e  (Fx  )  -  1 

|e2|*  Vi  J  ds  yj’1  y2  Fx  (e  ,  s)  eg  (Fx  ) 
a 
t 

1  E3 )  ^  Yi  j  ds  y:_1  y3  Fi  (e  ,  s)  eg  (Fx  ) 
a 

Fx  (e  ,t)  =  yx  (t)  (%,'  y3  -  y2  )  f(c  ,t) 


,  e2  a,  ,  3x'D 1  3xD"  ,  D’ 
-c()t  +  5-  +  D"  »  -D 


,3  P  r 

and  e  (y  )  is  an  operator  defined  by  e  (y  )=  exp(  I  y  (X  )  dX  ) 
Ot  Ot  J 


(4.4.  15) 

(4.4.  16  a) 

(4.4.  17) 

(4.4.  18) 

(4.4.  19) 


provided  that  the  following  conditions  are  met: 

(i)  the  roots  x.  are  real  and  distinct  (4.4.20) 

(ii)  fy  is  of  constant  sign  in  (a,b)  (4.4.21) 

and  further  the  roots  are  ordered  such  that  x3  >  xt  >  x3  if 

fy  >  0  and  x3  >  x1  >  x2  when  fy  <  0. 

The  derivatives  in  f(e  ,  t)  are  assumed  to  exist  and  the  interval  may  • 
be  infinite  if  the  integrals  in  (4.4.  16)  converge. 

The  condition  (4.4.20)  is  satisfied  if  the  discriminant 
q3  +  r3  <0  where  q  =  ,  r  =  - (4.4.22) 

The  validity  of  (4.4,21)  must  be  ensured  by  examination  of  initial 
conditions  and  oscillation  criteria. 
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Uniform  bounds  can  be  obtained  by  replacing  the  upper  limit  of 
integration  by  b  in  (4.4.  16  a). 

Simplified  bounds  can  be  obtained  by  noting  that  (n+l)l  2  2n.  Thus: 
®  n+i 

«p(„)-i=E  “—y 


a>  n+i 

<  s  V 
0 


~  .  a  n  a 

=  *  r  {  )  = 

0  1_2 


Hence  J  Ej  |  £ 

1  ~2 


if  a  <  2 


(4.4.23) 


where: 


a  -  J  Fj  dX 


Similar  bounds  can  be  determined  for  the  other  two  relations  in  (4.5.6). 

If  y j  ,ya  ,y3  are  bounded  by  Lj  ,  Mj  ,  Nx  respectively,  the  following 
relations  can  be  used  in  Theorem  3. 


\Ea  1 5  ^  YX  [  e‘  (F,  )  -  1  ] 


lEal^^  Yi  <t)0*(Fx  )  -  1  ] 


(4.4.16  b) 


o 

(|5  . 

Further,  in  the  case  when  — *5-  =  constant  the  bounds  turn  out  to 


be  simple. 


D(t)  =  (e  ,  t) 


(4.4.16  c) 


Fx  =  cu>  ~q'a(xa  -  *3  )  f(c  .t) 


(4.4.17  c) 
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4.5  Oscillatory  Case 


In  this  case  the  characteristic  roots  consist  of  one  real  root  and 
a  pair  of  complex  conjugate  roots,  symmetrically  placed  about  the 
origin  of  the  complex  root  plane.  Let  x^  be  the  real  root.  Equation  (4.4.3) 
can  now  be  written  as 


y(t)  =  Ayx  +  Bya  +  Cy3  +  f  h(t,  s)  fy  ds 


(4.5.1) 


where 


h(t.  s)  =  ya  (t)  (y,  &  -  ya  ft')  (s)  +  y3  (t)(fa  ft*  -  yx  &  )(s).+?x  (t)(£ft'  £  )(s) 


Now  if  yx  (t)  and  ya  (t),  y3  (t)  (ya  and  y3  are  complex  conjugates)  are 
bounded  from  above  by  Lq  and  Mg  respectively,  then 

t  t  t 

Jh(t,s)fyds<  J  |h(t,  a)  fyjds  £  f  |px|  )  f  |  |y  |  ds 

€  €  K 


(4.5.2) 


where: 


Px  =  Mg  { |  yx  -  y3  ft'  |  +  |  ya  -  y*  ^  |  ya  *»'  -  y3  ya  |  J  (*•  5. 3) 


Now  choosing  A=l,  B=C=0;  A=C=0,  B=1  and  A=B=0,  C= 1  in  that  order,  the 
following  result  is  obtained  by  applying  Gronwall's  lemma: 


i  t 

lyiU  |y.(l)|  +  j*  |  y^s)  Pa  f  |  exp(  |  |  Px  |  |  f  {  dT )  di 

'  I  f  fl 


(4.5.4)* 


i  =  1,2,3,  where  Px  is  given  by  (4.  5.  3). 

Some  simplification  can  be  achieved  as  follows.  Let  an  operator 
e^  (x)  be  defined  such  that: 


(x)  =  exp(  c  J  x  dX  ) 

4  i 


Hence 


and: 


e^(-x)  =  e*(x) 


e£  (x)  e3  =  es<x> 


(a) 


(b)  (4.5.5) 


(c) 
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The  approximation 


y(t)  =  D(t)  exp(  e  j  x  dX  )  =  D(t)  e1  (x) 


(4.5.6) 


y'(t)  =  (  —  +  e  X  )y 


Let  us  write  h(t,  s)  in  the  form: 


h(t,  s)  =  ya  (t)  (  yj  y3'  (s)  -  y3  y/  (s)  )  +  y3  (t)(y2  (s)  -  yx  y^(s)  ) 


+  yx  (t)  (  Ya  Vs  (s)  -  y3  Ya  (s)  ) 


(4.5.7) 


We  use  the  above  and  recall  that  D(t)  is  understood  to  be  the  magnitude  of 
(3xa  r  uj :  )  and  under  certain  conditions  (derived  in  Appendix  II),  is 
invariant  w.  r.t.  the  roots.  In  particular  if 


=  constant 


ui  3  (t) 
o 


4 

D  (t  )  =  U!X  ?  (t) 


y.(t)  =  D(t )  e^  (x. )  =  u);  ?  (t)  e^  (x.) 


(4.5.8) 


Using  (4.4.5),  (4.  5»  7)  can  now  be  written  as : 


h(t,  s)  =  e  D(t)  el  (x^  )  ^x3  (s)  -  x:  (s)^)  D3  (s)  e*  (x;  +  x3  ) 

^  »  g 

+  c  D(t)  e  (x3  )  (  Xj  (s)  -  xa  (s))  D3  (s)  ef  (x1  +  Xa  ) 
+  e  D(t)  e^  (x1  )  (x3  (s)  -  xa  (s)^  D3  (e)  e*  (x^  +  ) 
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Noting  that  x:  (t)  +  xa  (t)  +  x3  (t)  =  0  and  (4.  5.  5;  b,  c)  leads  to  the  simplification: 


h(t,s)  =  r  D(t )  w  Da  (s)  e*  (x:  )  {x3  (s)  -  xa  (s)^)  +  (xa  )  ( x3  (s)  -  x:  (e)^ 


+  (x3  )  (s)  -  xa  (s)^ 


(4.5.9) 


^  n  t 

Now  y(t)  =  y(t)+J  h(t,  s )  fy  ds.  If  e^(x.),  i  =  1,2,3  and  D(t)  are  bounded 
respectively  by,  say,  Lg,  M3  ,  and  R,  and  defining 

Pa  (s)  =  «■  R  f  a^1  (s)  L3  (  -x1  )(x3  (s)  -  Xg  (s)  )  | 

+  |  M3  u-V  (s)  (-xa  )(x3  -  x1  )  +  j  M3  u'*j  (s)  e®  ( -x3  )(xt  -  xa  )  j  }  (a) 


we  get  the  relation: 


(4.5.  iO) 


|y(0|^  | y <t) (+  J  jPa(  |f|  (y I  ds 


Cronwall's  lemma  leads  to  (4.5.4)  with  Px  replaced  by  Pa  . 

The  following  theorem  can  therefore  be  stated. 

Theorem  6.  With  the  following  modifications  theorem  3  will  hold. 


|y.|  =  iy.(t)|  +  E.  i=i.2, 3 


whe  re 


y.(s)  F j  f  e  (  P:  f)  ds 
l  s  4 


(4.5.  11) 


(4. 5. 12) 


Pi  =  Ma  {  |  yx  \b'  -  y3  yi’  |  +  |  y2  Vi'  -  Yi  y?  |  +  -^r|  y3  y3  -  y3  &)}  (4.  5. 13) 


provided  that  roots  are  distinct,  a  pair  being  complex  conjugates.  With 
the  above  modifications,  the  relations  (4.4.  11)  to  (4.4.19)  hold. 

K  is  any  point  in  the  interior  of  (a.b).  Lg  Ma  are  the  upper  bounds  of 
y.  and  ya  ,  y3  respectively. 
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(a)  The  characteristic  roots  are  complex  and  distinct  if  the  discriminant 

q3  +  ra  >  0  ;  q  =-ji  .  r  =  -  {4.4.22  b) 

If  qa  +  r3  =0  then  the  roots  are  all  real  and  at  least  two  of  them  are  equal. 

(b)  Uniform  bounds  are  obtained  by  having  a,b  as  the  limits  of 
integration  in  (4.5.  12). 

(c)  Simplified  bounds  can  again  be  determined  similar  to  the  case  of 
theorem  3;  i.e.  similar  to  (4.4.23). 

When 

3 

U)  1 

— =  constant 

u>  3 
o 

in  (4.  5.  12)  can  be  replaced  by  ,  which  is  given  by  (4.5. 10  a). 

4.6  Summary  of  Chapter  and  Conclusions 

Error  bounds  for  the  second  order  equation  are  derived,  first  fot 
the  noncanonical  form  and  this  is  done  using  Olver's  fundamental 
approximation  theorems.  This  reveals  a  difference  in  errors  in  comparison 
with  the  standard  LG  error  bounds.  Next,  Olver's  approximation  theorems  for 
the  self-adjoint  equation  are  rederived,  though  more  restrictively,  in  a 
direct  way  that  is  believed  to  be  new.  Further,  bounds  for  the  noncanonical 
equation  are  also  derived  using  this  method  and  are  shown  to  be  equivalent 
to  the  earlier  result  using  Olver's  theorems.  Simplified  bounds  are  obtained 
and  the  conditions  for  these  are  stated. 

The  third  order  1.  d.  e.  is  studied  next  and  two  approximation  theorems 
are  proved,  using  the  familiar  notions  of  variable  characteristic  roots  .  For 
one  class  of  equations  for  which 

3 

Oli 

— K-  =  constant 

U) 

o 

the  error  bounds  have  a  simple  form.  Simplified  bounds  of  the  type  obtained 
by  Blumenthal  for  second  order  equations  are  also  derived  under  appropriate 
conditions. 
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Though  the  approximations  obtained  in  Chapter  III  require  that 
u).(t)  be  independent  of  e  ,  in  the  error  analysis  this  limitation  may  be 
overlooked  and  the  more  general  u>.(e  ,t)  can  be  considered.  The  assumed 
condition  that  f(c  ,t)  is  a  continuous  function  of  t  is  not  essential.  The 
same  proof  applies  if,  for  example,  f(e  ,  t)  has  a  finite  number  of  discontinuities 
in  (a,b).  But  then  higher  derivatives  of  y(t)  may  be  discontinuous  at 
the  points  of  discontinuities  of  f(e  ,  t). 

The  conditions  under  which  fy  is  of  constant  sign  for  theorem  3  must, 
however,  be  determined  from  other  considerations  such  as  stability  theory 
and  oscillation  criteria  (Ref.  61,  46,  23).  The  bounds  for  the  oscillatory 
case  are  not  as  simple  in  form  as  for  the  second  order  equation. 

Besides,  it  may  not  always  be  possible  to  evaluate  the  bounds  in  a  closed 
form;  but  the  theorems  are  still  useful,  as  they  essentially  reduce  the 
problem  of  estimating  the  error  in  an  approximate  solution  of  a  differential 
equation  to  the  much  easier  problem  of  evaluating  a  definite  integral. 

It  is  felt  that  a  similar  approach  may  lead  to  useful  results  in  the 
case  of  thi  ci  order  noncanonical  form,  as  well  as  for  higher  order  equations. 
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CHAPTER  V 


EXAMPLES  AND  APPLICATION 

The  scheme  of  approximation  developed  in  Chapter  III  will  now  be 

applied  to  examples.  The  arrangement  is  as  follows.  The  first  part  of 

this  chapter  consists  of  examples  with  analytically  known  behaviors;  the 

examples  are  so  chosen  as  to  highlight  the  application  of  the  method  and  to 

afford  an  analytical  treatment.  Notice  that  even  though  the  coefficients  do 

not  always  completely  conform  to  the  ordering  assumed  in  our  approximation, 

considerable  information  is  obtained.  We  also  consider  a  special  equation  of 
th 

n  order.  The  equations  are  studied  in  the  noncanonical  form  directly 
without  the  necessity  of  transforming  them  into  the  canonical  form  which 
would  be  amenable  to  the  L  G  treatment. 

The  latter  half  of  this  chapter  is  devoted  to  an  actual  physical 
problem,  viz.  the  analysis  of  the  dynamics  of  VTOL  aircraft  through  the 
transition  from  hove  •  to  forward  flight.  The  problem  of  aircraft  dynamics 
in  unsteady  flight  has  been  treated  in  detail  by  Curtiss  (Ref.  27)  in  a 
recent  work.  Therefore,  at  present  we  will  not  go  into  the  many  aspects 
of  VTOL  dynamics,  but  will  emphasize  the  use  of  time  scales  to  provide  a 
uniform  description  of  the  transition  dynamics. 

5.  1  Examples  with  Known  Solutions 

1.  The  asymptotic  behavior  of  the  zeroth  order  Bessel  function  can  be 
recovered  directly  from  the  governing  equation: 


y"  +  -y'  +  y  -  0 


(5.  1.1) 


The  characteristic  roots  are: 


x=  -  it*1' 1  -tp* 


From  (3.4.21  c)  the  approximation  is  written  as 


V  (t)  =  (  1  -  ^p)  y*  exp{  'n  t  ±  i  J'  (1  -  dt  ) 
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1. e.  y  (t)  ~  t  5exp(  ±it)  (5.1. 

which  is  the  correct  asymptotic  behavior  as  t  -*30  (Ref.  50  ).  The  first 
and  third  coefficients  are  the  same,  while  the  second  coefficient  steadily 
decreases  as  t  . 

Consider  the  transformation  t  -  p  s,  Now  if  s  =  0(1  ),  as  t  —a> , 
e  Equation  (5.  1.  1 }  is  transformed  into 

sy"  +  y'  +  C2  sy  =  0 

with  -  >  >  1.  The  time  scales  formula  describes  the  correct  behavior  of 
the  solutions. 

2.  We  may  now  consider  the  confluent  hypergeometric  (or  Kummer's) 
equation,  which  has  a  number  of  engineering  applications  (Ref.  64).  The 
equation  is  given  by 

t  y"  +  (b-t)y'  -  ay  =  0 

where  a,  b  art  constants.  The  characteristic  equation  is  : 

xa  +  (-  -  1  )x  -  -=  0 
t  t 


For  t  -*®  ,  the  characteristic  roots  are  asymptotically  given  by: 


xi 


1 


+ 


a  -b 
t 


x 


a 

t 


From  (4.3.21)  .he  approximations  are  given  by: 

yi-  °i  (TQ)  exP  (  J  xt  dt  )  -  cl  e*  t^a 
Ys  =  Qa  <T0>  exP  (  j  xi  dt  )  ~  eg  t  a 
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These  describe  the  correct  asymptotic  behavior  of  the  solutions  as 
t  -*  »  ,  as  can  be  verified  from  Ref.  50,  63.  It  is  seen  from  the 
characteristic  equation  that  the  second  coefficient  approaches  a  constant, 
while  the  third  approaches  zero  as  t  -*®.  Notwithstanding  this,  the  correct 
behavior  is  obtained  from  the  approximation. 

Kummer's  equation  can  be  transformed  into  the  canonical  form  and 
is  known  as  Whittaker's  normal  form  (Ref.  63).  This  allows  the  application 
of  the  L  G  approximation. 


3.  We  will  now  consider  a  simple  third  order  equation  containing  a  large 
parameter  e  ,  viz.  : 


rlM 


y  =  0 


(5.1.3) 


This  is  of  a  type  which  is  fully  amenable  to  our  approximation.  The 
characteristic  roots  are: 


X  =  - 


l 

t 


_1 

t 


From  (3. 4. 21  c): 


y .  (t)  =  —  exp  (  e 
x 


h 


dt  ) 


i  •  c  •  * 


y  (t)  =  t(  cx  t  f  t 


,  14-1/3  4  ,  1  -i/3  , 

e  (  — 2 - )  f  (  — ) 


+  c3  t 


)  (5.1.4) 


The  exact  solution  of  (5.1.3)  can  be  obtained  in  the  form  y(t)  =  tm. 
It  can  be  easily  verified  that  (5.  1.4)  gives  the  correct  asymptotic  behavior 
for  large  e  . 

th 

The  spec:al  n  order  equation  of  this  type  can  be  studied  similarly 
and  it  is  discussed  at  the  end  of  this  section. 

4.  "Double  Airy"  equation 

This  is  a  third  order  equation  which  is  satisfied  by  products  of  Airy 
functions  and  is  given  by: 
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y"'  -  4  ty 1  -  2y  =  0 


(5.1.5) 


The  characteristic  equation  is: 

x3  -  4tx  -  2  =  0 

In  order  to  determine  the  asymptotic  behavior  of  the  roots  as  t  ~*m ,  the 
equation  is  written  as 

fix3  -  4tx  -  2fi  =  0 
and  is  studied  in  the  limit  fi  -*0. 

In  this  case  the  principle  of  maximal  balance  says  that  the  term 
x3  must  be  balanced  against  -4tx  for  large  t  (Fig.  12  ii);  i.  e.  : 

x3  -4tx  =  0  or  x  (x*  -4t)  =  0 

The  three  roots  are: 

xi  a  =  ±  2t* 

» 

x3  =  0 

This  can  also  be  seen  starting  with  the  equation: 

x3  -4tx  -2=0 

For  real  coefficients  the  sum  of  the  roots  =  0.  Hence  if  one  pair  of  roots 
goes  as  ±  Z/t  ,  the  third  root  must  remain  at  the  origin,  as  shown  in  the 
sketch  below. 


Hence  the  approximations  are 

y.  =  (  3x^a  -  4t)  ^  exp  ifxi  dt) 

4  _  1 

i.  e.  :  y1>3~  t  ®  exp  (±  —  t  3/a  )  ;  y3  ~  t  (5. 1.6) 

1 1 1 


* 


2  2 

Now  (5.1.5)  has  the  exact  solutions  Ai  (t),  Bi  (t)  and  Ai(t)  Bi(t) 
where  A i ( t )  and  Bi(t)  are  Airy  functions  (Ref.  50).  The  coefficients  of 
(5.  1.5)  do  not  fully  conform  to  the  conditions  of  the  approximation; 
the  second  coefficient  increases  as  t  — °°  ,  but  the  third  remains  finite. 
Nevertheless,  (5.  1.6)  describes  the  correct  asymptotic  behavior  of  the 
solutions  as  t-°°  . 

5.  Consider  the  equation 

t y 1 "  +  3y"  +  ty  =  0  (5.1.  7) 

The  characteristic  equation  is  given  by 

x3  +  -x*  +  1  =  0  (5.  1.8) 

t 


The  roots  of  this  equation  as  t-°°  can  be  studied  using  the  exact  formula  for 
the  cubic  (Ref.5"0,p.  17).  However,  let  us  write  (5.  1.2)  in  the  form 


x 


3 


+  R  -  x3  +  1  =  0 
t 


(5.1.9) 


and  study  the  roots  as  6-0.  Using  the  principle  of  maximal  balance  (Fig.  12  i) 
we  find  that  the  equation 


x  +1=0 

O 


(5.1.10) 


determines  the  leading  behavior  of  the  roots  as  t-®  .  Hence  we  can  treat 
(5.  1.5'  as  a  perturbation  problem  and  determine  the  correction  by  going  to 
the  next  order.  Thus: 


x  =  x  +  6  xi  +  . .  . 
o 


3x  x.  + 

O  t 


V  - 


=  0 


(5.1.11) 

(5.1.12) 


giving: 


Xi  =  ■  r 


(5.  1.13) 


This  correction  to  the  roots  can  be  verified  from  the  exact  formula  also. 
The  time  scales  approximation  is  therefore  written  as 
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y.(t)  =  c.  (  x.  ^  +  Xj_  )  2  exp  (  f  x.  dt ) 

1  1  l  t  J  l 


(5.  1.14) 


where  c.  are  adjustable  constants.  On  substituting  from  (5.  1.11)  the 
solutions  are  found  to  be 


r.  (t)  =  — -  exn  (  f  x  .  dt  ) 

i  t  ‘  J  oi 


(5.  1.15) 


where  x^.  are  roots  of  the  equation  (5.1.10);  i.  e.  the  solutions  are: 


y  (t)  =  ^  exp  (  -t)  +  Si.  exp(  )  t  )  +  D3  exp(  )  t)  (5.1.161 


The  transformation  u  =  ty  leads  to  the  constant  coefficient  equation 


u"'  +  u  =  0 


which  can  be  solved  exactly.  The  expression  (5.  1.1b)  given  by  the  time 
scales  theory  is  the  exact  solution  for  this  example. 

6.  Special  equation  of  order  n. 

The  equation  chosen,  agaiiy  conforms  to  the  conditions  of  the 
approximation.  Consider  the  equation 


m £  +  <  r  ’  v  =  0 


(5. 1. 17) 


in  the  limit  of  c-*00.  The  characteristic  equation  is: 


n  1  . 

F  =  x  +  —  =  0 

n  tn 

The  roots  are  given  by:  x  =  exp  (  —  i  )  (  ^  ) 


(5.  1.  18) 


(5.1. 19) 


The  time  scales  approximation  is: 
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Simplifying: 

Yi  (t) 


n-l 


2 

c.  t 
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(e 


tr  if  n 

e  ) 


(5.  1.20) 


The  exact  solution  can  be  determined  as  follows.  We  look  for  a  solution 
y  =  tm.  On  substitution,  m  is  found  to  satisfy  the  equation 

m(m-l )(m -2 )  .  .  .  (m-n+1  )  +  c  °  =  0  (5.1.21) 

which  can  be  written  in  the  roc.  locus  form: 

n 

1  +  G(m)  =  0  where  G(m)  =  — : - - - ; - — .  (5.  1. 22) 

m(m-l)  ...  (m-n+1) 

The  locus  of  the  roots  of  (5.1 .3,1)  on  the  complex  m  plane  is  plotted  in 
Fig.  13.  For  large  e  the  roots  asymptotically  approach  the  lines  which 
are  at  ang  les  1 80°/  n;  i.e.  m  e  77  ^  °  .  Thus  the  exact  solution  has  the 
asymptotic  behavior  for  e  "*  ®  as  predicted  by  the  approximation 
(5.  1.10) 

We  will  now  go  on  to  the  study  of  the  transition  of  VTOL  aircraft 
from  hover  to  forward  flight. 


1  14 


Transition  Dynamics  of  VTOli  Aircraft 


5.2.  Preliminary  Remarks 

The  problem  considered  in  this  section  is  the  longitudinal  dynamics 
of  a  VTOL  aircraft  through  its  transition  from  hover  to  forward  flight. 

The  method  of  multiple  time  scales  and  the  formulae  derived  in  Chapter  III 
are  employed  to  obtain  approximations  to  the  solutions  of  the  equations  of 
motion  of  a  typical  VTOL  aircraft.  The  point  of  view  adopted  is  to 
linearize  the  nonlinear  equations  of  motion,  and  treat  the  coefficients  as 
variable  during  the  transition.  This  is,  for  example,  the  approach  used 
by  Curtiss  (Ref.  27). 

The  assumptions  and  rationale  of  the  physical  problem  are  based 
on  Ref.  27  and  47.  The  main  contribution  of  the  present  effort  is  intended 
to  be  a  difference  in  approach  and  a  more  uniform  description  of  the 
phenomenon. 

An  independent  effort  in  the  application  of  multiple  time  scales 
aircraft  dynamics  was  reported  by  Ashley  in  a  very  recent  paper  (Ref.  52). 

This  deals  with  the  linearized  aircraft  equations  of  motion  with  constant 
coefficients.  Approximations  were  obtained  with  simple  time  scales 
using  linear  clocks,  the  objective  being  "a  heightened  rationality"  in  the 
study  of  the  subject. 

The  present  approach  differs  from  that  of  Ashley  in  the  following 
respects : 

(i)  in  regard  to  constant  coefficient  equations,  using  linear  clocks 

exact  solutions  are  obtained  instead  of  approximations  (Ref.  Chapter  II) 

(ii)  since  linear  clocks  are  inadequate  for  a  large  class  of  problems, 
nonlinear  clocks  are  introduced,  particularly  in  the  case  of  l.d.e. 
with  variaole  coefficients.  The  nature  of  the  clock  functions 

is  determined  from  the  equations  themselves  and  depends  on  the 
domain  of  interest. 

It  is  felt  that  this  dissertation  demonstrates  the  usefulness  and 
flexibility  of  a  general  nonlinear  c.ock  function. 
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5.3.  The  Sy stem 


The  development  of  the  aircraft  equations  of  motion  as  related  to 
VTOL  dynamics  will  be  briefly  traced,  for  the  sake  of  completeness. 

The  motion  of  the  aircraft  is  considered  with  reference  to  a  system  of 
body  axes  fixed  in  the  vehicle.  Fig.  14  describes  the  axis  system  for  a 
tilt-wing  VTOL  vehicle. 

The  equations  of  motion  are  obtained  by  considering  the  equilibrium 
of  forces  and  moments  in  the  various  degrees  of  freedom.  The  earth  is 
assumed  to  be  an  inertial  frame  and  the  atmosphere  is  assumed  to  be 
fixed  w.r.t.  the  earth.  The  aircraft  is  assumed  to  have  a  constant  mass. 

Only  rectilinear  motion  at  relatively  low  speed  in  the  vehicle's  plane 
of  symmetry  is  considered  and  the  effects  of  unsteady  flow  and  elastic 
deformation  are  assumed  to  be  negligible.  Under  these  assumptions  the 
longitudinal  equations  of  motion  can  be  written  in  conventional  notation  (Ref.  47) 
as  : 


u  ♦  wq  + 

g  sin  8  =  X(u,  w,  q  , 

^  rjs  * 

*  rr  ,  i  ) 

E  w 

w*  -  uq  - 

g  cos  8  =  Z  (u,  w,  q  , 

4  9 

6  -  .  i  ) 

E  w 

(5.3.  1) 

• 

q  =  M(u, 

w.  q  ,  fi  T  .  f  £  ,  iw) 

6  -  q 

These  equations  are  nonlinear  and  nonautonomous  in  general.  The  pri¬ 
mary  interest  here  will  be  on  the  d  -namics  of  VTOL  aircraft  during  transition 
from  hover  to  forward  flight.  In  a  tilt -wing  vehicle,  5  ,  f  £  ■  anc*  *w 

represent  the  control  parameters,  denoting  respectively  propeller  blade 
pitch,  pitching  moment  control  and  wing  tilt  angle.  Instead  of  dealing  with  the 
complete  nonlinear  nonautonomous  equations,  they  are  simplified  in  order  to 
allow  an  analytical  treatment  and  enable  qualitative  conclusions  to  be  drawn. 
The  equations  of  motion  are  linearized  in  the  usual  way  by  making  the 
following  assumptions.  The  motion  is  considered  about  a  steady  level  flight 
and  the  vehicle  is  fully  trimmed,  i.  e.  in  a  state  of  equilibrium,  with  all 
forces  and  moments  balanced  out.  If  the  vehicle  now  encounters  a 
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disturbance  such  that  the  resultant  motion  is  small  in  magnitude,  the 

motion  is  described  by  a  set  of  l.d.e.  with  constant  coefficients.  This  is 
done  by  expanding  the  aerodynamic  forces  and  moments  in  a  Taylor's  series 
about  the  prescribed  flight  conditions  and  retaining  only  the  first  order 
terms.  This  is  consistent  with  the  assumption  that  the  disturbed  motion 
is  small  in  magnitude. 

Thus: 

x,  v  v  cx  ax  ax 

X  +  A  X  =  X  +  -  A  U  +  7 -  A  w  +  - 

cu  aw  o  T 


oX 

ai 


w 


A  1 


W 


C  Z  ?Z  *  Z  *  Z 

Z  +  *z  =  Z  + a  u  +  a  w  +  -r —  a„  +  — : —  a  i 

5u  jw  a^_  i  a  i  w 

T  w 


(5.3.2) 


a  m 


e  M 


a  m 


M+iM  =  M+  —  Su  ♦  — 1  W+iT  +  5  ..,T  ■  8  ..rii  -  „ 


a  m  a  m  a  m 

A  +  Aiv  +r—  A  1 


w 


For  steady  level  flight  the  flight  path  angle  is  nearly  horizontal 
and  this  permits  further  simplification.  The  lowest  order  terms  are 
balanced  out,  leaving  the  equations  satisfied  by  the  perturbed  variables. 

Let  these  be  denoted  by  u,  w,  and  9  .  The  linearized  homogeneous  equations 
are  therefore  given  by: 

u'-Xu-X  w  +  gQ  =  0 
u  w  b 

w'-Z  w-Z  u-VQ  =0  (5.3.3) 

w  u 

ft  "  -  Mr  ft  -  M  u  -  M  •  w  +  M  w  =  0 
0  u  w  w 

For  a  conventional  airplane  at  cruising  flight  the  stability  derivatives 
arc  constants.  The  perturbed  transient  motion  can  be  determined  by 
solving  the  coupled  linear  equations  with  constant  coefficients.  For  a  VTOL 
vehicle  executing  a  transition,  the  flight  condition  varies  from  instant  to 
instant;  hence  the  aerodynamic  parameters  of  the  vehicle,  since  they 
depend  on  the  flight  condition,  also  vary  through  the  transition.  The 
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vehicle  is  still  assumed  to  be  continuously  trimmed  throughout.  Control 
required  to  trim  is  not  considered  in  this  analysis  and  we  shall  only  consider 
transitions  at  level  flight.  The  coefficients  of  the  linearized  equations  are 
therefore  treated  as  variable  if  the  time  history  of  the  trim  conditions  can  be 
predicted.  Further,  this  change  in  the  coefficients  is  assumed  to  arise 
primarily  from  the  change  in  flight  velocity,  although  in  general  they  depend 
on  the  wing -tilt  angle  and  power  setting  6^,. 

Qualitatively  the  following  observations  can  he  made.  At  forward 
flight  a  VTOL  vehicle  behaves  essentially  like  an  airplane  and  at  hover  like 
a  helicopter  in  regard  to  dynamic  motion.  The  forces  and  moments  produced 
by  the  propellers  and  the  wing -slipstream  interaction  largely  influence  the 
low  speed  characteristics  of  a  VTOL  vehicle.  Near  cruising  speeds  these 
effects  become  less  important.  The  stability  derivatives  have  constant 
values  corresponding  to  hover  and  forward  flight,  but  change  continuously 
from  one  to  the  other  as  the  vehicle  accelerates  until  it  attains  cruising 
velocity.  At  hover  the  characteristic  roots  consist  of  a  complex  conjugate 
pair  with  positive  real  part,  and  a  pair  df  negative  real  roots.  The  motion 
therefore  exhibits  oscillatory  instability.  In  cruising  flight  the  motion  is 
characterized  by  two  pairs  of  complex  conjugate  roots,  usually  with  negative 
real  parts.  One  of  the  modes  is  of  high  frequency  ( the  short  period  motion), 
and  the  other  is  of  low  frequency  (the  phugoid  motion).  The  transition  is, 
therefore,  from  a  helicopter -like  vehicle  to  an  airplane-like  one,  with 
accompanying  difficulties  in  the  analysis  and  control  of  the  vehicle.  For 
example,  at  hover  the  vehicle  needs  forward  stick  for  forward  velocity 
(stable  trim  gradient);  but  for  cruising  flight  the  trim  curve  in  as  shown  in 
Fig.  15,  necessitating  an  adverse  control  position  gradient  at  some  time 
during  the  transition  making  it  somewhat  difficult  to  fly. 

We  shall  now  consider  a  specific  example,  a  tilt -wing  vehicle. 

Ref.  43  contains  a  comparative  study  of  the  longitudinal  stability  derivatives 
of  three  tilt -wing  VTOL  aircraft.  The  vehicles  considered  are: 

(1)  The  VZ  2  Research  Aircraft 

(2)  Two  Propeller  Transport 

(3)  Four  Propeller  Transport 
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Based  on  experimental  data  and  taking  into  account  the  scatter  of  data, 
typical  stability  derivative  variations  through  the  transition  are  proposed. 
Variation  of  the  stability  derivatives  with  velocity  for  a  vehicle  of  this 
type  used  in  the  present  analysis  were  based  on  Ref.  4S>  Now  the  equations 
of  motion  can  be  written  as 

u'  +  a;u  +  aaw+g$  =  0 

w'  +  b1w+bau-V0  =  0  (5.3.4) 

0"  +C10,  +  CaU  +  C3W  +  C4W  =  O 

where  the  coefficients  correspond  to  those  in  (5.3.3).  The  functional 
dependence  of  the  coefficients  on  velocity  is  given  in  Table  III.  Different 
values  for  the  coefficients  of  the  equations  of  motion  are  possible  due  to  a 
different  choice  of  the  stability  derivative  variation.  The  values  used  in 
this  analysis  are  not  meant  to  be  representative  of  an  optimal  flight  vehicle, 
but  are,  rather,  typical  values  based  on  existing  aircraft. 

The  wing -tilt  angle  i^  is  in  control  of  the  pilot  so  that  any  variation 
of  i^(t )  through  the  transition  can  be  programmed.  The  dependence  of 
trim  velocity  V  on  wing  angle  is  assumed  to  be  linear  and  hence  V(t)  (Fig.  17) 
can  be  chosen  conveniently.  The  stability  derivatives  are  now  expressed  as 
functions  of  t  and  this  leads  to  a  set  of  time  -variable  coupled  linear 
differential  equations.  In  the  analysis  that  follows  the  stability  der  ivatives 
3^  and  ,  being  respectively  denoted  by  -ag  and  -c3  ,  are  neglected 
since  the  contribution  of  these  terms  to  the  dynamics  of  the  vehicle  is 
considered  to  be  small. 

5.4  Two  Degree  of  Freedom  Case 

Near  hover,  the  two  degree -of -freedom  approximation  is  employed, 
in  which  the  vertical  or  plunging  motion  is  suppressed.  At  hover  (in  which 
state  the  vehicle  can  remain  indefinitely),  the  damped  vertical  or  plunging 
mode  is  completely  decoupled  and  has  little  effect  on  the  other  two  modes. 

The  system  (5.3.4)  can  be  represented  by  the  following  set  of 
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equations  at  hover  and  low  velocities ,  since  the  terms  aa  ,  ba  ,  c4  are 
very  small.  (Ref.  27): 


u1  +  aa  u  +  gfl  =  0 


(a) 


w'  +  bj  w  =  0 


(b)  (5.4.1) 


0  "  +  C;  P  '  +  Ca  U  =  0 


(c) 


Th*'  w  mode  is  completely  decoupled  leaving  the  u  and  0  equations  still 
coupled.  On  decoupling  these  by  cross -differentiation  we  arrive  at  the 
following  equations  for  u  and  0  . 


u1"  +  (aj  +  Cj  )u"  +  (aj^Cj  +  2 a.^  1  )u*  +  (ax  11  +  a!x  c*  -  gca  )u  =  0  (5.4.2) 
0’"  +  (a:  +Cl  -  S2-I.  )Q  "  +  (ax  Cl  +  clt  -Cl  ^)0’  -  gca  0  =0  (5.4.3) 

ca  ca 

On  substituting  the  quantities  from  TablelHthe  equations  become 

( 1  +  0- 1 1  )u "'  +  (0.3  +-  0.  081  t)u"  +  (0.024  0.0l22t)u'  +  0.48u  =0  (5.4.4) 


(10  +  t)a  0  +  (4  +  0.  81 1 ) ( 1 0  +  t)0  "  +  (8.1  +1 . 83 1  +  0.  122  ta  )0 


a  i 


+  4.8  (10  +  t)0  =  0 


(5.4.5) 


.  .  ,0.1+0.071,  _ 

w'  +  (  i  +  0. 1 1  W 


(5.4.6) 


Equation  (5.4.6)  can  be  readily  solved  to  give 

f  ,  0.1  +  0.07t  ,  , 

w(t|  =C  exp).  J< - 1+  0-  1  t  )  dt  > 


(5.4.7) 


The  equations  for  u  and  0,  viz.  (5.4.4)  and  (5.4.5)  are  solved  approximately 
using  the  formula  (3.4.21  c).  In  orajr  to  have  a  comparison  with  the  exact 


solution,  the  equations  were  integrated  using  the  digital  computer  for  the 
conditions  given  below. 


Initial  Condition 

’-(0) 

u'(0) 

u"(o; 

6(0) 

8'(0) 

8  "(0) 

1 

0 

0 

1 

2 

0 

1 

0 

3 

1 

0 

0 

The  linear  combination  of  the  linearly  independent  approximate  solutions 
was  plotted  in  each  case  corresponding  to  the  above  conditions.  The  locus 
of  the  roots  and  the  approximations  were  plotted  in  each  case,  together  with 
the  exact  solution. 

This  enables  one  to  retain  to  some  extent  the  familiar  ideas  of 
the  analysis  of  constant  coefficient  1. d.e.  It  is  noticed  that  the  root  loci 
corresponding  to  the  variables  u  and  8  are  different.  This  is  because,  as 
a  result  of  decoupling,  the  coefficients  of  the  equations  for  u  and  0  contain 
additional  terms  involving  the  time  derivatives  of  the  stability  derivatives. 

A  constant  coefficient  analysis  of  these  decoupled  equations  obtains  the 
"  frozen  "  approximations  to  the  solutions  of  the  variable  equations. 
Further,  in  the  light  of  multiple  time  scale^,  the  fast  time  scale  shows 
up  as  a  quadrature  over  the  root  variation  and  describes  the  frequency  of 
the  rapidly  varying  motion.  The  results  are  shown  in  Fig.  18-27. 

The  frozen  amplitude  approximation  u  (T;  )Jt  ,  which  varies  on  the 
fast  time  scale,  seems  to  represent  the  frequency  of  the  exact  solution  quite 
well,  but  the  amplitude  suffers  from  errors.  The  present  approximation 
seems  to  represent  the  true  function  well  in  both  amplitude  and  frequency. 

The  agreement  is  seen  even  with  the  other  set  of  initial  conditions.  From  the 
above  figures  it  is  seen  that  the  approximation  is  insensitive  to  initial 
conditions. 


5.5  The  Three  Degree-of-Freedom  Case 

The  complete  linearized  equations  of  motion  coupling  three  variables 
u,  w  and  8  are  gv'en  by  (5.3.4).  Again  treating  a *  and  ca  as  being 
unimportant,  .  quations  are  written  as: 

u'  +  a*  u  +  g8  =  0 

w'  +  bj  w  +  ba  u  t-VS  1  =  0  (5.  5. 1 ) 


0  "  +  Cj  0  1  +  ca  u  +  c4  w  =  0 

In  order  that  the  time  scales  approximation  be  applied,  equations  descri¬ 
bing  each  dependent  variable  must  be  obtained.  Unlike  the  constant 
coefficient  case,  the  equations  will  not  all  be  the  same,  as  emphasized  by 
Curtiss,  so  that  the  time  histories  of  u,  w  and  0  will  be  different.  But 
decoupling  the  equations  is  itself  an  involved  task;  a  schematic  of  the 
procedure  can  be  seen  in  Fig.  ]6  .  The  equation  for  u  is  given  by 


u""  +  (ax  +  bx  +  cx  -  — )u"' 

C4 

t  c 1 

+  (ax  bx  +  1^  cx  +  cx  aj  +  Vc4  +  cx  +  3^  -  (a1  +  cx  )  )u" 

c  1  ' 

+(ax  bj  Cj  +  Vc4  ax  -  gca  — —  ax  cx  +  ^  c'j  +  2ax  (bx  +cx  -  -*-)+  3^"^' 

c4  c4 

/  / 

+(g(ba  c4  ~bi  c8  +  ^—  ca  -ca  ')  +  (bj  cx  +  Vc4  +  cx  '  -  Jj4-  cx  ) 

c4  c4 


+  ax  "  (bx  +  cx  -  —  )  )  u  =0 

C4 


(5.5.2) 


When  c4  =0  it  is  seen  that  this  reduced  to  (5.2.2)(with  ax  =  constant;) 
and  when  all  the  coefficients  are  constant  the  equation  becomes 
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u""  +  (aj  +  bj  +  Cj  )  u‘"  •»  (aj  ^  +  ^01  +  ax  +  Vc4  )  u" 


+  (a1  b.  Cj  +  Vc4  aj,  -  gca  )u*  +  g(ba  c4  -  bx  ca  )u  =  0  (5.5.3) 

as  can  be  verified  directly  from  (5.5.1). 

Now  substitution  o i  thf  assumed  stability  derivative  variation  from 
Table  III  leads, after  a  considerable  amount  of  algebra,  to  the  following 
equation  for  u. 

p*  u^4  ^  +  pa  u^  3  ^  +  p»3  u<3)  +  ft  u'  +  pQU  =  0  (a) 

where 


p4  =  (10  +  16  t)(l  0  +t)3  (-0.2  +  0.  01 75  t)t  (b) 

p3  =  (2-  1.35t  -  0.232ta  +  0.02643  t3  )(10  +  16 1){10  +  t)a  (c) 

(5.5.4; 

Pa  =  (60+13.  7t  -  10. 4875 13  +4.  19t3  -0.  9988t4  +  0.  058t?)(10+ l6t)  (d) 

Px  =(4-100.  78t  -  I2.523t3  +  1.6994t3  -0.  I23t4 +0.  0107t5)(10+l6t)  (e) 

Po=3. 22(10  +t)(30. 3  +  23.  07t  -  15.  7ta  -1. 4573t3  +0.  0034t4  +0.  0095t5)  (f) 

This  equation  has  regular  singular  points  at  t=0  and  U*11.42.  Near  t=0, 
therefore,  the  fourth  order  equation  is  approximated  by  the  third  order 
system  which  has  only  two  degrees  of  freedom,  as  is  usual  in  engineering 
analysis. 

The  other  singular  point  occurs  in  a  region  in  which  two  of  the 
characteristic  roots  coalesce  and  hence  the  approximation  via  time  scales 
fails  in  such  a  region.  We  may  recall  that  this  corresponds  to  a  transition 
point,  in  the  vicinity  of  which  the  short  period  mode  changes  from  monotonic 
solutions  to  oscillatory  ones.  Some  problems  associated  with  this  are 
discussed  in  the  next  chapter. 
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We  now  examine  the  solution  to  the  fourth  order  equation  characterized 
by  various  initial  conditions.  Figures  are  drawn  to  show  the  comparison  oi 
the  time  scales  approximations  with  the  exact  solution  obtained  by  numerical 
integration.  Different  variations  of  the  stability  derivative  with  flight 
velocity  were  studied  in  relation  to  vehicle  dynamics.  Figures  18  -37 
show  the  difference  in  the  root  loci  and  the  corresponding  responses.  It  is 
seen  that  the  nature  of  the  aircraft  motion  is  qualitatively  the  same,  though 
the  root  loci  are  quite  different.  A  slight  resonance  seems  to  occur  near 
the  transition  point  on  the  real  axis.  The  agreement  of  the  approximation 
with  the  exact  solution  is  seen  to  be  good  in  each  case,  near  both  hover  and 
forward  flight,  and  thus  provides  a  uniform  description  of  the  motion  of  the 
aircraft.  The  approximation  suffers  from  errors  in  the  vicinity  of  the 
transition  point  as  may  be  expected,  but  the  qualitative  nature  of  the  solution 
is  preserved.  In  general  the  ability  of  the  approximation  to  progress  through 
the  transition  point  depends  upon  the  choice  of  the  initial  conditions; 
use  of  the  digital  computer  on  the  other  hand,  for  the  approximation,  is 
likely  to  preclude  any  difficulty  with  the  transition  point. 

Also  appended  are  the  figures  depicting  the  approximations  obtained 
by  "freezing"  the  coefficients.  A  comparison  of  these  with  exact  solutions  is 
made  for  the  third  and  fourth  order  systems  (Fig.  £0,  21,  36,  37).  It  is  seen 
that  the  "frozen"  scheme  of  approximation  is  good  only  for  short 
times,  and  the  error  becomes  large  in  less  than  a  cycle  of  the  oscillation. 

The  approximating  function  grows  without  bound  and  does  not  represent  the 
true  nature  of  the  solution  anywhere  after  the  first  cycle. 

5.6  Summary  of  the  Chapter 

Examples  with  analytically  known  solutions  are  discussed  first. 

They  consist  of  several  second  and  third  order  equations  and  a  special  n**1 
order  equation.  Correct  asymptotic  behavior  of  the  solutions  is  obtained 
using  the  theory  developed  in  this  thesis. 

Transition  dynamics  of  VTOL  aircraft  are  studied  in  the  longitudinal  mode. 
The  two  degree -of-freedom  (hover)  approximation  and  the  three  degree -of - 
freedom  case  are  studied.  The  time  scales  approximations  are  compared 
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with  numerically  obtained  solutions  for  various  initial  conditions.  Errors 

are  found  to  be  less  than  10%  for  the  hover  approximation.  For  the  three 

degree -of -freedom  case,  a  uniform  qualitative  description  of  the  vehicle 

motion  is  obtained,  with  good  accuracy  except  at  the  transition  point. 

Different  variations  of  M  (V)  were  studied. 

w 
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CHAPTER  VI 


CRITIQUE  AND  EXTENSION 

6.  1  Extension  of  the  Method 

We  shall  first  outline  the  advantages  of  the  method  over  earlier 
methods  and  then  discuss  possible  extensions  regarding  open  problems. 

At  the  outset  one  may  consider  the  novel  point  of  view  adopted  in 
the  approximation.  The  method  of  extension  as  related  to  multiple  time 
scales  has  been  used  to  obtain  the  approximations,  both  for  short  and  long 
times.  This  is  done  systematically,  generating  suitable  clocks  on  which  the 
phenomenon  is  observed  for  both  short  and  long  times.  Criteria  of  validity 
are  presented  for  each  case.  For  the  short  time  approximation  this  shows  the 
breakdown  of  the  approximation  and  for  the  long  time  analysis,  upper  bounds 
for  the  error  are  obtained.  The  formula  (3.4.21)  for  the  general  case 
enables  one  to  write  the  approximation  for  a  given  equation  by  inspection, 
and  for  each  mode  separately.  Furthermore  this  affords  a  uniform  description 
of  the  phenomenon,  in  a  region  free  of  transition  points.  Also  since  the  result 
is  obtained  analytically,  it  is  useful  for  further  study  and  investigation  of 
related  problems.  In  this  connection,  one  may  mention  the  problem  of 
obtaining  approximations  to  the  solutions  of  a  system  of  coupled  equations 
without  recourse  to  decoupling  first.  Simple  "extensions"  of  coupled  equations 
seem  to  recover  the  "frequency"  of  the  solution  but  not  the  slow  amplitude 
modulation.  However,  different  schemes  of  "extension"  may  lead  to  better 
results,  and  thus  help  to  simplify  the  analysis. 

Further  work,  for  example,  may  lead  to  the  study  of  forced 
responses  of  time -variable  systems.  Consider  the  equation: 

£(y)  =  y^n^  +  u)n_1(t)  y^n  ^+...  +  uuo(t)  =  f(t)  (6.1.1) 

The  particular  solution  can  be  written  as 

t 

y(t)  =  J  h(t,s)f(s)ds  (6.1.2) 


126 


where  h(t,  s)  is  the  Green's  function  or  tht  time -variable  impulse  -esponse. 
h(t,  s)  is  expressed  in  terms  of  the  independent  solution*:  of  the  homogeneous 
equation  £(y)  =  0.  Approximations  to  h(t,  s)  can  be  obtained  using  the  theory 
develjped  in  this  thesis  and  this  can  be  used  to  study  approximations  to  mrced 
response.  The  choice  of  forcing  functions  is  dictated  by  "resonant"  and 
"non-resonant"  cases;  in  this  connection  one  may  refer  to  the  recent  work  by 
Feshchenko  et  al  (Ref.  24)  for  an  asymptotic  theory  of  forced  linear  systems. 

Another  aspect  of  the  approximation  scheme  becomes  apparent  as 
follows.  In  order  to  apply  the  formula  (3.4.21)  one  needs  to  know  the  roots 
of  the  characteristic  equations  as  functions  of  t.  For  systems  up  to  fourth 
order,  a  closed  form  of  expression  is  available  for  the  roots,  though  it  is  not 
simple  for  the  third  and  fourth  order  equations.  For  higher  order  equations 
in  general,  no  such  results  exist.  One  may,  however,  consider  the  approximation 
for  the  roots  developed  by  the  author  (Ref.  53)  as  a  Taylor's  series 
starting  at  the  instant  t  =  0.  This  technique  is  of  necessity  limited  to  the 
region  of  validity  of  the  root  approximation,  which  has  to  be  precisely 
formulated.  Nevertheless  for  smooth  variation  of  the  roots  the  result  can  be 
used  for  small  t.  With  regard  to  the  VTOL  example,  the  characteristic 
equation  is  given  by: 


F(x,t)  3  x4  +  u)3  (t)x3  +  u)a  (t)xa  +  u) i  (t)x  +  uuQ(t )  =  0  (6.1.3) 


If  t=0  represents  hover 

x(t)  =  x(0)  +  x  (0)t  +  (0)t3  +  ... 

where : 


x(t)  = 


ill 

'If' 


(6.1.4) 


Higher  derivatives  can  be  similarly  calculated.  For  example  with  reference  to 
the  VTOL  transition  problems,  two  expansions  for  the  root  can  be  made,  one 
near  hover  and  the  other  near  forward  flight  condition.  Substitution  in  the 
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formula  yields  the  approximations  valid  near  the  two  flight  conditions 
respectively. 

Parameter  sensitivity  can  be  studied  in  the  manner  discussed  in 
Ref.  53.  The  change  in  a  characteristic  root  caused  by  changing  a  particular 
parameter  can  be  computed  as  a  function  of  t.  This  change  in  the  root  is 
reflected  as  a  change  in  the  dynamics  of  the  system. 


6.2  Transition  Point  Analysis 


The  next  question  to  consider  is  the  breakdown  of  the  approximation 
(3.4.21)  in  a  region  containing  multiple  characteristic  roots.  These  points 
are  known  as  turning  points  or  transition  points.  The  problem  of  obtaining 
suitable  approximations  valid  near  such  points  has  been  ab  initio  difficult  to 
handle.  There  is  an  extensive  mathematical  literature  on  this  subject, 
which  has  been  studied  by,  among  others,  Langer  (Ref.  54),  Wasow  (Ref.  55), 
Erdelyi  (Ref.  44),  etc.  This  section  presents  a  brief  sketch  of  the  basic 
ideas  and  some  preliminary  new  results.  The  present  objective  is  mainly 
to  identify  and  outline  the  problem  areas  and  emphasize  the  need  for  further 
work  leading  to  a  more  complete  theory. 

The  simplest  equation  exhibiting  a  transition  point  is  the  Airy 


equation: 


y"  +  ty  =  0 


(6.2.1) 


For  positive  and  negative  values  of  t  the  nature  of  the  solutions  is  quite 
different,  being  oscillatory  or  monotonic  (as  used  by  Erdelyi,  in  the  sense  of 
having  at  most  one  zero)  according  as  t  is  positive  or  negative.  Thus  t=0 
is  called  a  transition  point,  to  describe  the  transition  in  the  nature  of  the 
solution  on  either  side  of  t=0.  This  can  also  be  seen  by  observing  the 
characteristic  roots  as  t  goes  through  zero.  For  t  <  0,  the  roots  are 
x  =  ±A  and  for  t  >  0,  they  are  x  =  ±iA }  and  the  two  roots  coalesce  for  t  =  0. 

For  a  more  general  equation 

y"  +  eauo(t)y  =  0  (6.2.2) 

the  asymptotic  approximations  for  large  e  ('  LG  solutions)  are  given  by 
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for  u>  (t)  <0.  Cj  ,  Cj  ,  c3  ,  c4  are  constants. 

These  approximations  are  valid  when  u'(t)  does  not  vanish.  Clearly 

when  u)  (t)  =  0  neither  of  these  forms  is  valid  and  transition  occurs  from  one 

type  of  behavior  to  the  other.  Two  problems  are  seen  to  emerge.  One  to 

find  the  connection  between  the  constants  ct  ,  ca  ,  and  c3  ,  c4  to  represent 

asymptotically  the  same  solution  of  (6.2.2)  for  both  positive  and  negative 

values  of  t;  and  the  other  to  determine  the  asymptotic  form  of  the  solution 

of  (6.2.  2)  near  the  transition  point.  It  must  be  noted  that  a  transition  point 

can  occur  also  if  the  coefficient  u'(t)  is  singular  at  a  point  t  ,  on  either  side 

o 

of  which  u)(t)  has  opposite  signs.  Coalescing  of  the  characteristic  roots 
therefore  generally  determines  the  transition  point.  The  approximations 
each  valid  on  either  side  of  the  transition  point  break  down  near  the  point  in 
question  and  a  different  form  of  approximation  is  required.  As  Langer  points 
out  (Ref.  56)  this  can  be  observed  even  in  the  case  of  an  l.d.e.  with  constant 
coefficients;  for  the  case  of  multiple  characteristic  roots  a  different  form  of  the 
solution  must  be  used. 

Two  methods  have  been  used  to  obtain  the  connection  formulae. 

The  one  used  by  Jeffreys  replaces  uj(t)  by  a  linear  function  (  t-t^)  sufficiently 
near  t  and  integrates  the  resulting  Airy  equation  in  terms  of  Bessel  functions 
of  order  ±1/3  ,  with  known  asymptotic  behaviors.  Comparing  these  with 
(6.2.3)  and  (6.2.4)  above  one  obtains  the  connection  formulae.  The  other 
method  used  by  Zwaan  (Ref.  57)  consists  in  integrating  (6.2.2)  on  a  complex 
plane  along  the  real  axis  up  to  the  point  t  on  either  side,  but  making  an 
excursion  into  the  complex  plane  along  a  semicircle  to  connect  the  two  sides. 
This  avoids  the  transition  point  altogether  and  obtains  the  same  connection 
formulae  as  before.  As  discussed  by  Erdelyi  (Ref.  40)  both  methods  can  be 
extended  to  cases  where  o>  (t)  has  a  zero  of  an  arbitrary  order. 

The  second  problem  is  one  of  more  mathematical  interest  and  it  is  the 

determination  of  the  asymptotic  forms  of  the  solution  near  the  transition  point. 
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The  works  of  Langer,  Olver,  Cherry,  and  Erdelyi,  referenced  in  Ref.  55, 
are  prominent  in  this  respect.  No  simple  elementary  function  seems 
adequate  in  representing  the  transition  from  oscillatory  to  monotonic  behavior 
and  higher  transcendental  functions  seem  to  be  needed.  The  works  of  the 
above  authors  deal  with  the  uniform  asymptotic  representations  in  terms  of 
Bessel  and  Airy  functions,  etc.  These  are  not  limited  to  the  vicinity  of  the 
transition  points  alone  but  are  valid  uniformly  in  the  domain  of  interest. 

With  reference  to  the  hover,  forward -flight  transition  of  a  VTOL 
aircraft,  the  characteristic  equation  is  seen  to  have  a  double  root  at  one 
instant  during  such  a  transition.  Equation  (5. £.4)  shows  that  this  occurs  for 
a  value  of  t  near  11.4.  This  shows  that  the  roots  which  eventually  correspond 
to  the  short  period  mode  change  from  real  ones  to  a  pair  of  complex  conjugates 
for  a  t  in  the  neighborhood  of  t  =  11.4.  The  solution,  therefore,  changes 
from  monotonic  subsidence  to  oscillatory  subsidence,  with  the  accompanying 
breakdown  of  the  approximate  solution.  In  order  to  use  an  approximation  from 
hover  to  forward  flight,  the  Stokes  phenomenon  (se  Chapter  I)  must  be 
investigated.  For  the  aircraft  problem,  the  precise  phase  of  the  solution  is 
relatively  unimportant.  Great  precision  in  the  knowledge  of  the  frequency 
and  damping  of  the  motion  is  seldom  required. 

The  amplitude  variation  as  given  by  the  approximation  grows 
without  bound  as  the  transition  point  is  approached.  Hence  proper  connection 
may  be  necessary  in  order  to  obtain  usable  solutions,  as  the  amplitude 
information  may  be  required  for  feedback  control  purposes. 

6.3  Shifting  of  the  Transition  Point 

We  shall  first  consider  the  second  order  1. d.e.  and  show  that  the 
approximations  derived  in  this  thesis  can  be  used  to  advantage  in  dealing 
with  the  transition  point  problem. 

In  approximating  the  solutions  to  the  noncanonical  equation 

y"  +  euJi  y'  +  c8  u>Qy  =  0  (6.3.  1) 

there  are  two  methods  of  approach.  One  is  to  transform  the  noncanonical 
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equation  into  the  canonical  form  and  then  obtain  the  LG  approximation. 

The  other  is  to  treat  the  noncanonical  equation  directly.  It  is  seen  that 
under  certain  conditions,  the  two  approximations  fail  at  different  points. 
Considering  the  first  case,  we  convert  (6.  3.  1)  into  the  canonical  form  by 
means  of  the  transformation 

(6.3.2) 

(6.3.3) 


!•  6*  i 


Z”  -.■{  iU- .  JllSo  +4J_)Z  =  0 
4  2c 

u:  3  4o)  a' 

For  — * — ^ - 0  +  >  0  ,  the  LG  approximation  yields: 

y(t)  =  Ax  (ujx  3  -  4wq+  exp(  -  |-|u)1dt  +  ^-j(u;x3-  4u>o+^pr  dt) 

+  Aa  (cux  1  -  4(u  +  )  ^4exp(  -  =-  1  u>x  dt  -  ~  f (u:  x  3  -  4u>  +  )^dt)  (6.  3. 4) 

Oc  C  J  C  J  o  c 


On  the  other  hand,  applying  the  time  scales  formula  (3.4.21)  to 
equation  (6.3.  1)  directly,  another  approximation  is  obtained  as 

y  (t)  =  Bx(<u  x  3  -  4u)o)"l/4  y(t)exp(-^-J  dt  +  J(u>x  3  -  4^)^  dt) 

+  Ba  (u>x  3  -  4uuo)  ^  y(t)exp(-|j’  0'4  dt  -  J(u  x  3  -  4u)q)^  dt  )  (a) 

(6.3.5) 

where  y(t)  =  exp(  -g*-  —-p/a )  (b) 

1  o 

Now  y"  is  unbounded  when 
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and  y  is  unbounded  when: 

u>,  2  -  4(u  =0 

o 

Therefore  if  ij  ^  0  and  (mj  3  -  4u>o )  and  u):  (t)  do  not  vanish  simultaneously, 

we  have  two  approximations,  each  bounded  when  the  other  is  not.  The 
transition  pointtfor  (6.3. 1)  and  (6.  3.  3)  are  given  by  a  value  of  t  for  which 

U)  ,  3  -  4(1)  =0 

*■  o 

and 

• 

-  2(u, 

(d  .  3  -  4(d  +  — *  =  0 

1  o  e 

respectively.  The  occurence  of  the  singularity  at  the  transition  point  is 

caused  only  by  the  use  of  the  approximation  and  is  not  intrinsic  to  the  original 

differential  equation  which  may  have  solutions  well  behaved  throughout  the 

domain  of  interest.  Thus  we  have  obtained  two  approximations  which  have 

different  transition  points.  Therefore,  in  effect,  the  transition  point  has 

been  shifted  from  t  to  t  1  where 

o  o 

u> !  3  (t  )  -  4(1)  (t  )  =  0 
o  o  o 

and: 

•  § 

u>xa<t  ’)  -4a,  (t  ')  +  11>^A  =  0 
*  o  o  o  C 

If  €  i*»  very  large  the  shift  is  small. 

Since  the  bounds  on  the  error  are  known  in  each  case,  the  idea  can 
be  used  to  shift  the  transition  point  by  a  desired  amount. 

An  alternative  view  is  as  follows.  Consider  the  canonical  second 
order  1.  d.  e. : 

Z"  -  e3fl( t)Z  =  0  (6.3.6) 
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L-  G  theory  yields,  for  Ci(t)  >  0  (nonoscillatory  case): 


Z(t)  r  Aa  exp(  ^-j/l^dt)  +  Aa  exp(  -  ~  J.fl^dt) 

Now  consider  the  transformation  (6.  3. 2).  The  equation  (6.3.6)  it. 
transformed  into  (6.3.1)  with: 

O  1  ,  3  .  2(|), 

JZ.=  — (o)l3-4o)  +  — J  ); 

4  o  c 

D  =  a)X  2  -  4u>o  =  4  (Cl-  p) 


(6.3.7) 

(6.3.8) 

(6.3.9) 


Using  (6.3.5)  y  (t)  can  be  written  as: 

y  (t)  =  Bl  (-^-p)  ^  y(t)exp{-  j  J  uj(  dt  +  J  [  4(12 -p)l  ^  dt  ) 
+  Ba  (/l.p)-1/4v(t)exp(- dt  -  |J[4(.ft-p)]*dt  ) 


(a)  (6.3.10) 


where  now  y(t)  =  exp  (  ) 


Using  (6.3.2): 

Z(t)  =  B:  (-^-p)  y  (t)exp(cJ(U-p)^dt)  +  Bs  (rt-p)  ^y(t)exp(-cj(^-p)^dt) 


(6.3.  11) 


By  Olver's  theorem  and  its  extension  (Chapter  IV)  the  errors  can 
be  computed  as  follows  (standard  L  G  theory): 


Z.(t)  =  Z.(t)  (  1  +  E  ) 

i 


i  =  1.  2 


where 


Sil  s  e*p  <  h  1  jp.  s?  ‘h’7')  |  <“> 


(6.3. 12) 


|e  s  exp 


(  2€  [  IT?4,  ^  ( -a"l/4 )  I dt 


)  - 1 
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Also  (for  the  alternative  approximation,  when  y  =  constant); 


Z.(t)  (t)(l+E.  )  i  =  1,  2  where 

1  l  '  i 

lEjl*  exP<  Jl  J  I  "  2D  5?  (  D  "^4  >  dt  >  “  1 

!’£3|^  exp  (  27  1  J  "  2°  dP  ^  D-1/4  )  dt)  -  1 

We  recall  that  the  condition  for  y  to  be  a  constant  to  this  order  is  that 

«>!  =0(  — ).  In  this  case  the  shift  in  transition  point  is  0(  —  3  ) . 

C  G 

On  substituting  from  (6.3.9)  and  simplifying,  we  obtain  (for  constant  y): 

z  (t)  =  Z  (t)  (  l  +E.  ) 

XI  1 

1=,!  4  l  1 i?  ( (n-p)'lA  >/d*>  -1 

& 

’£,1*  «p >  ldt  >  - 1 

Similarly  for  the  oscillatory  case  consider  the  equation: 


(6.3.14) 


Z"  +ca(2z  =  0;  £l(t)  >  0  in  (a.b) 


(6.3.  15) 


Using  (6.  3.8),  4-^=  -  (  ou x  2  -  4uuq  +  ) 


i.  e. :  4(il+  p)  *  D,  where  D,  =  -(yj,  3  -  4uu  ) 

*  *  a  O 


(6.3. 16) 


As  before,  denoting  the  L  G  and  the  time  scales  approximations  by  Z  and  2^ 
respectively,  y/e  can  write 

Z.  =  Z.  +  E.  (a) 


l  l  l 


t  a 

;re  |  ^  a  j  i  exp(  ±  —  J  ~qV*  ^  )  dt  )  ”  1  0>)  (6. 3.  17) 


Z.^.+E.^+p)-^ 


(a)  (6.3.18) 
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where  ,  for  constant  y  : 

|  Ei  a  I  *  exp{  “ 

Thus  in  a  given  equation 

y"  -  f  siTy  =  0 

suppose  that  il(t)  has  a  zero  at  t  and  fT(t)  >  0  for  a£t<0.  in  order  to  use 

o 

(6.3.14)  a  function  p  is  chosen  such  that  p  >  0  for  a^t<t  1  and 

o 

il(t  ')  -  p(t  ')  =  0.  We  require  that  t  1  >t  .  Hence  p(t  )  <0. 
o  o  o  o  o 

Similarly,  given  the  equation 
y"  +  eajQ.y  =  0  , 

let  -O-(t)  ^.0  for  t \tQ  ;  also  let  tQ^  ti  b  be  the  region  of  interest. 

A-,  t  is  a  transition  point  we  choose  a  p  sufch  that  p(tQ)  >  0  and  we  can  use 
(6.3.  18)  to  estimate  the  errors  of  approximation.  A  schematic  is  illustrated 
in  Fig.  39  for  -O-(t)  =  t. 

6.4  Choice  of  the  p  function. 

The  choice  of  the  shifting  function  p  is  governed  by  the  following 

considerations.  In  a  finite  domain  problem  a  proper  p  function  must  ensure 

that  a  transition  point  does  not  occur  in  the  domain  of  interest.  Thus  if 

i  1  (t  )  =  0  and  a  <  t  <  b  then  0-+  p  must  not  vanish  in  a  £  t  £  b.  If  time 
o  o 

is  the  independent  variable  as  is  usually  the  case  in  dynamical  systems,  the 
range  of  t  is  the  semi- infinite  domain  0  s  t  s  ®  .  The  transition  point  is 
then  moved  to  the  negative  t  axis.  The  function  p  must  be  chosen  such 
that -ft- +  p  is  essentially  the  same  as  fleverywhere  except  near  the  transition 
point  t^  where  p  is  non  zero,  p  is  therefore  a  peaked  function  near  t 
and  sharply  decays  to  zero  on  either  side  of  it.  p  need  not  be  symmetrical 
about  the  transition  point;  in  fact  an  asymmetrical  p  may  prove  more  useful, 
for,  on  one  side  it  must  shift  the  turning  point  while  on  the  other  it  must 
decay  sharply  to  havefi+  p  k*  a . 

Even  in  cases  free  of  transition  points  the  p  function  can  be  used 
to  advantage  in  reducing  the  error  of  approximation.  For  example, 
considering  a  finite  domain  problem,  the  error  of  the  LG  approximation 
may  be  more  than  a  specified  value.  Since  the  errors  of  approximation  are  known 


c2  P 


UW  -  a? 1  (a+p)'1/4  > 


dt 


)  -1 


(b) 
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through  Olver's  results,  a  p  function  may  be  chosen  to  keep  the  errors 
within  specified  limits.  This  might  be  approached  as  a  variational  problem. 
Consider  the  equation 

y"  +  e3  A(t)y  =  0  (6.4.1) 

in  the  interval : 

a  £  t  £  b 

We  wish  to  choose  a  p(e  ,t)  such  that  the  error  of  approximation  is  minimum. 
The  approximation  itself  is  in  terms  of  elementary  operations  and  functions 
and  reduces  to  the  exact  solutions  only  in  special  cases.  From  (6.3.  18)  it 
is  seen  that  the  error  is  minimum  when 


_i!_r 

(n+P)ya 


-  (n+.p)‘1/4 


dt‘ 


(  (/2+pf^4 


is  minimum.  Thus 

‘max1!  ‘(P.r.P.t)dt 
a 


where 

f(p.p.p.t) 


e3  P  _5_  (A  +  p)3  .  I  (A  +  V  ) 

(O+pF3  16  (n  +P)5/2  4  (TfTpF5 


(6.4.2) 


(6.4.3) 


(6.4.4) 


Application  of  the  Euler -Poisson  theory  (Ref.  5  8)  to  the  above  leads  to  an 
equation,  the  solution  of  which  would  yield  the  p  function  for  minimum  error 
of  approximation  in  the  interval  (a,b). 

It  is  felt  that  a  similar  approach  may  prove  useful  for  higher  order 
equations  also. 
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SUMMARY  AND  CONCLUSIONS 

The  main  results  of  this  dissertation  are  summarized  below. 

Approximations  are  obtained  to  the  solutions  of  linear  differential 
equations  by  suitably  extending  the  domain  of  the  independent  variable 
using  multiple  time  scales.  For  a  large  class  of  problems,  linear 
time  scales  are  found  to  be  inadequate  and,  therefore,  nonlinear  clocks  are 
employed,  on  which  the  solutions  are  observed.  The  clocks  depend  on  the 
coefficients  of  the  original  equations  and  are  determined  by  a  rational 
procedure.  The  Liouville -Green  (or  WKBJ)  approximation  is  obtained 
using  this  method.  For  the  noncanonical  second  order  equation  another 
approximation  is  proposed,  and  under  certain  conditions,  this  remains 
bounded  where  the  WKBJ  functions  become  unbounded. 

In  obtaining  the  approximations,  only  the  domain  of  the  independent 
variable  is  extended,  so  that  this  would  correspond  to  the  lowest  order, 
in  an  expansion  of  the  dependent  variable.  Specific  criteria  of  the  uniformity 
of  the  asymptotic  expansion  are  not  applied  per  se;  however,  they  are 
implicit  in  that  tne  "counterterms",  --i.e.  clocks--are  so  chosen  as 
to  cancel  the  nonuniform  parts  of  direct  perturbation  theory.  The  extended 
perturbation  equations,  therefore,  are  forced  to  be  homogeneous  equations. 

A  brief  discussion  of  the  criteria  of  uniformity  is  presented  in  Appendix  V. 

The  method  can  be  extended  to  obtain  higher  order  approximations  in  a 
straightforward  way. 

The  validity  of  the  approximation  scheme  in  different  intervals  is 
examined  and  criteria  of  the  failure  of  the  approximation  are  proposed. 

Error  bounds  of  the  Olver  type  are  derived  for  the  second  order  equation 
in  a  direct  way,  although  under  restrictive  conditions.  Similar  approximation 
theorems  axe  proved  for  third  order  equations. 

Applications  of  the  approximation  to  the  dynamics  of  VTOL  aircraft 
through  the  hover -forward  flight  transition  shows  good  accuracy  in 
comparison  with  solutions  obtained  by  numerical  integration.  For  the 
decoupled  equations  which  are  of  fourth  order,  one  may  expect  the  approximations 
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to  fail  at  a  time  which  corresponds  to  multiple  characteristic  roots.  The 
phugoid  mode  itself  is  not  subject  to  this  difficulty  and  may  be  isolated 
by  proper  choice  of  initial  conditions.  If  the  functions  are  computed 
using  a  digital  computer,  again,  the  difficulty  at  the  transition  point  may  be 
avoided.  Further,  different  variations  of  M^CV)  were  investigated; 
however,  no  appreciable  difference  in  the  nature  of  the  responses  was 
found  to  occur. 
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TABLE  I 


NONUNIFORMITY  IN  PERTURBATION  THEORY 

y"  +  e  U)  y  -  0 


U) 

~  dti 

o 

~  1  j  a)  dta 

Type 

of 

Nonuniformity 

1 

ta 

t3 

secular 

3  ! 

3! 

n 

n+2 

n+3 

t 

t 

t 

secular,  n>  -2 

(n+2)(n+3) 

(n+2){n+3 ) 

1 

? 

l  nt  -  1 

/,  n  t 

secular,  t-® 

singular,  t-*0 

1 

1 

? 

4nt 

t 

singular 
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TABLE  Ii 


EXTENDED  DERIVATIVES 


fV  T:  }  ;  To  =  t  ;  rx  =  ek(t) 


d  ..  a  x  *  o 
dc  ar  +ck  an 

o 


d2 

d? 


fro»  +  e(  kfn  +  £3  ^frl3* 


a3 


r_a  a1 


dt‘ 


3T 


+  e  ( *k  •— —  +  3k  i — — 

ati  at  arx 


+  3k 


;  a1 


aroaaTl 


+  ea  (3k  k 


ar; 


+  3k2 


aroaTx 


)  +  e  *  (  k3 


*T\ 


d4  .  a4  a  x  •••  a3 

dt  a t  aTx  st  art 

n  •  a  ■ 


+  6  k 


a ' 


3To  STl 


+  4k  d 


aroaTx 


•  •  *s  2  2  3 

51  p  0  /,  7*  (j  ,  0* 0 


+  c3  (  3k3  - -  2  +  4kk  -2 —  a  +  l2kk  - a 

a rl3  a t x 3  aTaTj3 


+  6  k3  y 


a  r  -  a  r-j 

A  * 


■a  ) 


+  e3  '  6k3  k  - —  3 
3Ta  3 


•  -s  4 

+  4  k3  0 


+  e‘(i‘‘  ITT1 


d 

d? 


a n  ,  .  t  n-l  -n-1  a‘ 

^p,  +  c  ()  +  ...+  c  (nk  -T 

o  o  i 


■n  n(n-l)*n.2..  3 

TT=1  +~? - k  k 


n-l 

t r- 


atx 


n 
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TABLE  HI 


STABILITY  DERIVATIVE  VARIATION  FOR  TYPICAL  VTOL  AIRCRAFT 


=  -  X  =0.2 


b,  =  -  Z  =  0. 1  +  0.  004  V 
1  w 


h  _  7  .  O'2*  V 

b3  "  "  u  '  10  +  V 


c.  =  -  M  =  0.  ]  +  0.  0034  V 

q 


ca  =  -  M  =  0.015  (-1  +  *—  ) 
u  150 


c.  =  -  M 


w 


Case 

c*  (V) 

V(t) 

1 

(-0.  02  +  0.  00025  V)  ~A 

1  D  U 

150  t 

10  +  t 

2 

°-02<-rb>a 

II 

3 

0.  005  +  0.  015  (-iJo-)3 

II 

0.005  +  0.  015  (ttr-)3 

1  D  U 


150  t 
20  +  t 


TABLE  IV 


EXAMPLES  OF  SOME  CLASSICAL  EQUATIONS 


Name 

Equation 

Asymptotic  Behavior 

1.  Bessel's  Eq.  of 
Zeroth  Order 

y"+  ^-y '  +  y  =0 

4 

t  exp(±  it);  t-*® 

2.  Confluent 

Hyper  geometric 

ty"  +  (b-t)y'  -ay  =  0 

t  a-b 

e  t  t  “*  00 

►  > 

-a 

3.  Euler's  Equation 

3 

y'"  +  p-  y  =  0 

i-i 

{ lY3  >} 

t 

■ 

4.  "Double  Airy" 

y"'  -  4ty'  -  2y  =  0 

1  } 

■ 

5.  Euler's  Equation 

(n)  ,  ,  e_.n 

y  +  (7-)  y  =  0 

Ar-  e  (  tri/ n) 
t  t  <  g  “*® 
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APPENDIX  I 

EXTENSION  OF  THE  nth  ORDER  DERIVATIVE 
With  the  two  time  scale  extension 

t  r  o.  t !  }  ;  Tq  =  t,  rx  =  e  k(t)  (1) 


The  derivative  operator  is  extended  as 


d_  t _  3 _  drT  9_  k  — 

dt  3T  dt  3TX  '  dT  E  aTx  (2) 

o  o  * 


<  Fro  *  £it  I?,1  {h  +  ‘Hr,1 


+  £  <  Hr, 


+  2k  %-r,  > 


(k3!— j)  (3) 

dTx 


Similarly 


a 

*To 


+ 


+  e(  )  +  Ca(  )  +  ...  +  ) 


(4) 


(5) 


Clearly  the  r.h.s.  contains  terms  due  to  the  binomial  expansion  of  the 
operator  and  those  due  to  successive  derivatives  of  the  clock  function. 

For  purposes  of  the  present  approximation  scheme,  only  terms  of  order 
/‘and  en  are  needed  in  addition  to  the  lowest  order  terms.  The  terms  are 
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n 


^  T  n 
dT0 


(a) 


r  +  l 

(n)  3 _  ,  r’1  /  n  \  Jn-r)  lrT>r 

'  TTl  r=l  (  r  '  k  ST0  3Tl 


(6)  (b) 


n-1 


,r,n-l  ?n  ,  n(n-l)  „*  ;n-2  lL,  S n  1  ,  , 

n(k)  5TTrTn''1  +  ~T~  |k>  *■  <c> 


n 


<k>"  h  “ 


(d) 


That  these  are  indeed  the  terms  is  proven  by  the  principle  of  mathematical 
induction  as  follows.  We  shall  prove  that  if  (6)  is  true  for  n  then  it 
is  true  for  n+1  and  show  that  it  is  indeed  true  for  one  n.  Letting  the 
derivative  operator  (2  )  act  again  on  (6),  : 


3T 


n  ,  ,  k(»)  i_  .  yf 1  /  n  x  Jn-r) 

ar:  '  l  stx  &i  r  aToraTx 


n-1  ,  _/,*  »n"l  3 


n 


n(n  -1)  .n- 2  "  a 


n-1 


+  •  •  •  +  c  <  n(k)  srsr,^  +  2  . 8T. 

* 


+  e 


n  (  (k)n 


n 


arx 


■n 


(7) 


The  various  terms  can  be  written  as: 


3  n”* 

aF^T 


(a) 


(n+1)  3  ,  , (n)  3 


3Ti 


+  k' 


dVr* 


+  k 


n+1 


&v  aTi 


n-1 


+  S  (*) 

r  =  l  r 


/  -  r+l 
(n -r+l)  a 


3  T  r  3rx 
o  1 


+  k 


(n-r)  3 


r+2 


(8) 


a'f7^-aTx 


l-i-: 


(n+1)  a y  ,  n+i  (n-r+1)  a  r  +  1 

STj  r=l  (  r  S  T  r  5  T  i 


(b) 


n 


•  ,  -s  n  n+l 

,,  ,n-l  ;<  a  ,  •  n  3 

n(k)  k  - - ^  +  kn 


;  ,n  o 


n+l 


o  T  i 


at  a r  1  u  +  n(k)  ar  arx  n 


n 


+  n(n-l)  *  n-1  j.  &_ 

2  (k)  BTi" 


/  ,  ,  „',n  5*1  *  ,  (n+l)n  *  n-1  "  a n 

=  (n+1)(k>  1PF7F7"  +  —2~m  ksn" 


(c) 


n+l 


,  _  n+l 

,n+i  a 

(k)  -,---n+~i 

0  7, 


(d) 


On  examining  (6)  and  (8)  it  is  seen  that  (8)  is  obtained  from  (6) 
by  replacing  n  by  (n+l).  Hence,  if  (6)  is  true  for  n  then  it  is  true  for 
(n+l).  It  is  easily  verified  from  (6)  and  (3)  that  it  is  true  for  n  =  2; 
thus  it  is  true  for  any  n. 
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APPENDIX  II 


CONDITION  FOR  THE  INVARIANCE  OF  THE  AMPLITUDE  FUNCTION 
W.R.T.  THE  CHARACTERISTIC  ROOTS 

Let  the  differential  equation  be  written  as  : 

y'"  +  e3  3u)j  (t)y '  -  e32uu0(t)y  =  0  (1) 

The  characteristic  equation  is  given  by: 

F  =  x3  +  3u) j  (t)x  -  2u)Q(t)  =  0  (2) 

Let  the  functions  s^  and  sa  be  defined  as  follows 

s,  3  =  a)  +  (u)!3  +  u)  3  ;  8S  3  =  w  -  (u>  3  +  u>  a  )^*  (33 

*  0  *0  o  1  o 

and,  therefore 

8i  8a  ~  ”^1 

taking  the  real  quantity.  The  roots  of  (2)  are  given  by  (Ref.  50,  p.  17): 

xx  =  (si  :  83 )  ;  xa(3  =  -  (  S*  *  8s  )  ±i  (»i  -  »a  )  (4) 

The  approximation  to  the  solutions  of  (1)  for  large  e  is  given  as 

y (t)  =  C  D(t)  exp (  e  l  x  ds)  (5) 

where 

D(t)  =(|^)"1/3  =  (xa  +  u)X  )"*/*  (6) 

0  x 

where  x.  are  the  three  roots  of  (2)  given  in  (4). 
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We  wish  to  show  that  under  certain  conditions  the  dependence  of 


D  on  t  is  invariant  w.r.t.  the  roots. 

Since  ujj,  is  a  given  function,  it  is  sufficient  to  show  that 

|  x3  |  .  i •  . 

—  I  -•  a  constant  in  some  limit. 

I  u>i  i 

Taking  xl  =  sa  +  s3  : 


x3  |  _  (Si  +  s3  J3  S! 

+  2  +  — 3 

U)  1 1  81  s3  ®3 

®i 

Thus 

1  —  -•  constant  if 

PJ 

-*  constant 

1  IDj 

>  sa  1 

1.  e. 


+  (UD.3  +U)3)^ 

4  o 

-  (,U>  1  3  +0)  3  ) 

o 


constant 


l.  e. 


2uj  3  +  uji  3  +  2u)  (  a;,  3  +in  3  ) 

o  o  *  o 

-g -  —  constant 


U)  • 


Thus  it  is  sufficient  that  liM  -•  constant  for  I  —  |  to  approach 

®0  '  11,1 

a  constant. 


Similarly  for  the  other  two  roots.  In  this  case  the  function: 


D(t)  = 


<*■  +  -« ^  7  +  i>* 


4 


=  (constant)  a>a  (t) 


(ii .  3  ( t ) 

Thus  when  — *a~/; T  =  constant,  D(t)  =  Uj  (t) 


APPENDIX  III 


CONDITION  FOR  CONSTANCY  OF  SIGN  OF  fy 

(1) 

(2) 

or  (ii)  Yi  (K)  >  y i  (€ )  o 

and  yx (C  )  ^  >  0 


In  equation 

y "  +  e  2  uu  y  =  0 

if  x  (t)  s  0,  the  characteristic  roots  are  real. 
Now  if  (i)  yx  (t)  0  ,  i  b 


then  by  corollary  to  Sturm's  theorem  y ^  does  not  change  sign  in  (a,b) 
and  is  positive. 

The  condition  for  constancy  of  sign  of  f(t)  is  examined  as  follows. 


f(t>  =  =  i4 : ] 


f(t)i0  if  -  -j-  (-)’  0 

uu  4  uu 


i  •  c  • 


4  (  z  +  za  )  \  5  za  ;  where  z  *  ■-  . 


(4) 

(5) 

(6) 


-  • 

i.  e.  4z  \  z  .  Thus  z  >  0. 

For  example,  if  w(t)  =  tn,  we  see  that  f  2  0  if  -4  £  n  <  0.  For  a 
more  general  ui(t),  conditions  must  be  similarly  established. 
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APPENDIX  IV 


PLK  METHOD  APPLIED  TO  A  SECULAR  PERTURBATION  PROBLEM 


Consider  the  equation  : 

y1  +  e  y  =  0 

The  variables  are  extended  a  la  Lighthill  as  follows. 
y(t)  <=^yo(s)  +  e  yx  (s)  +  . . . 

t  =  s  +  ett  (s)  +  ... 


(1) 


(2) 


dv 

Therefore:  r-L  +  £y  (  1  +  e  fJ  +  ...)  =  0 

ds 

(3) 

Order  by  order  the  equations  are 

^=0 

(a) 

i!L=  .  v 
ds 

(b) 

(4) 

^=-(7.  +y0t-.'l 

(c) 

and  so  on,  giving  : 


r  =  C  =  constant 
o  o 


Yi  =  -  C  s  +  Cx 
o 


Using  the  uniformity  condition: 
dtx  y, 

—  =  -  =  s  or 

ds  y 

'a 


From  (2), 
Solving  for 


.  =  s  +  e  {  |—  +  Ca  ) 
s,  we  have  a  =  -1  ± 


to  order*  e  . 
(1-c  (ca  r  -t)  )* 


(a)  (5) 
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When  the  constant  Ci  =  0,  s  =  -1  ±  (1  +  e t) 


+  e  t)1/3 


(b)  (5) 


The  exact  solution  of  (1)  is  obtained  as: 
y  =  A  exp(  -e  t) 

We  see  from  (5)  that  the  realtion  between  s  and  t  is  an  algebraic  one 
and  cannot  be  expected  to  capture  the  exponential  variation  of  y(t). 

Alternatively,  since  the  PLK  method  has  been  usually  applied  to 
singular  perturbation  problems,  one  may  be  prompted  to  convert  (1)  to  this 
form  and  then  apply  the  method. 

D  'lining  t  =  ~  »  (1)  becomes 


ll.  . 


c  y  =  0 


This  is  a  singular  perturbation  problem  as  can  be  seen  by  expanding 
y  =  yQ  +  e  yx  +  . . .  ;  whence: 

y  =  A  =  Constant 
o 


yx  =  -  -  etc. 
x 


Now  applying  the  PLK  method,  we  expand: 


y  «==>  y  (*)  +  eyi  (*)+.. . 
o 


x  —  z  +  €  Xx  (z )  +  . . . 

Substituting  (8)  into  (7),  the  resulting  equations  are  written  as  follows: 

«■  to.  „  o 


.»  to- , .  1 2x,  *to  . y , 

dz  *  dz  ’ o 


(b)  (9) 


»*  to..  -[  2Xl  ♦  ,„»  +  2*  -  v  la 
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Therefore,  y  =  C  =  constant 
o  o 

C 

Vl  =  -  ^  +  cx 


Setting  the  r.h.  s.  of  (9c)  for  uniformity: 


dxj 
^  dz 


2  Xl 


+  yi  =  0 


Integrating: 


xi 


(Ca  +  1)  z  +  C3  za 


(10) 


Now  from  (8)  :  x  =  z  +  e  Xl  ;  again,  this  being  an  algebraic  relation, 

it  cannot  be  expected  to  describe  the  exponential  behavior  of  the  exact 
solution; 

y(x)  =  C  exp  ) 

It  is  felt  that  this  simple  example  demonstrates  the  difference  between 
the  PLK  and  the  multiple  time  scales  approaches. 
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APPENDIX  V 


£  (t)  is  therefore  a  uniformly  valid  approximation  to  f(t)  in  a  region  J  if 
o 

and  only  if 

~  =  i  +  ou> 

o 

for  all  t  in  J.  The  ratio  of  f  to  fQ  approaches  unity  as  £  0.  This  form 

of  the  criterion  is  used,  for  example,  by  Erdelyi  (Ref.  40). 

2.  Difference  Criterion. 

From  the  above  expansion  for  f  : 

f  -  f  =  e  fj  +  . . . 
o 

Th’s  requires  that  the  difference  between  f  and  f  approach  zero  as  c  0. 

The  difference  criterion  can  be  misleading  when  dealing  with  very 
small  or  very  large  quantities,  as  can  be  seen  by  the  following  example. 
Consider  the  equation 

y 1  =  ( 1  +  e)y  =  0 
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with  the  condition:  y { 0 )  =  1 

The  exact  solution  is  written  as  y  =  exp  (  -(1+e  )  t).  A  perturbation 
expansion  yields,  in  lowest  order: 


y'+y=0  or  y=C  exp(-t) 
7o  ’ o  ’o  r 


The  difference 

y  -  yQ  =  exp(  -(1+e  )t)  -  exp(-t) 

approaches  zero  for  large  t  and  hence  according  to  the  difference 
criterion  one  would  conclude  y^  to  be  a  uniformly  valid  approximation  of 
y(t)  for  large  t,  to  order  £  .  This  is  severely  in  error  bece.use,  although 
the  difference  approaches  zero,  the  functions  themselves  are  vanishingly 
small. 

y  t  f 

However,  the  ratio  test  shows  that  =  exp(-et)  — 0  and  does 

yD 

not  tend  to  1 . 

For  the  equation  y'  -  (1+e  )y  =  0  ;  y ( 0)  =  1,  a  similar  analysis  shows: 


and 


y  =  exp(  (1+e  )t);  yQ  =  exp(t) 

R3  =  exp(et);  D  =  y  -  y  =  exp(  (1+e  )t)  -  exp(t) 
yo 


We  require  R-*  1  and  D”*0  as  t-*®.  The  actual  values  are  given  below. 


t 

R 

D 

0 

1 

0 

e 

e1/3  (e -1 )  ea  1.7  el//c 

£ 

1 

e1/6 

e^’le13-!).^8 

e3 

Thus  for  large  t  the  difference  error  is  much  larger  than  the  ratio  error. 
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Figure  6.  Extension  of  the  Domain 


Figure  7.  Function  Surface  in  Extended  Space 
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Figure  8.  Trajectories  for  Restriction 
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FIG.  lie  BALANCING  OF  TERMS 
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